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HYDROLOGY OF LAKE a 
| I. Morton, 1 and H. B. 


of water into Lake Ontario are discussed, Short and long term 
of flows and supplies of Lake Ontario are related to the factors causing varia 


made in relation to seasonal and areal distribution of precipitation, ‘tempera- 


_ ture and runoff. Flow records of many of the large tributaries of Lake roe ™ 
_ Ontario are analysed and the cir ested = long and short-term forecasting is 


The International ‘st. Lawrence River Board of Control is 
veloping plans for regulating the levels and outflows of Lake Ontario in on i 
_ junction with the St. Lawrence Seaway and Power projects. Thetechnical 
_ studies are being performed jointly by various governmental agencies in es 


United States and Canada. The Water Resources Branch of the Department of se 
Northern Affairs and National Resources is one of the Canadian agencies as- & 
; ‘sociated with this work. The problem : is international in all its ramifications _ 
since the International Boundary traverses not only Lake Ontario andthe St. 
Lawrence River to the vicinity of Cornwall, Ontario, but also the upper Great 
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tion in the water supplies. Detailed descriptions of the source of supply are z 
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conditions, tea water supplies utilized in testing the plans have been derived 
from historical outflow and water level changes and results beneficial to all 
interests have been achieved without the necessity for a detailed study of the 
i _ Sources of supply and the other hydrologic factors, If further refinements __ 5 
are to be made, operating experience as well mons more detailed | knowledge a of a 
eS _ The question may be asked as to what may be gained from a better under-— ae 
standing of the various factors contributing to the water supply to Lake a : 
Ontario. Firstly, the knowledge gained from such a study will contribute to 
Z the confidence of those responsible for the operation of the works when the 
actual water supplies to be regulated will assuredly be different from those a 
= which have occurred in the past which have been the basis for the develop- as | 
_ ment of the regulation plan. The hydrology of Lake Ontario is very complex mts 
a that a thorough understanding of it will be necessary for efficient, — } 
operation of the control works. Secondly, ‘such studies may enable the de- 
an of long-term and short-term forecasts of water supplies. This 
would contribute to the improved use of the available storage on Lake Ontario. ieee 
‘Thirdly, this knowledge can make possible an assessment of the possible ad- i 
_ vantages to Lake Ontario of the development of nme facilities on sonthe — | 
This study was prepared with these considerations in ha 
_ purpose was to investigate in a general way all the factors sathoencing the a 
— flows to Lake Ontario and to set them forth in a concise form, The effects ot 
Precipitation, temperature, lake areas and the physiography of the drainage 
_ basin on the character of the inflow are analysed. From the better under- .: 
standing of inflow characteristics gained from this analysis an assessment is 
er 
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warfare... The excess water from the le upper» Great Lakes, i.e. Superior, 
_ Michigan- Huron, St. Clair and Erie, p passes through Lake Ontario _ and thence a 
to the St. Lawrence River and the sea. This water which is transmitted to ey 


Lake Ontario from the ‘upper Great Lakes is of great importance to an under- 1: 


standing of the. hydrology of Lake Ontario since it 


basin of the 2 presents a summary of the the pertinent 
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recasts of water 
| =‘ The Great Lakes constitute one of the finest natural regulatory systems in e 
- a the world. At the outlet of Lake Ontario the total drainage area of the sys- 
a at ‘3 y 
& ee The characteristics of the water supply to the upper Great Lakes are sub- — a 
to the divergencies in climate, topography and geology that can occur in 
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glacial periods had a great ¢ effect on the entire region, ‘grinding and 


For comparative purposes the | outflows from the various lakes which are 


an area of 260,000 square miles extending 750 miles from east to west, and 
aa The region contributing to the water supply of Lake Superior and the “4 : 
‘northern shore of Lake Huron is generally rough, wooded country typical of — 
the Precambrian Shield. _To the south and east of Lake Erie lies the or 


9 _ mations lie the Great Lakes - St. Lawrence Lowlands which comprise a 


= greater part of the area draining into t the - upper Gr Great Lakes system. ‘The - 
_ denuding the Precambrian Shield, flooding the low areas and depositing silts, m 
and sands, while leaving morainic hills, ‘till plains, drumlins and eskers on 
the m more elevated parts and generally disarranging old drainage patterns. af ' 
a The upper Great Lakes region lies in the path of the prevailing westerly er 
7 a winds and is particularly subject to winter storms. In the lee of the lakes = ‘a 
areas of more moderate temperatures. and extremely high snowfalls. Gen- 4 
- erally, the temperature increases from north to south and the precipitation a: 
increases from northwest to southeast. The northern part of the d drainage © ce e 
- basin has a mean annual temperature of 31°F. with an average length of the — 
frost-free period of only 75 days. Near the western end of Lake Erie a 
mean annual temperature exceeds 48°F. and the average | length of frost-free 
eriod exceeds 175 days. average annual rainfall varies from 25 inches 
= _ west of Lake Superior, to cover 40 inches to the south and east of Lake Erie. — ; 
ei A Although: the inflows to the ‘upper | Great Lakes and the evaporation losses 
vary: widely from year to } year, from season to season and from lake to lake, 
q _ the diversified runoff characteristics and the great expanses of lake surface — 
contribute to a relatively uniform inflow to Lake Ontario through the he Niagara | 


x River. The > monthly mean ot outflows from Lake Erie » have v varied from a high | 


-e 254,000 cubic feet per second to a low of 121,000 cubic feet per second. ieee | 
shown: on Figure 2 and the outflows from Lake Erie just mentioned are re those 7 ae | 
. Certain adjustments would be 


and levels will include these adjustments. A net diversion into the — ~ 
ad Lakes system of 1900 cubic 1 feet per second i is presently in effect. T This in- is 


to east the greatest length of the drainage basin is 280 miles ‘and the greatest ie 


width from north to south is 240 miles. _ The lake has an average depth of fap- 


a average water level has been about elevation 2: 246 and the ‘maximum range of a oe 
onthly stages has been about 6.6 feet. : 


= 


&g 
thecessaly lO Make Wiese LiOWS CONSONaANL Wilh Present CONQILIONS OL 
UIC LUV CUNIL ct per Seu C al at — 
—— diverted into Lake Superior from wee 
- ne local drainage basin of Lake Ontario has an area of 34,630 square a 
4 ; including 7,540 square miles of lake surface or over 20 per cent of the el i 
_ drainage basin. Of the drainage area, 18,710 square miles are inthe 
— 


‘LAKE. ONTARIO 


The geological features which dominate the local drainage. area of 
= are the Great Lakes - St. Lawrence Lowlands, the Precambrian — te | F 
‘Shield, the Adirondack Mountains and the e Allegheny plateau. _ These features 
are highlighted on the contour map of the local drainage basin presented on 
ip <— Great Lakes - St. Lawrence Lowlands follow the shoreline of Lake ie 


- Ontario, striking inland for various distances. | The lower lands close to be 
s became covered with silts and — = the edge of this lower area — 


the Niagara escarpment forms another sennieny to this lower - lying. area. . In 

the higher areas the relief was formed by the shoving and depositing action =— ; 

of the glaciers which left behind the morainic hills, drumlins, eskers and aa 
“ plains. The glaciers also blocked off the established southward drainage of 7 


_ the Finger Lakes in New York State causing them to flow northward into | 
Ontario. 


ee The Allegheny plateau forms the northern edge of the Appalachian forma- a 


ir tion. - It extends across the southern end of the drainage basin in an east- west © 
- direction from Syracuse to Buffalo and is deeply indented by the Finger Lakes. 
_ The plateau slopes upward to the south and ranges in elevation from 1,000 ae 
2,300 feet above sea level. The terrain is very irregular and hilly and has 
‘been largely cleared of its original forest growth. _ The highest point that is ar 
achieved in this formation in the drainage basin is 2,550 feet above sea level. Fa a 
_ The headwaters of the Genessee River and the southern and western tribu- — 
taries of the Oswego River lie withinthis areaa 
a northern part of the drainage basin is located within the Precambrian — 
region. . Rock knobs « or hills with shallow soil depths predominate. ‘Little of 
4 the area is cultivated and it is still largely covered with forest. In the ex- i 
a treme north of the drainage basin elevations of over 1,500 feet above sea .” iv 
level are reached but the general elevation decreases 'to the south and east to 
‘the Thousand Islands Section of the St. Lawrence River. _ This region 
‘the headwaters of the Trent, Moira and Napanee Rivers. 
ang —_— Adirondack plateau is an elevated extension of the Precambrian Shield 


wooded with very rough topography. The highest hill in the drainage basin — 
_ reaches an elevation of approximately 3,700 feet although elevations of over _ 
aa 000 feet are reached farther east. The area includes the headwaters of the br 
_ Oswegatchie and Black Rivers and the northeastern tributaries of the Oswego 
River. valley of the Black indents the plateau ina striking fashion. 
ie Lake Ontario drainage basin lies along the track of many of the low ‘ 
; s areas which sweep across the northern part of North America from 


west to east. This results in stormy changeable weather with wide vertations i 
in n temperature | but remarkably uniform annual penne For the drain- 


a with tl the ne topography play a major part in the areal distribution of rainfall. The 
average ge annual precipitation varies from 27 inches in the area below the 
ee Niagara escarpment to over 50 inches in the Adirondacks to the east of ek 
lake. snowfall is even more with the winds picking up 
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quantity of ‘moisture as as th pass | over the open lake and depositing it 
e to the lee of the lake. The averse snowfall varies from 40 inches at the west i 


_ the lowest temperatures. The average annual temperature is about ASF wi ah th 
0 se 

extremes of temperatures of 105°F and - -40°F having been been recorded, 


analysis of the water supply conditions Ontario will be de- 


shorter term records with more comprehensive coverage. — tl 

a their importance to navigation, continuous water level records for the Great = fh 

S. Lakes and the St. Lawrence River have been kept since 1860. It has been 48. * a b 
| therefore, to construct records of the water ‘supplies and outflows © 

of the lakes for a period of almost one hundred years. Records of tributary 

and detailed climatological records, however, have been available for 

al 


t Total Wat t t 
ive ‘The net total water supply to Lake Ontario is defined as the inflow to the 


tary area plus the precipitation on ‘the lake surface less the evaporation re 

F the lake surface. The records of monthly net total water supplies have been > 
i based on the monthly mean discharge of the St. Lawrence River at troquois, B: 
sa" Ontario, with the water levels at the beginning and end of the month measure 
— selected gauges around the lake shore. A change of one foot in the eleva- 7 
tion of Lake Ontario in one month indicates that the difference between the < 


net total water Supply and the outflow has been 80, 000 cubic per | second, 


feet per — 
‘The net total water supply Lake ‘Ontario. be down into 
principal components, the outflows from Lake Erie and the net local water 
7 supply. The next total water supply has averaged 246, 000 cubic feet per : 
‘second with 210,000 cubic feet per second contributed from Lake Erie and Aa ‘al 


the balance from the drainage basin, The average monthly istration 


general the ‘net total water supply reaches a maximum in the ‘month of April 


supply. ——— The highest and the lowest lines are the maximum and minimum for ; 
3 each month regardless of the ‘year in which | they occurred, Generally, the 
‘ extremes: vary in the same manner as the averages. The e monthly mean net bh 
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is due tO EMect OF the Met Water supply. Figure o-b, 


total water supplies have varied from 416,000 cubic feet per second to 


145,000 cubic feet per second, ‘The crosshatched area in Figure 3 
a range of one standard deviation on either side of the mean. © The ‘Standard = -, 
deviation varies little throughout the year. 
Month- to- ‘month variations in the net total water er supplies are 
oj quite small in the late summer and early fall but can be very large during the ‘a 
ae and spring. An extreme value occurred in 1936 when the variation i in yh : 
the monthly means from February to March was 178,000 cubic feet per 
second, _ An example of extreme variation for three months occurred in 1913 
when the monthly mean net total water supplies for January, ee and — 
March were 324,000, 220,000 and 326,000 cubic feet per second. 
Ms _ The variation in the mean annual net total water supplies is shown ona 
‘ frequency distribution on Figure 4. ‘The for the yearly 
= values has been 93,000 cubic feet per secon¢ Pp 


oy The greater part of the excess water from the upper | Great Lakes passes 


a through the Niagara River into Lake Ontario. To obtain the total outflows eS > * 
ie: from Lake Erie the quantities of water diverted for power and navigation have - a 
been added to the Niagara River flows. _ These total outflows have also been -. 

adjusted for a continuous net into the the upper Great 1900 
. The flow in the Niagara River is essentially dependent « on the water levels 
of Lake Erie. Sr There are three categories of variation in the outflow from 
Lake Erie, i.e. transient, seasonal and long term. 
‘The transient variations are caused by winds and seiches (oscillation of 
_ lake waters due to changes in barometric pressure) which may affect the | 
ater level of Lake Erie at the head of Niagara River by more thanten feet, _ 
and by ice jams forming at the outlet of the lake at the head of the Niagara © vs 


ue River. A large variation in flow resulting from | a recent storm is — on 
ae variations in the outflows from Lake Erie are small due to to the a 
moderating effect of the storage on Lake Erie and the other upper Great ees : 
Lakes. _ Generally the outflow reaches its peak during the month of June and 
its lowest value during the months of January or February. The monthly | 
pattern of the average outflows is shown in Figure 3-A. 
‘The diversity of hydrologic factors affecting the water manatiii to the up- 7 
per Great Lakes and the retarding effect of the great volume of storage avail- — 
able, result in long-term variations in the outflows from Lake Erie. Short- 3 
term 1 variations in water supply a are absorbed by the storage with little effect 
F on the outflows of Lake Erie but long-term variations gradually change the . 
outflows so that the seasonal and transient variations are ain —: 


“Lake Erie from 1926 to 1935 is shown in Figure 6. Also shown 

a hydrograph of average outflows and a plot of the accumulated deviations 
from average of the precipitation on the drainage basin of the upper Great 

F Lakes. The flow gradually built up from much below average in 1926 to near =< 


b _ the n maximum nm recorded in 1929 and then n decreased to the minimum recorded My 


in 1934. These long-term variations were obviously in response to long term © 
variations in the on the drainage basin, ‘The figure way: clearly 


pail J 
iJ 
a 
ed 
= 
uc 
— 
x 


_NWET TOTAL 


IN MEAN WATER SuPPLES 


TO LAKE ONTARIO 1860-1964 


WATER SUPPLY 


COMPONENTS OF — NET TOTAL wer SUPPLY | ERIE me LOCAL WATER SUPPLY 


RANGE OF TWO STANDARD DEVIATIONS 


MONTHLY 


OUTFLOW FROM 


a LAKE ERIE 


4 


FREQUENCY DISTRIBUTION 
ANNUAL waTER 


1860 - 1954 


OUTFLOW FROM 
ERIE 


= 
= 


7 


LOW FROM LAKE ERIE-iOOO CFS 


ouTF 


N 


i 


| 
: 
—_ ~Bri 


FiGus 


HYDROGRAPH 

| FLOW OF NIAGARA RIVE 
DURING THE PERIOD OF 

NOV.7 TO10,1957 


EFFECT OF PRECIPITATION 
ON UPPER GREAT LAKES BASIN 
ACCUMULATED DEVIATION 


OF PRECIPITATION FROM \__ ON OUTFLOW F ROM LAKE ERIE 
LONG TERM MEAN 


"PRECIPITATION IN INCHES 


= 


_-WYDROGRAPH OF OUTFLOW 
FROM LAKE ERIE 


OUTFLOW FROM LAKE ERIE-I000 CFS. 


— 


: - much above or below average. _ Although the precipitation for 1932 and the 
Sy early part - of 1933 was average, the deficiency of water supply to the upper — i 
Great Lakes which occurred in the extremely dry period in 1931 had so de-- 
pleted the storage on that the outflow from Lake to 


had 


The trend. of these. extremes follows the trend of average ‘seasonal ell 
quite closely. The monthly mean outflows have varied from a maximum of 

y 258, 000 cubic feet per second to a minimum of 131,000 cubic feet per second, 
or a range of about 60% relative to the mean. The minimum occurred in bic z 


* 


variations in mean outflows from Lake Erie are 
- ‘Figure 4 in the form of a frequency distribution. The maximum variation of a Se 

the annual been 68, feet per second 32% the 


The net water supply to Lake is ‘the difference between 


ie aa of the land runoff from the local drainage basin of Lake Ontario — 
7 “3 precipitation on the lake surface less the evaporation from the lake sur-_ 
face. The seasonal distribution of the average net local water supply is 
id shown on Figure 3-A. The maximum monthly average was 90,000 cubic feet 
Ar per second in April and the minimum was 6,000 cubic feet per second in vt ; 


a seasonal distribution « of ‘the maximum and minimum net local water Bas 
supply for each month is illustrated on Figure 3- D. It 

sii ing the late summer, fall, and early winter, the minimum n net local water 
_ supplies are negative quantities. These negative quantities result when —= 

evaporation from the lake surface exceeds the sum of the precipitation on | a 

the lake surface and the land runoff. The net local water supply varies widely — 

with a maximum monthly mean of 188, 000 cubic feet per second and a mini 2g 

_mum monthly mean of | -32, 000 cubic feet per second or a range of about 610% 4 

“of the mean, The seasonal distribution of the extremes bears a close rela-— 
to the seasonal distribution of the averages. A range of one standard 
_ deviation on either side of the mean for the net local water ‘Supply | is also ’ ; 
a a shown. The standard deviation varies from a maximum in April to a mini- is bs 
| _ ‘mum in August. _ The standard deviations for the net local water supply aver 


ee The net local water "supply ¥ varies 8 widely from month to month, notably a 


he variation from February to iar was 138, 000 cubic | feet per second, 
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During the year I9IS the mean monthly net local water supplies varied rom 


114, 000 f. in n January, to 5, 000 « fis Ss. in February, 110, 000 cubic fe toa 


a Figure 4 shows the variation in the annual mean | net local water ote ‘as 
as a frequency distribution. The maximum variation of the annual means has 
been 39,000 c.f.s. or 109% of the average net local water supply. varia~ 
-¢ tion is 62% of the maximum variation in the annual outflows from Lake Erie. : 
;j The slope of the frequency distribution for the net total water supply is greater 


_ than the slopes for either the outflows from Lake Erie or the net local water 
supply which indicates that there is a tendency for both components to be high i 
or low together. However, the extreme values of both have fe 

Amore . thorough understanding of the \ variations in the net local water 
_ supply requires a more detailed study of the land runoff in the local drainage eo 


and the precipitation- balance on the lake 

_ Lands runoff is the inflow to Lake | Outario from the local drainage area ex- 

2 -_clusive of the lake surface. This area is 27,090 square miles. 


The major 


tributaries are the Oswegatchie, Black, Oswego, Genessee, Trent and as. 
_ Rivers which contribute more than half of the net local water supply to | 
: Ontario and form the basis of the analysis of land runoff in this study. mt. a 


_ The Oswegatchie and Black Rivers drain a large part of the Adirondack 
plateau within the Lake Ontario drainage basin. The fl flows of | both tributaries ir 


q waters. The Oswego I River is the largest tributary of Lake Ontario saiteaere ‘ 

. the ‘Niagara River. ¥ It drains the southeastern part c of the drainage basin oe 

with headwaters in both the Adirondack and Allegheny plateaus an and serves as_ 

- an outlet for the Finger Lakes. Its flow is regulated appreciably in conjunc- i 
tion n with th the operation of the New York State Barge Canal. The Genessee > ole 
River rises in the . Allegheny plateau of Pennsylvania and the south-_ 

central part of the local drainage basin. Its flows are regulated slightly by - 

i upstream storage reservoirs. ‘The Trent River rises in the extreme northern — 

a part of the ¢ drainage basin in the Precambrian Shield and drains the north- ae : 

‘ central portion of the local drainage basin. It has numerous lakes along its 

course an and the flows are regulated fi for power and for operation of the Trent 
Canal system. m. The Moira River has its headwaters in the Precambrian 
s Shield to the east of the Trent River and flows through a region with very _ ee 
little soil cover. The flows are not regulated, The period of record | for the i 
downstream gauging 1g stations on these rivers varies from 26 to 45 years. PAsSee 
_ None of the major tributaries arise in the Great Lakes - St. Lawrence 2 
wae to the west of the lake or in the northeastern part of the drainage ve 
basin, Two smaller rivers, the Credit and the Napanee, are considered to me a 
be representative of these areas. The Credit River rises above the Niagara — a 
_ Escarpment and flows into Lake Ontario west of Toronto. The Napanee drain ed 
a. age basin is similar to that of the Moira. The flows of the Credit and Napanee 
‘Rivers are substantially unregulated. The periods of record at the down-— 
_ stream gauging stations are 10 years for the Credit River and 21 years -* > 
‘The density of stream measurement stations on the local drainage oe is = 
- _ adequate f for the determination of the areal distribution of ‘runoff. Figure 7-D | 


pice A 2 a map of the drainage basin showing the distribution of mean annual sll 


in inches, The area of maximum is of Ontario in the 
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of the the Black and © Oneida Rivers. Here the. » runoff ii is as 
2 high as 35 inches per year. The area of lowest runoff follows the shoreline - =e 
of of the western half of the lake where the average runoff is as low as 10 a i a 
_ Generally the runoff increases with distance from the lake. The reason for _ 
_ this distribution is apparent from inspection of Figure 7-B showing isohyetals — 
of mean annual precipitation. From a comparison of Figures 7-B and 7-D, = 
- may be deduced that the annual average water loss is generally about 20 inches E 
yr, although it is somewhat less to the east of the lake. It may be seen that the 
; ‘= of greatest precipitation and runoff are coincident with the highest land. Re 
In addition, the prevailing westerly winds, after crossing Lake Ontario, are 
4 affected orographically by the Adirondacks. This is especially true in the : 
- winter time where the winds sweeping over | the open water pick up large 9 i 
: quantities of moisture and drop it in the form of snow in the lee of the lake, 
‘The distribution of the mean annual snowfall is shown in Figure 7- a 
‘The seasonal distribution of average precipitation on the drainage areais 
- shown on Figures 8A, 8B, 8C and 8D. The areal distribution of precipitation ra 
during any one season generally follows the same pattern as that for the year, 


- @ remarkable feature illustrated by these figures is that the precipitation en 


does not vary widely from season to season. _ Figures 9A, 9B, 9C and 9D show 

areal variation of average temperatures for the months, January, 
and October respectively. The decrease in temperature with both lati 

tude and altitude is apparent. The range of areal variation of temperature is. 


g oa greatest in the winter. The seasonal and areal variations in runoff are shown . 


on Figures 10A, 10B, 10C and 10D. During the winter season the effectpi 
- areal variations of temperature on the runoff is very marked. In the northern 4 
a sections of the basin and in the Adirondacks, where the temperatures are del ae 
lowest, the precipitation is stored in the form of snow so that the ‘runoff 
a as “ much less than the precipitation. In the lower and more southern part of the 
a: _ drainage basin a combination of higher temperatures, small evaporation a 
transpiration losses contribute to a runoff that is almost as great as the pre- ® 
; ( _ Cigihation, During the spring months the situation is reversed. In the areas 
f _ where the precipitation was stored in the form of snow the runoff is either = 
_ greater or slightly less than, the precipitation, whereas in the other areas it 4 te 
is considerably less. During the summers the large evaporation and trans- A, 
- plration losses result in a generally low runoff. In the fall months the water 7 


losses decrease so that the runoff becomes g generally higher . The areal vari- 


cantly. Thus, , although the precipitation does not vary much from season to 
season the effects of temperature and topography ‘cause large seasonal fluc- . 
. Extreme high and low monthly mean precipitation and temperatures at five 
stations are shown on Figure 11. range of one standard deviation on either 
4 side of the average precipitation is also shown. The stations at McKeever, hill 
_ Alfred, N.Y., and St. Catharines, Apsley, and Kingston, Ontario, are con- one 
sidered to be representative of the Adirondack plateau, the Allegheny plateau, E 
the western end of Lake Ontario, the Precambrian Shield and the eastern end — a 
not vary much from month to month. The higher summer average pre- 
ee | Sa at Alfred indicates a somewhat different climatic condition atthe 
Lt _ extreme southern end of the local drainage basin. ‘The greatest variations in 
late summer or autumn, po possibly que to 
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ee only the maximum and minimum monthly mean temperatures 
are shown, the degree of variation of these quantities is an indication of the | ae 
variation that can take place within the month. The winter period exhibits the ‘. 
widest variations in mean temperatures, — - The minimum and maximum month - : 
- temperatures for all the stations follow the same trend. The similarity i in a 
variation between two widely separated stations i.e. McKeever 
= Apsley may be noted. Apparently the effect of increased elevation = 
McKeever balances the effect of the northerly location of Apsley. 

_ For a more detailed study of the seasonal distribution and variation al the 
land runoff it is necessary | to examine more closely the behaviour of the ae 
dividual tributaries. On Figure 12 is presented a hydrograph of the mean ~ 
_ monthly flow in inches for various tributaries and for the net local watery 
supply. Also shown is the percentage | of the mean annual flow for each month. 
The much higher unit runoff from the Adirondack plateau into the Oswegatchie 2 

7 and Black Rivers is immediately obvious. The lack of a deep soil cover and = 
+ storage reservoirs in the Moira and Napanee drainage basins has a marked ae 
. effect on th the distribution of of runoff. Much higher percentages of the annual ~ an 
pes _ flow occur in the spring and much lower percentages in the late summer than 5 
with the other tributaries. _ Moreover, on these tributaries there is a smaller — 
iz percentage of the annual runoff during the winter which would indicate - 7 
5, greater accumulation of snowcover, This would tend to increase the spring © 
- runoff. Throughout the greater part of the local drainage basin of Lake . 
4 Ontario the peak flow generally occu 
of Credit Rivers in the western part of the basin the spring season ‘seems 5 to - x 
- more advanced, The Trent and Oswego Rivers have the most uniform dis- | 4 
tribution of average monthly flows because of their large size and of the ef- hE = 
fects of regulation. With the exception of the Genessee and Credit Rivers the | 
4 average flow for January for all the tributaries is higher than for either : 
" - December or February. This would give some credence to the local tradi- 
tion of a “January thaw.” The runoff pattern for all the tributaries is, in : 
q spite of these differences, generally similar to that for the net local water _ 
supply. The difference can generally be attributed to the effect | of ait hid 
a -precipitation-evaporation balance on the lake. The maximum and minimum _ 
x for each month are compared to the mean flows for each month on 
‘Figure 13. The monthly distribution of the standard deviation is also shown, 
The figure illustrates that the range of variation is smallest in the summer 
_ months and greatest in the winter and spring. The reason for the small — 


range in variation of runoff is that the evaporation and transpiration abe. 


_ increase with the availability of water. The decrease in the effect of ol BS 
_ losses in the fall allows a much greater range of runoff. The effect of large a 
: (ustedes ae on the snow that has fallen and the small water los- 
ses tend to make the variability of runoff during the winter extremely high, ee 
_ The variability during the spring period results from the variation of a 


_ ‘Snow ¢ cover at the end of winter as well | as the | spring precipitation. Generally fer &§ 
te standard deviation follows the same pattern as average monthly runoff. eae 
_ It may be concluded from the foregoing that the variability of runoff is de- oO 
pendent not only on the variation of precipitation but also on the variation of ae 
temperature, A specific example. in which the precipitation and 
are related to runoff is shown on Figure 14. Deviations in precipitation and © - 
temperature at Belleville for the period 1941 to 1950 are compared graphical - 


ly with the deviations in flow of the Moira River. Belleville is at the me: of of 4 
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“the river ond it is believed that the deviations of the precipi-_ 
¥ tation from average reflect conditions throughout the watershed with cataciend 
Although there are many interesting occurrences shown on this plate ne 
¥ three periods will be examined closely, These periods have been chosen to | 
_ illustrate the factors influencing the variability of runoff which have been dis-_ 
: cussed results of the examination can be | summarized 
a 


‘Precip: Temp: | Runoff 


June transpiration losses reduce the 

1945 August Decreased evaporation 
September 01.61 -piration losses allow the increased 
October rainfall of September and for 


aa “0.01 4 Some snow melt, 
8 42.36 All snow melted 

March -0.% <-4 -79 Accumulation of snow, 
February and March. 


‘It may be analysis | of runoff 
- the local drainage basin of Lake Ontario lies in a transition zone with respect — 
_ to winter and spring runoff. To the north of the basin the snow cover is con- 
_ tinuous throughout the winter so that the winter runoff is usually low and 
_ form and the spring runoff high and irregular. South of the basin, continuous — 
Bs accumulation of snow is rare so that the runoff varies more directly with 
- variations in precipitation. In the drainage basin of Lake Ontario itself, ¥ 
* _ either of these conditions can occur depending on temperature variations 
_ during the winter months. It may be concluded that, because of the varia- — ee 
bility of temperature during the winter, the use of snow surveys for the ai we 
diction of runoff may have a limited use only, 
- ica Although instantaneous and daily peak flows are of little importance ina _ 


the drainage basin is of interest. Figure 16 illustrates a plot of the peak 

flows of the gauged tributaries against the ‘drainage area. The envelope of i 
_ peak flows and its equation are : also shown on the plate. —‘Itis: interesting to is 

note that the peak quarter-monthly net local water supply lies slightly above “f 


re the envelope of extreme flows for the tributaries. This may be explained by pe 
ie a the fact that the period of record represented by net local water supply is 95 


years which is much longer than the period for any of the tributaries of the — 5 
basin. The “maximum monthly 1 net local occurred in 


study of the hydrology of Lake Ontario, an indication of the flood potential oo 
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a a is derived from the 2 precipitation falling ont the surface as modified by ¥- 
4 evaporation from the surface. Since there is no direct method of measur- 
= either precipitation o1 or evaporation © on such large | lake e surfaces the clas 
_ precipitation- evaporation balance can only be found from a water budget study. 
Such a study has been made for the period from 1934 to 1952 inclusive. _. 
_ total land runoff for each month has been derived by adding together al all the 
; "measured inflows, with the exception 1 of the Niagara River, and estimating the 4 
of the land runoff on the basis of the known unit water yields, 
_ estimated land runoff was then subtracted from the net local water supply to 
“obtain the precipitation- -evaporation balance of the lake. Since it represents 
_ the difference between two large quantities, the precipitation- evaporation — i | 
eS for any one month is subject t to large percentage errors. . However, ie. 
“ee the average for each month over 19 years shows the average distribution wt * | 
i precipitation- evaporation balance with sufficient ee for the pur- . 
_ The average distribution of the precipitation-evaporation balance is shown. 
on Figure 15. This figure would indicate that the average contribution of the — 
_ surface area of Lake Ontario to the net local water supply is negligible and 
that the only effect of the lake area is to increase the net local water supply _ a 
a. during the seasons of high land runoff and decrease it during the seasons of 
low land ad runoff. average cumulative effects of the seasonal increases 
decreases are substantial, being of the order of 60,000 cubic feet per 
second for one month or three quarters of afootonthe lake. = © a 
There is no method of accurately analysing the two components of the | 
precipitation-evaporation balance. However, by using the precipitation 
Ps 5 i records of stations near the shore of the lake the apparent precipitation on ae 
the lake has been determined, By subtracting this from the precipitation- 
‘evaporation balance, the apparent evaporation was deduced. 
z _ Tt may be noted from Figure 15 that the average apparent precipitation — 4 
varies very little throughout the year. ‘The average apparent evaporation, 
however, varies considerably reaching its maximum value in October and its 
minimum value in June. is relatively uniform from Decem- 


because of the large of the lake and the constricted outlet 
conditions the changes have been modified and delayed. Thus, although a 
a range of one foot in the elevation of Lake Ontario can absorb 80,000 c.f.s. ce 
fe for one month a a one foot change in the water level only changes the outflow ue 


section nue resulted in overall changes i in the levels of Lake Ontario. “The | 
water which) will be us in discussion are the recorded valu 
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‘The hydrologic factors influencing the water supply to Lake Ontario have 
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— i i i 
stream from Prescott, Ontario. Previous to the construction of the | se: 
__Lawrence Power Project, this was made up of a natural rock ledge across" vay 


LAKE ONTARIO 


and for the diversions. During the winter, ice is. 
r, formed at and above the control section restricting the outflow and causing _ ne 
_ increases in the water levels of the lake. Temporal changes in the outflows - x 
: and water levels of Lake Ontario are caused | by winds and seiches. ‘ Because 
_ of its location and greater depth, the effects of winds and seiches on Lake _ ” = 
Ontario are not as greatason Lake | 
> To obtain a more detailed knowledge of the interaction of the outflows sien “i 
Lake Erie and the net local water supplies and their effects on the water ae 
h levels of Lake Ontario in the past, a further investigation was made. It was = 
_ first assumed that the outflows from Lake Erie remained constant at their 
- average value but that the net local water supplies varied as they did in the ; 
past. These hypothetical net total water supplies were then routed through q 
Lake Ontario using standard routing procedures and taking into account the 
~ actual effects of ice formation that had occurred. The same procedure was po 
used with another set of hypothetical supplies which were prepared by as- nite 


suming that the net local water supplies remained « constant at their average _ 
Tan but that the outflows from Lake Erie varied as they had in the past. The z 
obtained were thus directly comparable with the recorded adjusted 
= levels of the lake and are summarized in the aes og syd 
Kater 


from ‘Lake Local Water 
constant Net constant Outflows from 


_ ‘The » marked decrease in the . average : annual r range nge of stage, ‘compared iwith 
- recorded adjusted conditions which would have occurred had the net local 
5 water supplies been kept constant, indicates that the variations in the rota 
local water supply cause most of the annual variation in water levels. Con- a 
versely, the slight increase in the average annual range of stage that would — 
have occurred had the outflows from Lake Erie been kept constant shows that a 
- the seasonal variation of the outflows from Lake Erie is not in phase with on =i! 
Seasonal variations in the net local water ‘supply. effect of the long term 
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Maximum range of stage with respect to recorded adjusted conditions that 
‘he | would have occurred had these been kept constant, The relatively greater de a 
in maximum range with constant net local water supplies shows at 


"1959 


the 1s of net local water supplies. have had a greater on the 


variation of water levels than have the variations in the outflow from Lake “ay 

_ = Erie. Since the variation ¢ of annual average e water s supplies is smaller for the — 

net local water supplies than for the outflows from Lake Erie, it can be con- — 
cluded that the e seasonal \ variations of the former are of considerable impor- 


tance in 1 producing the maximum range of stage. 
a is difficult to evaluate the effect of the variations of the different sources 


ditions can vary in 1 length from three months to heen years. In the pm bi 
which have taken place it has been found that the critical recorded adjusted Ai 
levels can be used as an indication of critical conditions for regulation. With — 
a this criterion and from the results of the foregoing routing study, it may be — bey 
ae: concluded that the variations in the net local water supply will cause most of ; 
hig the difficulties in regulating the outflows and water levels of the lake. It can ° 
also be concluded that regulation of the upper Great Lakes which would mere- | 
oo ly smooth out the seasonal fluctuations in the outflows from Lake Erie would a 
a not assist in the regulation of Lake Ontario although if they could be changed © 

radically to offset the seasonal variations of the net local water supply some a 
improvement might result. A system of regulation which could smooth 
the long term variations in the outflows from Lake Erie would benefit con- = 
jae the interests concerned with the regulation of Lake Ontario. Con- > 4 

_ sideration of the vast volume of water involved in these long term variations — 


_ indicates that an appreciable improvement would be very difficult to a 


ae _* The regular seasonal acai in the outflows from Lake Erie and ae 


we 
- delayed response of the long term variations in outflow to the variations in 


i precipitation on the drainage basin of the upper Great Lakes indicate that it i 


is possible to forecast the outflows from Lake Erie for a long period in -— 
vance with a considerable degree of accuracy. However, with the extreme san 
F variability « of net local water supply and its dependence on variations of tem- - 


7 ‘perature as well as of precipitation, accurate long term forecasts of this ‘a ice, 
component could not be expected. With the net local water supply such an — a 
important factor in the regulation of Lake Ontario, it would seem — 
ble to devise a method of regulation that was based entirely on forecasts. — 
Forecasts may be of use, however, as a supplement to regulation. “Such 
forecasts may be of short or long term, The accuracy of short-term Such 
casts of water supply would depend upon the accuracy of temperature and _ _ 
precipitation forecasts. Long term indications of future supplies may be de 
termined on a statistical basis utilizing the observed, but not invariable, ten- 
_ dency of the annual mean outflows of Lake Erie and the annual mean net local | % 
water supplies to be high or low at the same time. Because forecasts of suf- — 
ficient accuracy would be of too short a period and the long term m indications nf 
would not be sufficiently accurate, the use of these techniques in the regula- _ 
at tion of Lake Ontario would require limitations on the use of storage estab- 
lished by testing over the period of record. odt 
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corm 
components. The outflows s from Lake Erie ‘provide | a a very large and ‘stable -. 
base flow to Lake Ontario which increases or decreases predictably inre- | 
sponse to long term changes in hydrologic conditions on the upper Great ‘a — 
Lakes. ‘ Although the ratio of the - variations is small in terms of mean out- = 
flow, the long term variations do represent an extremely large volume of - 


_ water. , The net local water supply, although § small in proportion to the out- 


supply and water levels of Lake Ontario. Although the average seasonal 
tribution of precipitation for the drainage basin is very uniform other factors. 
cause the net local water ‘supplies to vary widely. - These factors are the vai vari- 
ation in the areal distribution of precipitation on the local drainage basin, the | eu 
seasonal variation of water losses, the effect of the winter temperature vari- 
ations on the accumulation of snow and the effect of the seasonal variations of 
_ the precipitation-evaporation balance on the lake which intensify the effect of a 
7 seasonal variations of the land runoff. In addition to yielding an understanding ~ 
: of the water supplies to Lake Ontario the results of the study have indicated 7 an 
that the possibility of making accurate long range forecasts of water supply is 7 
rather remote. It has also been n indicated in the study that the regulation of | oy: ‘ 
the upper Great Lakes in . such a manner as to obtain an . appreciable benefit to — 
‘interests on Lake Ontario would be very difficult, 
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of ‘the American’ Society of of Civil Engineers 


EXPERIMENTS IN A TRIANGULAR CHANNEL* 


Ralph V Ww. W. Powell,” M. ASC ASCE and ‘E and Chesley J. Posey, M. ASCE = 


A 400- ft adjustable - triangular flume was tested smooth and also 
_ roughened with small rectangular battens at various spacings. wee rand the 


is made of Manning’s formula with some more recent suggestions, and the — 


velocity distribution, secondary currents, and other phenomena a are reported 


This paper reports the first phases of a of 
to flow in open channels being conducted at the Rocky Mountain Hydraulic ¢ : 
Laboratory, Colorado since 1954, major part of the 


which we were received through the Colorado State with Maurice L. 
_ Albertson, 3 M. ASCE as co- director. Support has also come from the Water 
_ Resources Division of the United States Geological Survey, the State Univer- 
‘sity of Iowa, the Ohio State University, the of Kansas, and 
7 a The triangular shape of flume was adopted to maintain geometrical simi- "4 
larity at all depths. Several others, such as Rogers, (5) varwich, (9) Owen (3) a 
and Moeller-Hartmann(15) have tested the same shape, | but their work was 
_ done indoors w: where the length was necessarily limited, The authors feel “= 
with a short flume it is very difficult to tell when the flow is at normal depth, _ 
na since the upstream and downstream depths may be the same as nearly as can 
be measured, for quite a te a range of de depths. (This point is discussed more fully 
| Role Discussion open until October 1, 1959. To extend the closing date one month, a 3 e ; 
= written request must be filed with the Executive Secretary, ASCE. Paper 2018 is aa 
: _ part of the copyrighted Journal of the Hydraulics Division, Proceedings of the ood 
American Society of Civil Engineers, Vol. 85, No. HY 5, May, ~~) One a &§ 
a, Presented at the October 1958 ASCE Convention in New York, N. Y. Mies “a t 
Emeritus Prof. of Eng. Mechanics, Ohio State Univ., Columbus, Ohio, 


= of Iowa, Iowa City, lowa. ety caer 


> Research Foundation, Colorado | State Univ., Fort Collins, 


— 
4 
} I 
Wy 


tater.) 3 Therefore this flume was | built o outdoors and made 400 feet long. sf It is 
a believed to be the longest channel systematically tested since Darcy and ne, 
¥ Bazin(1) and quite certainly the longest one of variable slope. As described - 


‘The purpose of. an 1 open- -channel friction formula is is to express the sai = 
between the channel shape and roughness, the discharge or velocity, 


’ ~~ of flow, and the bottom slope for the condition known as steady rsaier 


~ flow. Since the yearly cost of open channel construction in the United States — 
undoubtedly totals many millions of dollars it behooves us to make sure that 
we evaluate this relationship as accurately ; as possible, introducing minor 
variables such as the temperature, if necessary, and even increasing the 
complexity of the computation if that would yield a better result. a 
i: ‘Despite the wide use of open-channel friction formulas, it should be noted 
: _ that the condition of steady uniform flow is seldom encountered, either in the 
or in the laboratory. The depth of flow for Steady ' uniform flow, 
as the normal depth, is an abstraction, an extremely useful tool without which — 
_ flow profiles in open channels could not be predicted with any y satisfactory 
— determining normal depths in the laboratory, however, ‘itis essential 
that steady uniform flow be attained just as closely as possible. Unsteady 
flow may differ imperceptibly from steady flow in all other respects but the 
_ depth of flow. In the present experiments, not only was the discharge intro-_ ay 
duced into the flume held as constant as possible, but the water-surface pro- 
files were repeated as many times as necessary until they reached equilibri 7 
um, as far as that could be determined by the precision of setting and reading 


the point gages. For conditions of flow, this would take ‘minutes 


After the flow is is even more difficult to be sure that it is uni- 


- form” Backwater curves of the M1 and M2 type, corresponding to slightly e- 7 
ferent tailgate settings, may be very nearly parallel to the bottom, and yet ail 
_ have an appreciably different depth of flow. That the flow is uniform cannot 
ol safely be determined by reading gages | at each end of the flume, for the on 
in setting and reading only two gages, especially if the flume is short, may ee 
seem to indicate uniform flow whereas actually the flow may be well along an 
Ml or M2 curve, and the apparently uniform depth far from t the desired nor v4 
ssn lesser degree, the same thing is true for flow on steep slopes. Although 
» 4 the S2 and S3 backwater curves are not nearly as long as the I M1 and M2 — oe 
: mi curves, there are. apt to be standing w: waves on the surface that make ‘accurate 
getting of point gages more difficult. 


a deter mining the normal depths, therefore, profiles were obtained both 


A. &e 


at 


or : 


+4 


a 


= and above the normal depth, each profile for a steady flow 

are > obtained from of the tail at the lower end of the 
ia , and those of Fig. 2 from different settings at the entrance. Notice th a 


‘ 
i 
» 
fi q 
> 
— oF 
we 
— 
— 


a 

= 


5 


= 


discharge | too great, 


> 


- 


ate 


| 
1 
y | | Water -Surraces 
ougn TIME | peerr ) VERTICALLY As 


ey 
mM | 
~ 


tak 


pos 
smesth | (ne bettems) 


June26, Juv6,57BATTENS cTRs, 
PROFILES ARE DISPLACED VERTICALLY AS SHOWN BELOW ae 
£06 -/: 308 3.38-358 
10:56-125 | -1.300 
-0.500 0.00797 
10:37--AS 
1:25-1:44 
22 


> 
@ 


no 


co 


= 


tae 


| 


CHANNELS 


: fairly well aeetiat by a an M1 run, an M2 run, and two runs “ens near the 
- normal depth. This is about the maximum depth that could be obtained in the 
_ flume at the setting of that date, for no normal depth could be determined for 
slightly greater discharge of 9:26 a.m., July 15,1957, 
‘The 1:54 - 3:06 run of June 26, 1957, shown in Fig. 2, gives a normal depth | 
; of 1.45 feet which is about as great as could be determined at this roughness _ _ 
~ and slope. Of course, for the smaller depths ona steep slope, the S2 and S3 ; - 
curves are comparatively short, the normal depth is soon reached, and it was 
= sufficient to to measure only the S2 curves starting f from the critical 


 Whena test is being started in the field, with the depth for uniform flow 


” shooting. “To be sure, if it is well within e either ri range the characteristics ode | 

of the flow will immediately show which one it is. Flow on the critical boun- 

_ dary range, however, does not have any well- -defined characteristic except ia 

surface instability. Hence, whenever there was any question, profiles were — ig 

_ taken with different headgate settings and also with different tailgate settings. a 
7 This was not done with intent to explore the interesting backwater curves <7 
this vicinity, but rather to establish t the depth of uniform flow as reliably a 


' Fix Water is diverted from the North St. Vrain Creek into a forebay, through 

a metal flume, and intc a weir box 42 feet long, 8 feet wide and 7.5 feet deep. ‘nil 
: "See Fig. 3. The flow can be controlled at the point of diversion and there is z 
ss emergency spillway for protection against high water in the forebay at the 

_ side of the metal flume, which latter can be barricaded against extreme high 
water. Flow is ' stopped. by inserting the needles at the point of of diversion and 


_ wasting leakage over the emergency spillway. It is started again by closing eas 
the emergency and removing enough needles to obtain approximately 


4 


___ Discharge from the weir spilled into a deep tumble bay f from which i 


tumble was 6 feet the bottom of the flume. 
‘The variable - -slope flume had a triangular cross- -section n with central angle 
G ‘that did not vary more than 0.0025 radian from ninety degrees. The bisector 
_ of this angle was set vertical; in other winds, each side was set at 45° a 
i the horizontal. _ The walls of the flume were made of four-foot by eight - -foot 
3 panels of one- -half inch fir plywood, exterior grade, stiffened at one-foot cen- bs 
ters by fir two-by-fours. _ The two- -by-fours along the | bottom edge were i "4 
7 set slightly so that the reinforced panels nested and formed a a tight joint, — ; 
‘ which was held together by 1/4 inch bolts at one-foot centers. When the oe % 
3 were fastened the correct distance apart by struts and heavy wire ties spaced — . 


two foot eight inch centers -foot section o of flume was was formed 
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needles in a wasteway in the side of the weir box. Only a small amount 
lischarge was permitted at this station, the operator of which couldseea 
| of temporary hook gages fastened inside the weir box, and adjusted to 
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Plan view of the variable-slope flume 
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TRIANGULAR CHANNELS 


which could be joined its neighbors by splice plates. Fabrication was with 


The flume was supported on heavy timber bents at twenty- -four foot audio. a 

“7 The upper surfaces of the supporting timbers were set to a common plane at e 
an inclination of 45° with the vertical. Bearing blocks under the reaction stif-_ 
 feners of the “bent -plate ” girder formed by the flume sections rode up 

_ down the 45° inclines, as controlled by jacks resting on the inclines. _ The 

flume was, therefore, approximately to a true grade whenever it was in align- 

ment, Final adjustment perpendicular to the incline was made by driving 

ae shingles under the bearing blocks. Before water was introduced additional — 

bearing blocks were inserted under the flume at the third- -points between he 

~ inclined supports, and wedged tightly above vertical logs carrying the load . 7 
to wooden plates set directly on the ground. To prevent tipping, the ends 


upper flanges of the bent-plate girder did not have sufficient stiffness perpen-__ — 


| dicular to the plane of the side, they had to be braced at the third points, This , 
also prevented the flume from tipping. 

eg Setting the flume accurately to a new slope required three days at first, but | 
= improved methods the time was cut down to little more than two days. ‘ aa 
‘The flume was quite stiff, and raising one jack only half an inch would lift it~ . 

* from supports 24 to 48 or even 72 feet on each side without causing distress — 

_ or noticeable distortions at the joints lifted. It is hoped that when the hydrau-— 

. lic tests have been | completed, the flume will still be in good enough — : 
to permit structural tests of the bent-plate girder tobe made. 

The plywood panels were given a coat of penetrating sealer and one coatof 
forest green en paint when fabricated at the hydraulic laboratory of Colorado &§ 
. State University in Fort Collins early in the summer of 1955. When it was - 2 

- desired to start testing in 1956, it was found d that t the inner surface of the | 7% 

a flume, which had been exposed to the weather, was beginning to check, lose 
paint, and even suffer some ply separation. Since the bearing blocks unier 

; = intermediate stiffeners were not of oak, as had been specified but found 2 _— 
to be unavailable, and the intermediate splice | plates not completely fastened 
due to a shortage of screws, the flume had developed a “sway-backed” con- | ; 

3 dition with sags at the intermediate points ranging from 0.01 to 0.03 foot. a 
Dare a were also slight local bulges (less than 0.02 foot) over over the main nop 
ports, caused by placing wedges under the stiffeners instead of under the 

While the flume was being realigned and the missing screws inserted, it 
_became apparent that the painted surface was suffering rapid deterioration. 

the altitude is 8000 feet above sea level, day- -time temperatures in 
the intense sunlight in the flume went as high as 130°F, while frost formed 


S on the surface at night. Several solutions to the problem were tested. = 


sealed with on wes Cr caulking compound, and, after | painting, 
covered with black vinyl chloride adhesive tape. The tape alone would have 
dd sufficient to seal the cracks, and perhaps better by itself as the asphalt — 
did not have as much flexibility, ; and weakened the adhesion of the tape. ens). 
= After a preliminary test with a maximum load of water, 
the bottom ee and the wires holding the top struts uts tight t were tightened. 
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CHANNELS 


4 During subsequent tests ne no leaks greater than a few slow drips were re observed. § 


: ae In July of 1957 a new method of bracing the flume that was being tried es 

an effort to cut down the time required for changing slope » overstressed the a 
- ties, and during a run with the flume nearly full, 160 feet of it collapsed. The 
merits of the design became 2 apparent | when it took k only a little over two weeks — 
tor rebuild, using mostly salvaged material. During the rebuilding greatly im a : 
‘proved methods of fabrication were developed, and a simple and fast 


_ of bracing and aligning the flume was discovered which solved that hitherto — * 


Specially -built point gages were installed at twenty-four foot spacing along _ 


| tte centerline of the flume. These gages were held vertical by adjustable 
| supports fitted with grooved guides and a atte so that the gage could easily 
a In order to read the. point gages } which were set to indicate the water | sur- 
face and special level rods which were used to measure the bottom elevations, 
three levels were kept set up alongside the flume. Intermediate bench marks — 
were used to establish the instrument heights accurately, and in addition the x 
pak gages midway between two instruments were usually read from both. bes 
‘The “headgate” of the variable-slope flume, which had to be used when the 
-velocity was ‘supercritical, consisted of a vertical bulkhead with rounded bot- 
tom surface made from a log and sometimes from two or three additional logs 
placed so as to decrease further the entrance depth. It was fastened in at the © 
_ upper end of the flume, just at the end of the curved entrance section. Due to 
the dissymmetry of the tumble bay the surface waves below the headgate were 
_ often unsymmetrical. Usually they could be reduced by fastening a baffle =a 
_ board at a location below the headgate that had to be determined by trial. . The 


F. unsymmetrical surface waves seemed to disappear within a hundred feet of mad eh at 
‘The “tailgate” of the variable-slope flume was simply a set of cone. - 

could be fixed in place as desired. 


4 the. practice to wait ‘until the flow became steady, as indicated by. the point 

_ gages, before recording any readings, and to adjust both head and tail condi 
a tions if there was any doubt as to which should govern, _ Occasionally itwas 
- found, upon repeating a set of profile readings, that changes had taken place, i. 
¥ indicating that the flow had still not become quite steady. _ The readings were = 

Calculation of Flows 

In general flows over. 4 cfs were with a 2 full- width s 

“crested weir and those smaller with a 90° V-notch weir, 
‘The full-width weir was of stainless steel 24 in wide and 7.948 ft long. ¢ petite 
was about 0.25 in thick but the upper edge was beveled at 30° to the vertical — et 
om. for a flat top machined to 3/32 in width the same as Schoder and li <f 
- Turner’s(3) Series D to P. % At first every time the zero of the gage was . 
termined, levels were taken at each foot along the crest. It was found that — 

the north end was 0, 004 ft higher than the south end, and the average of all | 
the readings was used to compute the discharges. But it was found that this oe 
; average was apt to be the same as the height at the middle, so that sometimes eats 


| ¢ 

| 

4 

a 

"od 

hen it — — 

— 
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only this was taken. The upstream face was as nearly vertical as it could be 
set with a carpenter’s level, and the weir plate was held in place by counter- : 
sunk screws, so that the entire surface was smooth, © The end of the w weir box | 
a asl below the stainless steel weir was of concrete and quite smooth and vertical, 31 
in the plane of the downstream face of the weir 
--' The 90° y- notch was cut in a sheet of 1/4 inch | aluminum -painted tempered 
iy masonite the same length as the full-width weir and 4 ft wide, which was ae 
fo | _ placed with its lower part resting against the stainless steel weir plate and 
its upper part supported bya wooden frame. The bottom of the notch was 
i. a about 0.3 in above the crest of the steel weir. The edge ¢ of the V was formed — 
_ to the same shape as the full-width weir. It was found to keep its sharp edge ei 
a wang well, a The ¢ ends and bottom w were sealed with tape SO as to make — 
i - head on the weir was read on a hook gage attached to a plank which a | 


was spiked t to two trees beside the weir r box. The hook gage well was a lucite 
i = cylinder about 5 in. in diameter which was was connected by about 4 ft of 3/4- inch 


box. The hook gage well was drained and refilled from the weir box whenever 
_ it was found that the water in the well was as much as 1.5°C warmer than the - 
ny water flowing over the weir. The “zero” of the | gage was determined by s ‘setting 
ha the hook to the cross hair of an engineer’ ’s level, and then reading a rod held — 
on the weir crest as mentioned above. . For the V-notch a special target was i 


“4 constructed which could be placed in the notch with contact along both edges 
_ When a run was being made the hook gage was read every few minutes and | 
- recorded with the time. These readings were plotted against time, and from — 
t this graph the average head for each run was determined. The median varia- 


a 0. 001 ft. But there may have been another 0. 001 ft error in the “zero reading.” 
= The tlows from the V-notch were jane magn by the formula given by Lenz(6) 

not than 2.485. R 20-5, 1.5 /v and Wis be 2/¢ where gisthe 
- surface tension in lb/ft. “Since the laboratory is located at an elevation se ce 


8000 ft above sea level, the value of g is less than at most laboratories. 
_ Upon the advice of Professor Weikko Aleksantri Heisknen of the Ohio State 
iJ University, it was taken as 32.136 ft/sec?, ONE 
— was originally planned to compute the flows over the full-width weir by / 
_ Schoder and Turner’s formula (D). 3) But when a velocity traverse was made 
in the weir box opposite the the hook gage piezometer t to determine Va and Vb * 
was found that the Vb was quite a bit larger than Va that is, the highest ap- 
proach velocity was near the bottom of the weir box. Schoder and Turner’s: - 
_ data show that when Vp is larger than Va, the flow is about the same as given 
Francis’ ’ simple formula, or sometimes even less, while formula 
always makes it more. As described elsewhere, (17), a of Schoder cal 


Torner’ 's data was made, to formula 
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= (3.08 + 0.44h/P +0. 15 r)(b - 0. 003)(h + 0. 003)!-5 


For our weir, b was 7.948 ft and P was’ 5.5 ft at the weir in bane and es, | 
i ‘and this value was used for runs 1 to 108. _ However, the floor was not level _ ; 
and the average depth of the floor opposite the piezometer connection below 
_ the weir crest was only 5.24 ft. Perhaps this should have been used for P for “ { 
these runs but the change in Q would have been a very small fraction of on 
«per cent. Before work was started in 1958 a new floor was installed with _— a) 
t P = 5.2 ft at the piezometer connection and 5.0 ft ft at the weir face. » Pv was . 
a _ In the formula, r is the ratio of the velocity of approach above the level of 
crest to that below. course this was not measured for each run, ‘From 


_ traverses at five different heads and from Schoder and Turner’s results, it 
was decided that r could be taken as 0.435 h + 0.870. (This was for 1957 and 7 
et) after additional baffling had been installed at the upper end of the weir = 
- box. For the data of 1956 r was taken as 0.435 h + + 0.199. This is a ; 
than the formula published in (16), so that the full- width discharges given in 
that publication have been revised.) leads to the 
Q = 7.945 (3,110 + 0. 145h) (h +0. .003)1- 
1987 7.945 (3. 210 0. 145n) +0. 003)-8 


(3.210 + 0. .150h) (ns + 0. .003)1. 


R. E. Glover of the U.S.G.S. made some diffusion by 
- the salt velocity method during the summer of 1956 and kindly supplied us my 
_ With two va values of the discharge which had also been — by 7 full- ai 
width weir. The results are as follows: 


USGS. 


gave tl the following rough che checks of the weir ‘formula: 


Date of in cis 
y weir 


24, 1956 No 17, -0.( 02. cfs .1% 


19 17.64 O7 cfs 6.5% 
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by weir formula by traverse 


Aug. 15, 1957_ 14, cis 0.41% 
17.76 0.14 cfs 


Determination of of Bottom im Stopes 


For each cole and roughness a number of runs were made covering a aie 

_ ‘range of discharges and a similarly wide range of depths of water in | ME aie 

fo The weight of water in the flume varied from a few pounds per foot 7 
as much as five hundred. Since the supporting system was necessarily sub- 

—_ to deformations from the varying load, and the = come of the main a 


# of the test runs, with the water flowing as near the normal depth as possible. 
To permit this, special steel level rods were made which could be held in ei 
_ position even during the high velocity runs and read from the three levels set 
up along the length of the flume. Readings were taken at eight foot intervals 
along the bottom. ‘The r rod probably interfered slightly with the water- -surface 
Bs profiles, but water surface -profiles were not taken while the rod was in place, 
i es _ and the rod was in the flume only part of the time. _ For some positions of ‘ag 
2 flume and levels, , the | supporting bents or other obstacles interfered with the 
line of sight, in which case the readings were simply omitted unless it was 
possible to include t them by “making” an additional level setup without too much 
Because of the time required, bottom -profile data were not taken for every 
a water-surface profile run. Usually bottom-surface profile data were taken 
. corresponding to the highest and lowest runs of a given series. The differ- - 
hy ence between the mean bottom slopes and elevations determined at these ex- 
i tremes were usually negligible or very slight, so that the values for inter- _ 
er mediate runs could be assumed to be the same as that of the nearer run, or 7 
_ The profile data were plotted in pencil on matte - surfaced acetate sheets, _ 
using a scale of 20 feet to the inch horizontally and an adequate vertical scale | 
= depending upon the slope. After each of the observed elevations had been inked 
in, using small circles centered on the plotted points, the straight line of best | 
_ fit was determined by shifting a straight line drawn in ink on a separate sheet 
z around under the acetate until its fit was judged to be as good as possible, — = “a 
B _ Then a fine pencil line was drawn on the acetate above this inked line, using 
a steel straightedge, , and its elevation at each end of the flume seas 4 
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obtained | by 50 individuals from i identical prints of the graph of data ¢ of 3:00 4 Bag 
ted 
ea the numerical data. He found the mean of the graphically-determined i 
slopes to be 0.008000, with standard deviation of 0.000030. Least — 
F solution gave a slope of 0.007998, with standard deviation (corresponding to — 
_ the scatter of the data) of 0. 0000144. _ It seems reasonable to conclude that ia 
the slopes to within one-l -half of one per cent. 


pee at twenty-four foot intervals along the flume were plotted on the same 
acetate sheets as the bottom- -profiles, and after the plotted points were inked 7 
the line best representing the normal depth was drawn by a similar technique = 

‘ taking care to make it parallel to the appropriate bottom-profile line. For the 
_ 1956 tests, this was done in the. drawing room after the runs had been com- 

7 pleted, but for the 1957 and 1958 tests the pencilled plots were made atthe _ 

- game time the observations were taken, the slope and normal depth co a 

tions being made later. All of these _ drawings for the 1956 tests are cette 
duced at reduced size in the R.M.H. L.’ s Report No, 21, and it is expected that 2 


“report. ‘Ther results are summarized in Table I. 

Brief History of ‘Search for a | for Resistance 


> 


recognized is V = = CV RS , Proposed | by 'Antione Chezy in 1775. He - evidently 
_ supposed that C would be constant for any given channel, but this proved not 
to be the case. In 1869 Ganguillet and Kutter presented a resume of porecngl a 
_ ly all of the measurements of open channel resistance that had been made up 
to that time and from them derived the formula which bears their name. fact 
was very widely used for many years, but its cumbersomeness and the fact 
- that it is not dimensionally homogeneous has led to its replacement by other 


ve a b i in 1897 he Cc Neither 


“these formelas seem to have been used in this country. 1868 Philippe 


Gaspard Gauckler, a French engineer, proposed the formula 


‘The first formula for the resistance to flow in open channels that is si Es : 


“ 
<J 
Pic. 


| 


4 but for some reason he limited it to slopes less than 0. .0007. This formula 
_ seems to have attracted no attention at the time. Hagen proposed the same 
formula in 1881, without the limitation, and Manning 2) repeated itin 1889 
without acknowledgement. But he was not satisfied with the formula and. 
instead recommended a much more complicated one. In 1894 Crimp and 


Bruges suggested V = 124 R2/3 s1/2 for sewers and water mains, a formula 


which i is in In Cc. H. ‘Tatton | proposed to the Boston 
ith | thn n to be the same 
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we today call the ae formula, it seems to have never attracted attention. 
ity It was not until books published by R R. B B. . Buckley in 1911 and P. A, M. ‘Parker — 


all other formulas for estimating channel resistance. 
ie After the work of Prandtl, von Karman and Nikuradse had put a new light i. 
2 the resistance to flow in pipes, , and had shown the important difference be- 
‘~— smooth pipe and rough pipe flow, it was natural to ask if channels could 
- not be treated similarly. _ Probably the first to propose formulas on this basis 
was Keulegan.(4) He was soon followed by Johnson,(7) Powell, (8, 10, 
i Kirschmer, (9) Robertson and Albertson(12) and others. However, none of — 
these developed formulas that are satisfactory for universal use. 2 fact it a 
_ became | obvious that when a channel acts as smooth arse is any roughness 


‘ - draulic radius and then use the Weisbach formula and the Moody « diagram for 
both pipes and channels. It is, of course, allowable to use the formula, but it a 
‘must not be expected that f will be the same at the same Reynolds number, 

ibe a and the same wall roughness, The free surface and hemi currents would 

be.expectedto increase the resistance, 

ae noteworthy treatment of the resistance to flow in rough channels is that “ 

¥ of Morris. (14) Although the formulas he presents seem of little practical use, 
om some of his generalizations too sweeping, he seems to have shown that 
there are at least three types of rough channel (and pipe) resistance, namely ei 
-_isolated-roughness, wake-interference, and quasi-smooth; and that itis 
foolish to suppose that an “equivalent sand grain size” with the Nikuradse 

pipe formula « om bet used to represent them all. 


Discussion of Results 74 

The results of tests are given in Table 1. It must not, 

= be assumed that values are correct to the last significant figure. 

_ Suppose that about 21 cis were flowing in the flume at a depth of 2.40 ft ‘ie a 
_ of 0.00200. This would be measured by the full-width weir with a head 
of about 0.850 ft. The measured head might be in error by 0.002 ft which . i, >. 
makes a percentage error of 0. 235%. Since Q varies as the 1.5 power of the 

head this indicates a percentage error of 0.35% in Q. But since the formula — 
which is used in computing Q is probably not correct, it seems that the 
a vs probable percentage error in Q might be as much as 1%. The slope may be 

4 assumed to! have a probable error of 0.00001 or 0. 5%. The error in deter- nal 

ia mining a. might be as much as 0.005 ft or a probable percentage error of 4 


0. 21%. f{=8SR = 90. 9s s S yn°/Q?, the probable percentage 
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TRIANGULAR CHANNE 


ai But now suppose that 0.043 cfs were ‘flowing : at a depth of 0.250 ft with a a 
slope of 0.00045. With the same assumed error as before, the percentage 
error in the slope would be 2.22% and for the depth 2.00%. The flow would = a 7 
be measured by the V-notch weir with a head of about 0.200 ft and probable 7 
percentage error of 1.0%. Since the discharge varies as the 2. 5 power of the 


head, this makes | a probable sa ge ang in 1Q of 2. 50%, or allowing for 


(5 x 2.0)* 
‘This indicates that 1 we e might expect percentage er errors in f in the - range from - 
2 to 12%. When the flow is critical or shooting the error in measuring Yn 
might be even more than assumed above. Since C and n vary inversely 
the square root of f, the expected perenne errors in them would be asi 


ns When the values of f as found from the runs on the smooth flume (calles - 

“fy in what follows) were compared with the values of f for a smooth pipe with — : 
t _ the s same Reyuaids number computed by Prandtl formula 
=? log Ry f) - - 0.80 PE: 


"(called fp in what follows) the ratio Y=f Sty varied fi from 0. 778 to 2. 248. we 


\ _ Runs 1-16 a and 102-108 where ‘the flow was definitely tranquil, wit with: a maxi- 7 


_ mum Froude number (defined as =) of | 0. 0.85, gave an average W wo o 


P 0.930. In other ‘nants tin: f was only 93% as much as for a smooth pipe at ip 


= 

for tranquil flow in smooth triangular channels. The average | 
values of f computed by this formula from the observed values was 7.1%. — 
‘Since the square root of 1,071 is 1.035, this corresponds to an error of 3.5% 
in Vor n, The average n was 0.0092 and the average | discrepancy from 
was 0.00032 or 3.5%. Therefore, the logarithmic formula has no advantage 
over the Manning formula for representing these data. 


ay Turning now to the smooth channel data in which the flow was critical or a 


1. b= 1.15. This to : 
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_ that the Manning formula does not fit these data as well as the logarithmic. wl ie 
But the difference is not great, and neither formula gives satisfactory agree 


The used in the  deacribed in 
as of rectangular strips of redwood 3/16 in. thick by 3/8 in. wide tacked to the 
7 sides of the flume perpendicular t to the flow. ' They stopped | about 1/2 in. 2 in, , Short 4 
; of the bottom, so that the flume could drain dry after the day’s runs were ys a 
_ finished, but it was felt that this would have very little effect on the resistance. 
‘s The strips had been sanded on the three exposed surfaces and the 3/16 in. We a 
_ dimension did not vary more than a few thousandths of an inch. - ‘The strips — 
_ were used without any surface treatment except that they were soaked in water 
_ before being applied to avoid their buckling out between brads when the water — 
7 was turned on, Four longitudinal spacings were tested, 24 in. center to center, 
12 in., 6 in., and 3 in, This last required 3200 battens. The longitudinal = 
sonra to center in feet is represented by the » Symbol A as in ey 
In searching for a formula to represent the rough channel results it was -: 
found necessary to divide the tests into two groups - those in which the flow — 


3 noted that these Froude numbers do not include the velocity ‘distribution 
For the tranquil case three formulas jm tested. The first was Manning’s. 


rm 


Py 
; The w weighted average » discrepancy in values of n n for non 52 runs was 2, Ob. 


The same runs were competed on the of formulafor 


= 


fox shootine flow in sm riar h Is. The average percentagze 
“9 : ee the surface is so irregular that there is too much error in measuring the oa i 
x 
— 
— fi 
— 
Was definitely tranqui Froude number less than 0.70), and those in Which 
a 4 _ flow was critical or shooting (Froude number 0.70), and those in which ae had Va 
Al 
te 
a Th 
mae 
4 
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this formula € was for each run and average value oe 
spacing and the average percentage diecrepency in f computed ‘by using 


Batees spscing” rage € Averag rage % 
24 00459 ft | 2. 


The average error in was as 4, which corresponds to 


on which for our case reduces to to oe 


= +43 - 145 (4, 88 - 
data fairly well, (This indicates that the maximum value of 
-W would occur when A = 0.2049 ft = 2.46 in. or 13.1 times the batten height. el 

_ This agrees quite well with Johnson’s results. See Fig. 14 on page 558 of aa 
fe "reference (8). The formula is not expected to hold beyond the maximum.) al) a 
3 Values of f were compute for all 52 runs by this formula and the percentage a 
found, The mean was 4.53% which to 2. in 
mets 
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ny TRIANGULAR CHANNELS 531. 
Substituting 2 Rfor rg this reduces 
2 
ter 
a 
oe 


: e fact that _ decreases with the spacing seems to indicate that evenat | 

"these spacings the roughnesses are not entirely isolated, , and that especially a 
a at the 3 in. spacing there is some wake interference. Also, the fact that none 
of the Cp’s approach the given n by Morris , throws : some doubt on his” 

_ three methods are 2.05%, 2.18%, and 2.24%, none is outstanding. Because 

4 of its simplicity and familiarity Manning’s would yr be chosen, but ie 


: aa Turning now to the runs with Froude number 0.70 orn more, there are only 
three spacings, as no runs were made with a 3 in. spacing at a steep enough vi 
slope to give critical or r shooting flow. The average values of Manning’s 
and the average percentage discrepancy from the average for spacing 


2 


ght 
‘The Same runs were computed on the basis of of Nikuradse’ s rough pipe for-— 


‘for + cael flow, but the scatter was so great that the result seemed worth- ry 
less. The Nikuradse type formula has a slight edge over the Manning formula 

for shooting flow on the basis of these tests, but the discrepancies from both — _ 
_ formulas are so large that the whole question remains in an unsatisfactory a 


oe _ state, In spite of the experimental errors present in these tests, it is | 
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v al ey 5. Ce WO SHOOUNE 110 in trl Bular Chanhneis 
definitely greater than for tranquil flow for corresponding conditions (i.e. 
at the same relative roughness for rough channel flow and the same Reynolds 
umber for smooth channel ane sults found for 
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rectangular channels 1.(8) (As has been elsewhere the 
4 terminology in that paper was incorrect as Table 1 (a) included both comet 


and shooting flow and (b) was ultra-rapid flow. Ultra-rapid flow was not ob- 


tained in the tests reported here. To obtain it in a triangular channel poke 
a much higher Froude number than in a rectangular channel. 
43 
_ the “capillary rise” of water wives the sides of a rectangular flume increased | 
- greatly when the velocity in the flume became supercritical. The “fillet” at : 


i = water’s edge, normally very small, was apparently invaded feral on 
_ Which so modified the balance of forces that the “fillet” enlarged to - in, on 


phenomenon undoubtedly occurred in the present tests in 
it was noted that water from standing and transitory waves caused by rough- — 
i ness and turbulence “rode up” the sloping side and fell back with ‘considerably 
decreased momentum. In spite of the dryness of the climate at the laboratory, 
-. the flume was kept wet for a distance above the normal shore line, ‘sometimes — 


i e. in. or more above it. High speed photographs show this condition and the 
a _—_ water in the waves, Clearly, this would add to the effective wetted 


af 


“value for the average, perhaps as many as twenty readings have to be 


a ie taken. For this reason the traverses, a sample of which is shown in Fig. 6, . 
_ cannot be considered to show an accurate average velocity at each point. POM, 
They do, however, provide a fairly reliable check on the discharge and value — 


: 2 of a@ since the average of many observations is used in arriving at these = 
_ Of more general interest are pitot measurements made to obtain the be van 
fia ation of velocity perpendicular to the mid-point of the side ‘of the flume. Close 
Pee to the flume wall and in the middle of a long straight portion of the wetted a 
perimeter, the variation of velocity ought to be free of any effect from the 
Be distant water surface, or other segments of the cross-section _ . 
boundary, this relationship could be established in terms of roughness and 
flow parameters, the mutual effect of the different boundary segments could — E 
be studied by means of interaction ‘diagrams or other suitable techniques, thus 
yielding a more general means of evaluating the effect of cross-sectional | cee 
shape. Use of the hydraulic radius implicitly assumes either constant or oe a 
mutual effect of velocity gradients generated at the | different boundary seg- 
ments, ‘While this causes little difficulty in comparing cross-sections ia 
E are not too dissimilar, ‘results | are known to be grossly 3 in error + when | the ~ 
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re 
| 
at _ wetted perimeter was taken into account. Possibly this effect could be — | de 
_Pitot- tube velocity 
different batten roughnesses. Despite the fact that more than one reading d 
was seldom taken at a given point in the cross-section, the process was time-_ | 
_ When two or three readings were taken at a fixed point, the vari- 
ae 
| 
if 
&§ 
4 
a 
oh ¥ 
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— im 
4 ( ee sub-areas and the hydraulic radius, average velocity, and discharge of each | + bay 
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_ realized that the pitot indications ‘might not be reliable ane the tube was held 
_ too close to the side. Four or more readings were made at each location, ‘but — 
as can be seen from Fig. 7, the variability is so great as to suggest that © 
twenty or more m measurements would be required to establish the average 


value accurately. * ‘Logarithmic a and semi- logarithmic f plots s of the data v were 
= and empirical relations “" either type could be fitted fairly well. ‘The 


. = ent roughnesses, 2 as shown in Fig. 8. In each case the fit is fairly | good ne near — 
the wall, and diverges near the middle of the cross- -section where the effect — 
of the the other becomes appreciable. 
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“During the 1956 tests, i it was that floating | debris thrown into the 
‘channel tended to become concentrated near the centerline as it moved down- 
seemed to indicate the existence of secondary flow. Moeller- 

7 -— Hartmann(15) gives a complicated pattern of secondary currents deduced from 


the velocity distribution in a cross- -section. No instruments were available © & 


- that) would permit the relatively weak secondary currents to be measured in 
° such turbulent flow, but the relative strength of the current near the surface 
_ was roughly measured by the following scheme A large number of small 
7 Be ga of wood, , approximately 3/8 in. in diameter by 3/4 in. long were soaked 
; 4 in water long enough so that they barely floated. Then, for runs with different 


roughnesses and very nearly the same depth, particles were released | from ae 
4 rest, one by one, at the shore and their distance from shoreline noted after ee 
they had travelled 16, 32, 48 and 64 ft downstream from the  aneeed eal release. 5 


_ The results are shown in Table II and id are plotted in Fig. 9. eal ted 


Depth N Number Average off- shore distance in feet | 
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Aug. 23, 1957 12 in. 2.53 


ow 


the of this waite seem that for 
urate flow in the triangular cross-section there is a large weak secondary 
- _ spiral that is directed inward at ihe top, from the shoreline to the centerline. 
4 (Although the floating-particle measurements were all made from the north 4 

a. bank, there would presumably be another spiral, symmetrical about the _ 4 
as — by the 1956 


ie a The velocity determinations along the perpendicular to the mid-point of the ; 
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the spiral is strongest for the. 12 i in, batten but 


‘Backwater Curves 
Although the long variable- slope flur 
_ some of the cases of the backwater curves, the objective of the present ag 
- was the evaluation of losses in uniform flow, ¢ and the only data that might be " 
: useful to test backwater curve theory were those incidentally obtained while 
Berm: uniform flow. . After the practice of plotting the profiles” at the 


_ time observations were made was started very little data of any possible Sead 
value for backwater comparisons was taken, but before that some interesting — 
profiles were obtained which seemed to justify analysis as backwater curve 4 
These computations were made by Mr. K, P. Singh (20) as a master’s thesis 
_ project in the Department of Mechanics and Hydraulics of the State University 
of Towa. Fig. . 10, reproduced from Mr. Singh’ s thesis, _ shows a typical com. 
arison of observed data and theoretical backwater curves. In general, ‘Mr. Z 
Singh found that the assumption ordinarily made in computing backwater peo ee 
. curves, that the rate of friction loss is the same as it would be for cence 
_ flow ata given depth, gives very close agreement if velocity-head changes — 
: are correctly evaluated. He compared a number of recently-proposed back- 
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ae flume. The length of the flume was found to be sufficient for “é 
the determination of hormal depths up to approximately 2.5 — 
the runs classed as critical flow along with those classed as shooting, 


results may 
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ter to the velocity was given with name mean of 2. the 
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where fp is the f for flow ina smooth | pipe with D=4R the same 


nolds mi number, as given by ‘Prandtl’ formula. This also applies (but with a 


A infinite), “The authors were unable to derive a ae epee for critical 


7 oe loss for tranquil flow in the smooth channel than for the corresponding smooth 


- ing pipe. - Also contrary to to the presuppositions of the authors, it it was | found that 
__Manning’s s formula gave about as good results as a logarithmic one, both for 


Pitot- -tube traverses indicated discharges that agreed fairly well with the 


The values of a varied from 1.03to 1.09, 


_ Study of the velocity distribution for different roughnesses showed that 2 


the variation of welocity near wall could be represented by exponential 


The existence of un secondary currents with surface direction from shore 
_to centerline was inferred from measurements of the motion of nearly- eae 


_ submerged floats, and the strength of these currents compared for the dif- 


curves incidentally obtained in me normal 


National Foundation support and U.S.G.S. and institutional | coopera- 


_ Viduals w who have assisted in carrying these tests to completion, sometimes | 
under most difficult circumstances: “Richard Tice, H. J. “Koloseus, x j 
_— ‘Fiteh, J. C. I. Dooge, W. Donahue , Preston Farrell, Leigh Nason, Cornelius — 

Shih, Strelkoff, Richard Warnock, Bertrand Barnérias, PK Mohanty, 
Hau ‘Wong Ho, G. Aron, Ram Garde, Peter Ellsworth, Peter Landweber, 

Donnan, B, Lettau, Palmer Sekora, Barton Bickel, and James B. Posey. 

Emory W. Lane ne and Ralph Parshall contributed valuable advice during the - 


= and conduct of the tests, and R. E. Glover supplied useful data ob- 


a, tion have been acknowledged in the introduction. It remains to list the indi- 1 


of c ‘cross-section in 


coefficient of drag 
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a constant in Bazin’s formula; total area 
Ban's 
length of weir crest: a constant in Bazin’s formula, 
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ig. 12, Rapid flow, battens at 12-inch centers. 
shows ride-up on shore, Height of wetting not visible 
“shore, but can be seen on near shore at bottom of pict 
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of joe in ft/sec? 


P height of crest above floor. 
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= discharge in cfs. 
- hydraulic: radius in 


Reynolds number awe weir = 


of energy 


‘mean n velocity of of approach above weir event a 
mean velocity of approach | below weir crest. 
or weir = pg? 
lor weir = pgh 


Nikuradse roughness in ft. 3 
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‘PROTOTYPE STUDY OF A PLUMBING DRAINAGE SYSTEM? 

M. Dawson2 M. ASCE and E. C. Lundquist? i 


q to express the performance of a plumbing drainage system. It is also penal a 
oem that data obtained from tests on a scale | model can be correlated with | 
those the prototype in terms of the dimensionless parameters 

idl 

“The results of the study described in this paper can be grouped into two pala: 


parameters required to describe the performance ofa drainage 
system. The data obtained from experiments run on test drainage systems 
_ were put into the form of ¢ dimensionless parameters and plots were made to _ 
determine the functional relationship existing between these parameters. 


‘full- -size drainage system; the other was a half-size model of the ‘same sys x . 
tem. _ The performance curves obtained from the tests on these two i 


_ indicated that the performance of the full-size system could have beenpre- 
dicted on the basis of observations made on the model. However, it chowe be 

Note: Discussion open ail October 1, 1959. To extend the closing date one month, . 


a - part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
_ American Society of Civil Engineers, Vol. 85, No. HY 5, May, oie Os 


a The research leading to this paper was performed by D. A Gyorog o on a a 
a fellowship at the State Univ. of Iowa sponsored by the Crane Co. 
Instr. in Univ. of Wisconsin, Madison, Wis. 
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| size drainage system could be predicted from the — 
Prof. Mech. Eng. State Univ. of lowa, lowa City, Iowa, 


pointed out that only steady were in this otedy. 


4 Pie. _ A numerical illustration is presented at the end of this study which shows — 
_ how the pressure at a typical point in a proposed system « can be —— 


from data obtained from a scale model of the system, 


plumbing drainage is composed of of piping carries 


Ey eateahail foreign matter away from plumbing fixtures in a building. a Au 
os plumbing fixtures are equipped with traps which prevent odors and filth in the 
drainage system from flowing back into the open portion of the fixture. Most 

* traps a are only a few inches in n depth, and pressures in the drainage system — 
must not deviate positively or negatively from atmospheric pressure by > m2 

* amounts greater than the pressure head in the traps. If the pressure devia- = 

_ tions are greater than this, the traps may be siphoned or blown out, thus des- 

_ troying their effectiveness as barriers against transmission of odors. Vent 
4 ‘pipes opening to the atmosphere at the roof of the building and connected to 24 
_ the « drainage system are used to prevent large pressure deviations from oc- 
curring in drainage systems. If the vent pipes are of the proper size and are 
Ss to the drainage system at the proper points, the pressure through- 
out the drainage system will remain nearly _ atmospheric even when water be 
from the fixtures is being discharged. 
_ Laboratory experiments on full-size drainage systems are difficult an ; 
to conduct because of their large size. would be advantageous, 4 
_ therefore, if measurements made on a scale model could be used to predict 
and pressure conditions in a full-size system. 
aan the research to be described two simple drainage systems were — q 
structed; One a full-size system, or prototype, and the other a half-size at, 
geometrically similar scale model of the full-size system. Figure 1 shows 
q the geometrical design of both the prototype and the model drainage systems 
_--—- The prototype system employed a 4 inch diameter vent stack, andthe 


model a 2 inch diameter stack. plastic pipe was used so 


= 


Tak 


connected to pressure taps on the vent ‘stack. A constant pressure was 
: tained at A during a given test run through control of movement of the inverted 


calculate the air and water velocities through the system. All measurements — 
and observations were made only after steady flow conditions had been 
=i The Sutation. system was defined as that portion of the apparatus bounded ~ 
ag by points A, B and C as shown in Fig. 1. Point A was designated as the en- 
e trance to the vent stack of the system, In a real ul system this point would my 


ee open to the atmosphere, but for reasons to be explained later, the flow of air 
: into or out of the experimental system, and therefore the pressure at point m 


ee: js During the experiments the pressures existing at typical points in the sys- 
> tem (such as D, Eor F) for v various air and water flow rates were studied 

closely. . Under dynamically similar flow conditions in in the 1e model and 
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= ‘the p: pressure > difference between the to = 
Point in the stack at which it is wished to determine the (such | as at 
D, E, or Fin Fig. 1), Ib/ft? | 
i Va» the air velocity ; in the plumbing stack, ft/sec 

it pag The independent variables that must be considered are: 

the the water velocity in in the plumbing stack, 
D, “the « diameter of the stack, ft were ne 
the acceleration of gravity, , ft/sec? “indicate 
density o of the air, Ibs /it3 Reyn 
density of th the x water, Ibs/tt8, a ana & counter 


; 

Po, the pressure difference between the and exit of the Reynok 
plumbing stack (at points A and C in Fig. 1), lb/ft2, 

. = It will be noted that two dependent and six independent variables are in- ee 
_ volved. The theory of dimensionless analysis was employed to combine these i 
The 


«eight variables into five dimensionless parameters. | Furthermore, if no vari- 


ables other than the eight defined above are pertinent in the definition c of the the A was ( 
problem, the physical relationship existing between variables can be ex- Bz) point w 


pressed as a function of the five independent dimensionless parameters that § tank. 1 
can be obtained employing this theory. In general, | the theory of entediink- | © measur 

oe less analysis is useful for reducing the number of parameters required to ae § The ste 
express the function ‘defining the physical phenomenon, The actual form of betwee! 

_ the function must be obtained from eee or from other known facts as- movem 
~The 

and the 

pheric, 


functional exists between these five parameters. q which t 


rate of 
the inv 
7 Wat 
watch, 
if all the parameters in the have values for 
the prototype and the model. _ Therefore, if the form of the function is. deter- . 
mined experimentally with a model system, the results ‘should be applicable — jm The 
function expressed by Eq. (1) may be further simplified. If water and 
} a air are the fluids to be employed in the model, as they are in the prototype, 
a the ratio y ty will be virtually constant over the range of | conditions likely 


be encountered in either the prototype or the This j is because 


the pressure and temperature variations that are apt to occur in either are 

not of sufficient magnitude to produce appreciable change in the density of ° 
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If the stocsaiiinn of the two fluids” had been admitted as independent v varia- 
a forms of the Reynolds’ number would have appeared in the functional | ee 
relationship. 1 However, t the effects of variations of the Reynolds’ numbers ‘ae 
were neglected in these experiments. Experience in other fluid flow problems 
indicates that variations of pressure coefficients caused by small variations = 
in Reynolds’ number are usually | not great unless a critical region is is en-— = 
countered. ke, In this “study it is is argued that the viscosity e effect and hence the 
Reynolds’ number is of little importance unless the lh of the model is very 


A was designated as the entrance to the plumbing stack, The pressure at this a 
point was controlled by the magnitude of the counterweight on the inverted air o ar 
tank. The rate of air flow into the vertical stack past point A was obtained ” ‘ 


The stack exit at C was open to the atmosphere, and differences in pressure a a 


between the entrance and exit (P,-Pc) were created by controlling the rate of 
movement of the inverted air tan 

=» #§ The pressure at the entrance to a plumbing vent is ordinarily atmospheric, 
and the pressure at the drainage exit may be something other than atmos- a a 
pheric, as determined by the pressure conditions existing in the sewer to 
which the drainage system is connected. In this investigation the situation — 
Was reversed for simplicity. The drainage end of the system (C in Fig. 1) 
was open to the atmosphere, and the pressure was varied at the stack en-_ 
trance (A in ‘Fig. 2). In this manner, , the same effective pressure difference 
(P,- A-Pc) was obtained. In addition, it was more convenient to measure an 
rate of air flow through the system by measuring the rate of displacement of - 


e 

4 


ee Water flow into the system was controlled by a valve upstream from B, Pat 
and the rate of water flow was measured with a weighing tank and a stop i 
watch. Pressures at various | stack were measured 
“ih The velocities Va and Vw were tonen as the air and water flow rates ae 
respectively in cubic feet per second divided by the cross sectional area of we : ae 
the stack in square feet. These velocities, then, were only pr proportional 4 me 
the true velocities since the Grainpipe was never 


Eq. (2) « contains four variables, and it is necessary to employ two families 
of curves as shown in Figs. 2 and 3 to express the relationship between these_ ; 
variables, In these > figures the quantities | and 
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il this pressure difference would be determined by those factors governing the 


i the model fall on common curves within the accuracy of the measurements. 7 
| This provides experimental verification of similitude in performance b between a 


are plotted against for constant values. of Vv w/VDeg . These curves 


are applicable for | only one location on the vent aack (point D, Fig. 1) at which a 
the pressure is P, Curves are not shown ; for other points along the the stack, but aa ; 
they were found to be of similarform 

Figs. 2 and 3 show that test points obtained from both the pr prototype and C. 


the prototype and the model. 


Fig. 3 may be used to predict the pressure in a similar drainage system at 
a point corresponding to that for which Fig. 3 is drawn (point D in Fig. 1). In 
addition, the required venting air flow rate may be computed if the water ond ime 
rate, the pressure difference between the stack entrance and exit, andthe = 
diameter of the stack are known, ' This presumes, of course, that the system | 

for which the calculation is being carried out is geometrically similar to that 
An example will serve to clarify the use of these data. Suppose a prototype it, a 
plumbing system with a stack diameter of 6'' and geometrically similar to : 
} that used in this study were to be set up. | The problem could be that of deter- 
| mining the pressure at point D (the point for which the curves of Fig. 3are yf 
B arawn) and the air velocity required in the stack for a given water rate of aa ose 
flow and for a given pressure differential between t the entrance and ex exit of the | 


system -Pc). If the at the stack entrance, Pay is atmospheric, | 
4 


pressure Pq, The flow conditions and discharge connections to the — ystem 
downstream from point would govern the pressure PQ, 
or this example, suppose the water flow rate is such that 


322 


= 


| 

a — 
— 205 — 
433) 


= 


— 


DATA 


OLID POINTS FROM PROTOTYPE DATA 


2 


to express the of a plumbing system, The technique was applied 

to a simple system, and the performance curves obtained are not intended to | 

be of value as specific design data, However, it has been demonstrated that 

data from different size systems (so long as they are . geometrically similar) da 

can be reduced to such a form that common curves will define the periormmcs 

of both systems, _ Specifically, ‘it has been shown that four parameters are suf- a 
ficient to define such performance so long | as the difference in the size of the e 
systems is not so large that possible variations due to —— number need 7 
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RESISTANCE PROPERTIES OF SEDIMENT-LADEN STREAMS 

1 M. ASCE, and George N. 


= Experiments were made to determine the relative effects on the friction — 
factor of a flow of the suspended sediment load and the configurations which © ys: as 
form o1 on n sand beds < of streams. Results showed th that the friction factor of a a 7 
_ stream c carrying g suspended sediment is less than a comparable on one without 
| They also showed that the reduction of friction factor due 


} ‘The res resistance to flow of clear water in a channel with fixed boundaries 
B has been studied extensively and can be predicted with a satisfactory degree ‘ 


as the or r Darcy- -Weisbach have been determined and are 

given in reference works on hydraulics as a function of the texture or rough- 

ness of the channel walls and the Reynolds 

7 However, when a stream has a movable bed, and ‘sediment is being trans-- 
ported, _ the problem of determining the resistance is much more complicated = 
than in the simple case ¢ of clear water flowing in a channel with fixed walls. Ss n 
Observations have shown that the friction factors for such sediment-bearing ee 


or alluvial streams vary over appreciable ranges. Be conten thar problem has ap 


Note: Discussion open until October 1, 1959. To extend the closing date one month, bya 
written request must be filed with the Executive Secretary, ASCE. Paper ag me. 

a part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 

_ American Society of Civil Engineers, Vol. 85, No. HY 5, May, 1959, pene 

1, Prof, of Hydrs., California Inst. of Technology, Pasadena, Calif. para 

2. Specialist— ~Scientific Engr., Aviation Corp., Farming - 
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for a long time. Despite this, the sein is not discussed at any length in 
_ the literature, and textbooks are particularly silent on this point. This leaves 


all but the expert in this field with the erroneous impression that the friction 
factor of streams is constant and can be determined once for all time. The 
: lack of discussion in the literature of the resistance of alluvial streams may — 
; be explained by the fact that the problem is not well understood, and informa- 
~ tion on it is often conflicting. For instance, some evidence presented shows — 
: _ that the friction factor of a stream increases as the sediment load increases, 
other evidence supports the opposite conclusion. q 
_ The variation of roughness of sediment-laden streams is caused by two 4 
distinct processes: (1) appearance of dunes and bars on the bed which may : 
a _ increase the roughness several fold because of the additional form resistance; 
and (2) the damping effect of the suspended sediment on the turbulence in the | 
a stream tending to reduce the hydraulic roughness. _ ‘Because of these two op-— 
: ——- posing processes, there has often been confusion about the effect of sediment 
- on roughness. Indeed, it is very difficult to separate the two effects be- 
cause, for streams in which a significant amount of bed material is carried — 


_ in suspension, the two processes are always acting together, but with differing 

relative effects. For example, in a sediment-laden stream flowing at high +a 
velocity over a flat bed of sand, , the friction factor, or roughness” coefficient, 

xn _ will probably be lower than that for a clear flow over a fixed bed of compara- 

_ ble grain roughness because of the damping effect of the suspended load. ~ On- 
the other hand, at a lower velocity with light suspended load , dunes may form, 
and the net result will be an increase in the friction factor over that for ge 

fixed flat bed of the same sand roughness. 2 
d ie: This paper presents results of experiments in which an effort was made to 
separate the two processes in order to determine qualitatively their relative 
order of magnitude. _ To accomplish this, runs were made with suspended load 
and a movable bed, and then with clear water but with no change in the bed q 
configuration. This necessitated the development of a procedure for solidify- 
_ ing a natural bed configuration in place so that the clear water could be made _ 
_ to flow over exactly the same bed as the sediment-laden flow. A comparison 
c of friction factors then showed directly the damping effect of the sediment ae 
_ Separately from the effect of changing bed configuration. In the following id 
_ pages a brief review of available information on the problem is given, and 9 


Resistance of Sediment- Stream 


ae Early, Studies of the Resistance Problem 


One serious obstacle to the of the alluvial 
; "streams by the hydraulician of the latter part of the last century was the 
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— started out with the logical idea that the loss in energy head between any two lk 
stations on a stream is merely the difference in elevation between the sta- 
= tions, It was then argued that this drop, or head loss, was used up in two ~ ae 
a # ways: (1) in overcoming the hydraulic friction of the flow, and (2) intrans- | 4 n 
porting sediment. According to this concept, the more energy used in 


. eee to the Chezy formula, the velocity of a stream is propor- tf 
ee tional to the square root of the > energy slope where the energy slope is om , 
the loss of e energy to overcome friction | per unit length of stream. | Now, if a 
ae turbid stream uses part of the total available energy to transport sediment, it 
will have a smaller energy slope and hence a smaller velocity than it would a 
: have if it carried no sediment. ‘Thus, it was concluded that a clear stream | 
Le would have a larger energy slope and hence, according to the Chezy formula, -o 
would flow faster than a comparable one carrying sediment. = — 
ar In view of present knowledge of flow problems, it is clear that some of the 
ideas upon which the above conclusion was based were incorrect. To ial Ie 
‘much of ‘the energy to transport sediment comes | from the turbulence 
«A "which itself results from the expenditure of energy in ‘overcoming the friction. 
5 _ Therefore, this energy of turbulence is no longer available to overcome fric- 
- in , and the fact that ‘some of it is used to transport sediment s should not nec- . 
7 affect the dissipation of the main energy head. Also, apparently it aa a 
was not recognized that the friction, or roughness factor, of sediment- -bearing 
_ streams varies widely with load and flow ‘rate, and that because of this the 
mean velocity « can vary even 1 without a change in energy gradient or depth, . “t 
a The energy concept held by the hydraulicians of the last century, and the pa 
deductions resulting from 1 it, W were so firmly established that they influenced bans 
the thinking of many workers of that era, and some vestiges of the ideas are “ 
be found even in relatively recent literature. It appearsfromhisearly 
publications that G. K. Gilbert(1, 2) adhered firmly to this concept. This, n eS 
doubt, had an important influence on the thinking of his time because of the a “ 
=— reputation he had amply earned in his active and productive career. ef 


* Latham, in 1886,(3) reported that observations of a storm sewer four | sell 
in diameter in Croydon, England, over a period of years showed that the oe ; 
_ mean velocity for a given water depth was always greater when the water we 
clear than when it was turbid. He concluded that the retardation of the turbid 
at water was due to the load of sediment it was carrying. _ Hooker(2) p. 288, in 
bi commenting on Latham’s conclusions, noted that Gilbert also reached the __ 
sa same conclusion, | Hooker also > went on to state, “The law of conservation of 
energy will not admit any other decision in this matter, though the statement _ 
af has sometimes been made that such a retardation of velocity does not exist. 7 
One statement contradicting the belief that a stream at a given stage would © 
ae flow faster when clear than when turbid, was made by McMath, and was a 
cluded in the paper by Hooker(2) footnote p. 289. McMath(4) was apparently 


Pp. 229, expected to confirm the idea _— a clear stream flows with less re- 


loaded flow. the opposite result was found in eleven of the | sixteen 

- pairs of runs selected for the comparison, The depth measurements in these a 

; runs were apparently not as reliable as desired, and it was thought that this 
_“May have introduced a 1 systematic error into the result. Therefore, weet 
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ic placcical flume eynerimoants in codimant trancnortation Gilhert(d) 
| experiments which had the same depth and slope except that one of them car- 
ried sediment and the other did not. The bed in both cases was composed of 
a pavement of sand grains which did not move. According to the then preva- Pag 
theory 
two 
a- 
ns- 


a for the clear flows than for the loaded flows. However, it was also noted that ie 
_ the reduction in mean velocity due to sediment transport was greater for Be tind 
small load than for large, for gentle slopes than for steep, for low velocity 
than fi for high, and for large depths than for small. . As a result of these data, Lil 
- Gilbert was forced to admit that the theory was not supported and that there — 
probably some physical law which escaped his analysis. 
AY valuable contribution to the resistance problem was made by Buckley(6) — 
- in his work ot on the Nile. He made elaborate continuous observations of slope, 
mean velocity, and cross-sectional area at several stations on the river : 
7 which showed an appreciable fluctuation of the coefficient of roughness. rc A 


Typical values of these are shown in Table 
TABLE 1. Measurements at the Khannaq Nile 


rates 


Data applying to this station on both dates: 


data show ‘clearly that at given the mean of the 
on is larger on the rising than on the falling stage, while other channel charac- a 
teristics remain unchanged, The reduction in mean velocity represents an j 
increase in the Manning friction factor of 11 per cent. A further pertinent -, 
- point to be noted is that the sediment concentration, like the velocity, >. eat 
aa _ Fig. 1 shows the ‘measurements made by Buckley(6) in 1921 at Beleida wait 
ia - Station on the Nile, which is about 28 miles upstream from Cairo and about 
| miles upstream from the Delta Barrage, or dam, which diverts water to ia 
the: canals of the delta system. The variation in slope | shown is due to the ~— 
change in level of the water at the Delta Barrage required to feed the canals, 
| The value of Kutter’s roughness coefficient is seen to decrease appreciably 
as the discharge and mean velocity increase. | For example, the following ex- 
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4 in channel the flood must have been due 
. primarily to smoothing out of dunes and bars on the river bed, with the damp- 
ing effect of the : suspended sediment being secondary, be s¢ 
.: In discussing Buckley’s paper(8) p. 231, Lacey pointed out that in the re- 
port of the Indus River Commission covering the period 1906 to 1910, it was” 
= stated that the velocity a and discharge for a rising river were higher than for 
a falling river for the same stage, and the reverse was true for the water | 
is Kantlack(6) p. 265, cited evidence which disagreed with that obtained by r 
i _ Buckley. He called attention to the result of research on rivers by Swiss and 
every case as the sediment concentration increased. According to 
j this had been demonstrated prior to 1870. He also stated that in his experi- a 
"a ence in the Punjab Canal 1 system of India there was never an excess of flow 
i= the time of high s sediment load, as’ 


on the Nile and as was | other workers in India, 
6 


. eae the opening paragraph of his paper(®) Buckley also called attention to _ 


ome laboratory work by Duponchel which showed that a sediment load tended sh 


to increase the mean velocity of flow. 


re 
Studies of the Resistance protien 

+. 7 Interest in the resistance problem in | in the past two or ci decades 


been stimulated by both laboratory and field studies of the sediment 
z One of the main objectives of research in sedimentation is to develop a theory 
a which will enable one to predict the sediment load of a stream asa function z 
ia of the water discharge and the stream channel characteristics. Once the dis- es 
¢ charge is given it is necessary to calculate the flow velocity and depth, a = 
a. requires that one know the roughness coefficient of the channel. Thus ae 


- bed showed that a flow with suspended load had a smaller friction factor than 


— 
= =86©6©)—é<“‘<«séséS@Etperrimeenntts: iby Vanoni(7,8) in a flume with fixed, artificially roughened 


E a clear water flow of the same depth in the same channel. These experiments 


a rectangular | pipe showed a a decrease in 1 friction factor with increase in load 
in in some cases, but the opposite trend in others, thus duplicating the conflict- — 


4 4 


Ae os ‘Fig. 2 shows graphs of velocity profiles measured at the center of a flume 
for two flows of the same slope and depth, one with clear water and the other 
. = a mean | Sediment concentration of 15.8 grams per li liter of 0.1 mm sedi- 
- ment, It is clear that the sediment-laden flow which has the higher velocity — 
has the lower friction factor. The right-hand graph, in which the relative ant 


: Soe: up from the bed is plotted }on a logarithmic scale, ‘shows that the SOF 


making the roughness increase instead | of decrease when sediment was added ed 


need for knowledge of resistance of streams arises naturally in dealing with 


were made first with clear water and then with various small amounts of ‘Pee a 
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= = dimensional flow, profiles can be expressed by 4 
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d=0. 295 ft 
S= 0.0025 
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U, feet per sec 


Ts a flow 0.295 ft deep and 33.5 in. wide with clear water and > * 
with a heavy suspended load of 0.1 mmsand,. 


at the profile, 7, is the shear stress at the bed, , is the 
sity of the fluid, d is the flow depth, 


uy and are e the velocities 2 at levels | ‘respectively. Sin since 
To can be calculated from the measured slope and depth, one can n calculate k .: 


from a pair of points on the velocity profile by Eq. 2. The effect of the sus-_ aa 
a sediment in reducing k, as indicated by Fig. 2, isa general one which 
occurs at all times(8) in the laboratory as well as in the 
Recent experiments by Brooks(10) demonstrated clearly how the friction 
¢ factor is affected by the dunes which develop on the sand bed of a stream. He > ae ‘ 
= friction factors for dune-covered beds which were as much as seven = 


times as large as for flat sand beds of the same material. Studies by Ali and — Ep 


Albertson, (11) and Barton and Lin(12) also showed clearly the importance of 


oun Recent studies of rivers have contributed appreciably to the iil 7 


s of the resistance problems, Leopold and Maddock(13) presented detailed ob- a 
servations of streams which show that at a given discharge high sediment — 
loads are associated with ‘high | mean velocities and low water depths. ‘This 
means that high loads are also accompanied by low friction factors, as was 
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Eden(14) made a s n” as a function 


of river stage for the lower Mississippi River. He found that n varied from 
about 0.040 at low stages to about 0.030 or less at flood stages. Carey ol 


y Keller(15) have recently investigated the bed configuration of the lower aie 
Mississippi River by means of a sonic fathometer and found that sand waves 7 
were present at all times and that the largest waves were formed at the high- 


est flow rates, They concluded that these waves must be a major factor in | 
: determining the friction coefficient of the river, although they were unable to 
7 ts measure it during their investigations, een and Keller are of the opinio 


discharge. The minimum roughness coefficient obtained by him was 


ich 7 
‘showed that the Manning iver near with increase in 


size > of the bed material, go > it seemed reasonable to conclude | that dunes 
on the he bed at lower wer discharges. buy 


_ The equipment and techniques used in the experiments are ecccatialy — 
_ those described by Brooks(10) and therefore are described only briefly here. : 


The sis?) as wel were carried out by Nomicos and are described in his Be ; 
7) as well as ina the Corps of Engineers, U.S. 
i= _The flume | dailies in india studies is 40 ft Jong, 10-1/2 i inches wide, and 10 in. 
“deep. It is equipped with an axial-flow pump with a variable speed drive, Siti | 
venturi meter, transparent section in the return pipe, and electric immersion ~ 
heaters to control the temperature of the water. The open flume section is » aad 


paint. Along the top of the channel are e mounted steel rails and a movable 
instrument carriage which can measure vertical elevations relative tothe 
‘flume to the nearest 0.001 foot. The entire flume is mounted on a truss, the oid 

_ slope of which can be easily adjusted even during the course of a run, “e° 
though the channel itself is 10 in, deep, experiments were made with depths a. 


shear distribution of the flow. 


The e elevation the surface sand relative to the on the 
a flume was read with the point gage. - To get the bed elevation once the water be 
surface had been measured, the flow was stopped and the water drained out. el 
The bed was then leveled in reaches of 3 or 4 ft with a scraper carried on the qi 
instrument carriage rails. The elevation of this leveled surface was taken as 2 

e mean elevation of the bed. ‘The water depth at : any y point was taken as the it 
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difference between water gurface‘an and bed elevations and the local elevation of 
the energy grade line was calculated by adding | the he velocity head based a2 


Profiles along the flume of the water surface, bed, mean depth, silent energy ny 
grade line were plotted for each run to determine if uniform flow obtained, 
: Usually several trial runs were required before uniform conditions were es- 
tablished, _The water surface and energy grade line were not always parallel 7 
} to the flume itself , even with uniform flow. _ Therefore, the actual slopes were _ 
obtained by adding the slope of the flume and the slope of the water surface or 
energy grade line relative to the flume. 


= - ‘The discharge was measured with a venturi meter with a throat diameter 

of 3.00 in. placed in the return line, which was a 4-in. standard pipe. The 
mean velocity was calculated as the discharge divided by the mean depth times 
the flume width, Point velocities were measured with a diam. Prandtl 


a 


P 


The mean sediment discharge was determined 
~ sampling the flow in the suction well of the pump at the downstream end of Biers 
' the flume. The ‘Sampler was a brass tube | 0. 302 in. in inside ¢ diam. with a — 
half-loop at the end so that when the tube was held vertically « during sampling 
the inlet faced upward into the flow which was downward. =  — ser i 
Point concentrations, c, in the flume were determined from samples taken a 
M ‘with a 3/16 in. . outside diam, brass tube flattened at the end so the tip o opening E> 
“was in. high by .217 in. wide. The tube was bent in a right angle much 
_ like the common pitot tube, and was mounted on the instrument carriage with 7 : 
__ The velocity at the inlets to both tiger s was ds Sere same as the local 


runs was maintained at 25° + aieatnie wi 


a Since the width-to-depth atte of the flc flows in the experiments were small 

‘compared with natural streams, the effects of side-wall friction are relatively i 
more important. To reduce the laboratory results to correspond to those _ a ae 
which would obtain for much wider streams a correction is made, This ae 
lows a method presented by Johnson(19) using a procedure developed by} 


— (20,18) This yielded an hydraulic radius rp, shear velocity U,,, and 
y-We aback friction factor fp applying to the bed itself aa free _ 
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least three one-liter bottles of the sediment-water mixture. The concentra- 
sion § tions were determined by filtering, drying, and weighing the sediment collected ae 
he accomplished with immersion heaters in the return pipe which cou 
0 watts in convenient 
nd . 
tier 


al influence. “The f friction factor tion factor fp is is s defined 1 by the the Darcy- “Weisbach ' 


a S is ie alae of the energy grade line, V is the mean flow velocity in 

a é the cross section, and g is the acceleration of gravity. To calculate the _ 
mean value, f, of the friction factor, one substitutes the hydraulic radius we ‘ 


the section into Eq. (3) in in place of rp. 
Cumulative size anon curves for the two sands used in these exy experi 4 
—_ are shown in Fig. , 3, and t the size characteristics are summarized i — 
im Table 2. _ These sands were prepared from white foundry sands with well- ‘?@ 
= grains. The material was almost entirely silica, so that the mean 
< specific gravity of the grains was 2. 65. The mean settling velocity, w, of the 
“4 various sands used was determined from sedimentation diameters of sieve a 
a fractions measured by Vanoni(21) where the fractions were obtained with 7 
next of sieves with openings of adjacent sieves in the ratio” These— 
sedimentation diameters are shown in Table 3. 
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if 


_ destroyed when the bed profile was taken, When the bed was reestablished, — a 


the flow’ was as stopped and the water carefully drained off. _ The bed was then 


solidified by spraying with chemicals, - explained in detail in the following - 


| paragraph, «After t the chemicals had set, experiments w with the stable bed were 


made, first using clear water ; and then with var ying small amounts of sedi- é 


or. ment ‘added to the system. 


| TABLE 2. f Characterist of Used in 


‘Mean sed. fall 
diam. i velocity 


at 25'S. me 


rr ‘The method of ‘stabilizing the sand bed was oe by | tests of samples 7. 


: int trays | and in actual flume tests. The most satisfactory 1 method involved 


sodium aluminate and sodium silicate : as stabilizing ng chemicals. The 
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hd run with unlorm was Iirst esta isne e usua with 
of depth, water discharge, sediment discharge, and slope. 
3 After the final determination of the bed profile by the leveling i 
| flow was again started in order to regenerate the bed configuration which was — 
— 
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the stabilizing process are a — 


hours so that the moisture content is about 10 per cent by wei weight. 


aes _b) Spray the sand with a mixture of the following three materials: ee ' 


; a: 1, 68% by volume of a 2%solution of sodium aluminate, 


_--2,: 22% by volume of a solution of sodium silicate of specific gravity ¥ ' 


retained next 


7 “ae showed ed signs of saturation, ‘This required an amount of the mixture 
i i about equal to 10% by weight of the sand being — 


One of the most difficult problems in the stabilization | process was the 


sand to form a waterproof surface. 


ply an excess of chemicals which would flood the surface, and the solutions — ‘9 
were applied gently with a paint sprayer aimed parallel to the bed instead of — 
= the bed. Although the results appeared satisfactory visually, there was, — 


no method of testing whether the grain of the 


Diameter Based on Vanoni (21) for r Nevada White \Sand 


Sieve No. opening 


‘Tyler ‘pening Mean Mean Fall ral. 


opening’ sieve sedimen- velocity velocity 4 


0.088 00.096 0. 107° 
.088 0.074 081 096 0,024 
ae 250 0.074 0. 061. 0675 0.086 0. 020 0. 0220 
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STREAMS 


Outline of 


— 
‘The general plan of the may be seen Table 
‘shows some of the principal quantities measured. The experiments are a 
vided into sets and individual experiments. . The flow depth, ‘discharge, and 


| mean velocity are kept the same in each set of experiments, but the slope is - 
varied, In Sets III and IV the depth, discharge, and velocity are kept ‘constant — 


throughout, Sets I, Il, ‘and III were made with a .091 mm sand and a .148 mm 


_ Each set contains runs with (i) movable bed, (ii) stabilized bed and clear — - 
water, and (iii) stabilized bed with varying amounts of loose sand added to the - 
system where the material added is the same as that in the bed. In Runs 2, a 
6, and 8, the bed was stabilized and no sediment was being transported ai all 
flow, i,e., the flow consisted only of clear water. 
Table 4 shows the most important measured and calculated quantities 
tained in the experiments. Where blank spaces occur in the table, the items = 
were either not measured or not calculated. Measured values of the exponent 
-z in the suspended load equation (Eq. 7) are listed in Fig. 8 a 1, 3, bs 4 


7, which were those with movable 


Observations of the bed configuration and the movement of the configuration. 
_were made visually and by means of systematic photographs during the ex- a 
periments. _ Fig. 4 shows a side view of the bed for Run 1 and a plan view wall aS 
—s + where the bed configuration for Run 2 is the same as that for Run a 
except that it has been stabilized by the method described above. Fig. 5 es 
proms a similar set of pictures for Runs 3 and 4, Fig. 6 shows the side view | 
and plan view of the bed in Run 5, . which has been described as “flat” in . - 
Table 4, It will be seen that the bed is actually not entirely flat but that 
there are small ripples near the wall. The pictures of Fig. 4, 5, and 6 ius 
trate the meaning of the termes used to describe the bed condition i in the = h a 
column of 4, 


Velocity profiles 


measurements fitted law very well, be seen 


in Fig. 7 which shows the profiles for Runs 1, 3, 5, and 6. The quantity m 


ag 
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7. Measured velocity profiles at the centerline of the flume 


at Station 24 for Runs 1, 3, 5, 6. 
voaKarman constant for the centerline was determinedfrom 


where e the centerline shear ve velocity tach is given 


where Ued is the mean velocity calculated from the “measured profile at the 

- centerline and f is the friction factor for the channel obtained from Eq. (3) or 
using the hydraulic radius of f the section r instead of yp The values of f koh, 

which are listed in Table 4, consistently decrease as the concentration 

"Sediment Concentration Profiles 

The suspended load equation for distribution of AEE Oe c of pa a 


7 
= ey | | o257 [ozs 
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q 
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q 
q 
— 
a where y is the distance up from the bottom, d is the "ae _ 
at the level y =a, and z the exponent is given by & 


1989 


ot where p is the r: ratio | of the ¢ diffusion coefficients for sediment and momentum, 


4 on the longitudinal centerline of the flume. It is seen that these ‘sediment pro- 

a files follow Eq. 7, the well-established equation for distribution of suspended 
‘. load. The values of the exponent z calculated by Eq. 8 using k,j from Table 
4 p= 1, and u, = =uxcl, (the centerline value defined by Eq. 6) are very 
7 close to those appearing in Fig. 8 which are determined from the slope of the 

_ lines in the figure. It is of interest to note that the profile for Run 1 agrees — 
with the theory as 5 well as the other 1 runs i in Spite of the fact seem sit 


Table 4 it seen that the bed | fp is smaller for 


sediment - laden and clear-water flows are in Table 5. This 

a] ‘shows clearly the effect of moving sediment in reducing the friction factor, 
which in these experiments amounted to from 5 to 28 per cent. The results : 
ae shown in Tables 4 and 5 also indicate that the larger decreases in the friction 
factor are associated with larger mean concentrations C. 

Attention is called to the fact that fp - -values for the “identical” runs, 

a - Runs 8 and 8B, with clear water and stabilized bed, are not the same. Run 8 
oa a _ made once and then repeated on the following day, and the results agreed 


> “4 closely. However, when the water had been in the flume for longer than a | . 


ee total time of about 10 hours, some changes occurred, giving the results of © 
A “5 Run. 8B. To further investigate this problem, Run 8B was repeated three ae 


times, once on each of three days, drying the — - after each experiment. 
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STREAMS 


On the second ox the experiment was completed in as short a ll 
‘ possible and the pump stopped, leaving the water in the flume. After about — 
_ three hours the run was repeated. The flow characteristics and friction fac-_ 
of the flow were changed by approximately the amount as they would 
y have been if the water had been running continuously. The run for the third a 
day was a repetition of that for the first, with approximately the same results. . 
From this it was concluded that the water was in some way affecting the var af 
nih coat applied to the sand when the water remained in the flume for long 
periods. 8 is considered to give the valid results, with Run 8B being 


ry Discussion of Results 


‘The effect of the : sediment load on the friction factor may be seen in Table 4 
51 by comparing the results of clear and sediment-laden water for each set of 
‘runs. Sets I, I, and Ill were made with the same sediment, so the results 
we directly comparable. Set IV was made with a coarser material. It is 
: clear that by adding a sediment load to a flow the friction factor is decreased, 


"since in all four sets of experiments the sediment-laden flows had smaller 


Table 5, it is seen that in Set I the friction factor fp for : a flow with a con- | 
"centration of 3. 64 grams per liter was 5 5 per cent lower than in the compara 
* ble clear flow. In Series II and III the reductions in fp were 25 per cent and 
i 28 per cent with concentrations of 4.60 and 8.08 grams per liter, respectively. , 
From these figures it appears that rate of change in de- 
creases at high values of the concentration. 


‘TABLE Comparison of Friction Factors for Sediment- Laden 
Flows with Those for Clear Water over Stabilized 
Sand Beds of the Same Configuration. Ae 


0.2 at 


O- 
le 
“2 
3 d 
on 
— 
ed friction factor seems to elatlve 
— 

— 
2. 


The reduction in n fp ma may ~— attributed to alterations | of the turbulent ‘strue- 


to attain the high vabes of 0. 106 in Run 1, compared with the values for all q 
. runs which were less than 0.030 because of predominantly flat beds = 
From these data one can conclude that in flume experiments changes in fric. -"q 


tion factor due to dunes can greatly exceed changes due directly t to the a 


study of the Mississippi River, observed that at Vicksburg 


= 


_ There were several systems of dunes of varying wavelengths, the largest oe 
which had lengths of as } much 2 as 3 two miles and heights of the order of 30 ft. — 


“files showed that at low some of the longer waves disappeared but that 


- short waves of the order of 20 to 50 ft in length and 1 to 5 ft in height re- Re j 
- mained. Fi Fig. 9, reproduced from the work of Carey and Keller, shows pro- 
files at Donaldsonville gage for two stages. For the high stage which occurred 
on ‘April 10, 1956, it is seen that there are waves of considerable length. ‘For- 
the lower. stage which occurred in July, some of the longer waves are no nt 
longer present. Despite the increase in dune height with stage, Eden(14) _ 


found that the Manning « coefficient | a decreased as the stage increased. The e 


long waves actually do not contribute much to the bed roughness ; and that the 
shorter, steeper ones, which seem to predominate at the low stages, may 
a actually result ina higher roughness than the larger ones which have been 4 
observed at the highest stages. 
Data presented by Leopold and Maddock(13, 22), an example of which is _ 
“Fig . 10, for the spring flood of the Colorado River at the Grand Canyon (1941), 
show that for a given discharge the velocity is higher on the rising than on the 
falling flood means that the friction factor is lower on the rising flood, as- a 
suming insignificant changes in slope and slow rate of change of stage. This” 


= the same pattern as observed in; flume experiments w where flows with 


‘dominant factor in effecting the change. 
Leopold and Wolman(23) report that in degrading portions of the Rio —«_—ts 
_ Grande, the river channel characteristics follow a pattern s similar to that of 
the Colorado River ‘shown in Fig. 10. . In the lower reaches of the river where 
the bed was aggrading, the velocity and sediment concentration were lower on 
= ‘ the rising than on the falling flood, so that the friction factor was not highest 


on the rising flood. _ The lower friction factor is still associated _ high — 
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Effect of Bed Configuration on the Friction Facto 

—_ : "The effect of dunes and ripples on the friction —_— 

— reguiariues Only On the rising part OF the highe age. E largest OF & 
7 ij p had heights much as 34 ft. Carey and Keller(15) observed bed configu- 2 | 

_ ns of the Mississippi River near Baton Rouge by means of a sonic fatho- | a 
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Effect 

7 Bi _ factors, although it is not clear from the field data whether the damping ef- _ 
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Bed for ‘Mississippi River (Louisiana) 

for some reason the sequence of change has been 
These same authors(23) also report that the characteristics of streams in i 
sub-humid areas, such as in the eastern United States, are much more stable > 
than those of the arid western part. In these streams ‘the sediment loads are . 
much lower and the beds apparently much more stable, so that the changes in 
the bed are a much ‘smaller than in streams whose be beds are eed 


eroded, 


_ Effect of Sediment Load on the Ve ocity Profile ise 


ing The effect of suspended gultiaai al on the velocity profile o of a stream is il- 7 i 


lustrated by Fig. 2, which shows profiles for a clear flow and sediment- laden a 
- flow with the same depth and slope. It is clear that in addition to having the 
higher v velocity at ¢ any level, the sediment-laden flow also has a higher veloci- 
ty gradient. Fig. 2 shows that very near the bed the velocity of the two flows _ 

is about the same, so that at some distance from the bed the flow with the _ 


higher » velocity gradient will have the higher velocity. By differentiating the 


equation for the velocity 1, with respect to one obtains 
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STREAMS” 
"Since in the ae flows the — velocity us is same, the gradient for 
- value of y is inversely proportional to k, and since the gradient is seento 
increase with load, k must decrease with the load. _ Thus it can be seen that 
the effect of sediment on the velocity gradient can be studied by observing ke 7 
“The fact that suspended sediment increases the velocity gradient was ob- " 
served by Gilbert(1) page 229, although the full significance of this was ap- t 
parently not realized. Flume ’ experiments(7, 8) have shown that this is al 
systematic effect, and observations on the Missouri River by the Corps of Ae 
Engineers(24) showed that this effect was also present in natural streams. “uss 
The mechanism by which the sediment increases the gradient, and hence the ; 
mean velocity, has been explained in ter in terms of the ‘damping effect of of suspended . 
= To keep sediment in suspension, i.e., to prevent it from settling due to Pat - 
| gravitational force, work must be done on the sediment grains. _ This work Dore a. 
f must come from the vertical components of turbulence fluctuations and must 
result in damping of the turbulent motion. To investigate this problem, _ se — 
Einstein and Chien(25) correlated k with the ratio of the power, P, ;, to sus- 
| pend sediment to the power, Pr, to overcome hydraulic 1 resistance to the ke 
flow, using data obtained from measurements on the Missouri River and ec 3 
laboratory experiments by Ismail.(9) A similar graph was presented(8) for ro 
flume experiments made with 0.1 mm sand. - Buckley(6) and Chien, (26) in dis- ei 
cussing the effect of sediment on the flow, postulated that the main comer 
effect of the sediment on turbulence occurred near the bed where | the concen- — ‘ 
tration was highest. — ‘This appears reasonable since it is in this zone that the my 
shear rate is highest and most of the turbulence is produced. | With this idea 
‘in was correlated the ratio of the power P' to suspend 


"sediment 3 ina . thin layer n near the bed to the power Py required toovercome _ 


frictional resistance. The e power to ‘support sediment lll 
_— with unit horizontal area and depth d equal to the stream depth is © 


ae 


wd 


where y y and Ne are the specific weights of the water and sediment respec- 
tively, c is the mean concentration over the depth in weight per unit volume, 
and w is the settling velocity of the sediment. e The anes Ps toovercome — 
the friction of the prism of water is F 


where S is the slope of the energy grade line and u is the mean velocity st 


the depth. The ratio of the two powers 


- The: power x P' to support rt the sediment suspended ‘between levels y=y 


and y=y2 ina prism of unit horizontal area is, 
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where cy is the n mean sediment concentration in per unit volume 
the levels yy and The ratio of P, to is 


4 


plotted against k for flume data ‘obtained 
with wiien: sand of 0.1 mm size, where v1 and y y2 1, have been n taken as .001 ms 
and .01 times the flow depth respectively. Fig. 11b shows a graph of the _ 
same data(8) plotted according to the method of Einstein and Chien(25) i. “7 me fl 
k is plotted against P,/Py. It is seen that the correlation of the data in Fig. a 
. is better than in 11b. This tends to support the ideas of Buckley and vel 
_ Chien and may also explain why the velocity remains logarithmic ina = 
_ sediment- -laden flow, as seen in Fig. 2, even though the turbulence has been 
modified considerably. If significant ‘damping of the turbulence occurred 
° over the entire depth, one would expect a distortion in ne distribution of the | 4 
turbulence, and hence | in the shape of the velocity profile. However, ifthe 
‘principal interference with the turbulence occurs near the bed, the turbulence 
_ energy will be reduced but it may still follow the same laws in diffusing = - 
through) = en: and one might expect the sa same shape of acpanteateadh distribu- § 
Referring : again to Fig. 2, it is seen that the two flows have the same a a 
a nd hence » hydraulic radius, and the same slope. “ee The one with the higher — 


ee _ The decrease in f, shown in 1 Table 4, for the sediment -laden flows, re- 
_ sults from the damping phenomenon discussed above A comparison of the — a 
k-vaiues listed in Table 3 will show that the turbid flows have much the lower — 


values, thus that has occurred. 


‘The present experiments have shown that a sediment-laden stream acts in 
_ two ways to change its roughness coefficient: | (1) a smoothing-out of the bed 
es tends to reduce the coefficient as the velocity increases; damp- 
_ ing of the turbulence by suspended sediment tends to increase the velocity “~ 
a wanes so that the difference between the velocity near the bed and near the ~ 
7 surface increases along with the mean velocity. In the present laboratory — 
_ experiments with 0.091 mm sand, the change in the bed friction factor fp re- 
sulting from change in bed configuration alone was as much as five-fold in 
_ going from a dune-covered bed to a flat bed. On the other hand, the maximum 


reduction in fp due to the ‘Sediment | load was 5 only 28 per cent (Run Since 


“friction i. e., 5A, also had the lowest value of k as 
. from the centerline velocity profile. From these few data, it appears that the 
; 2 _ change in bed configuration has a much larger influence on the roughness co- 
L hi efficient than the damping effect. From this one would judge that the large 4 
reductions in coefficients reported by Buckley for the Nile, Eden 
i - for ‘the Mississippi, and Lindley for the Indus River, could not be the result | 
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Gunes on the bed. Judging from the flume experiments one would | expect the: 
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0.295 0.0025+-—4— pa = 


Reduction of von Karman constant ii in sediment- -laden 
ry 
dunes to diminish as the velocity increased, while, at the same time the load 34 
increased. Since both of these effects tend to reduce the roughness coeffi- Fs cl 
cient, the field observations are in agreement with flume studies and the pa) “ 2 


process: can be ‘be explained at at least qualitatively by the combination of the two 

The explanation of the observations that for the same the rough - 

} hess coefficient of a stream on a rising flood is smaller than on a falling Fo 

: flood presents some difficulties. Buckley and Lindley noted that for a given 

| flow the concentration was higher during the rising flood than during the fall- _ , 

flood, and the reduction in roughness by the effect of 
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a te be accounted for by damping effect alone, and some readjustment of the bed 
4 must be involved in order to achieve the changes observed. - During the rising 
_ stage the stream bed probably becomes smoother because the transportation - 
‘ of a high sediment load requires a lower depth and higher velocity compared | 
With the falling-stage flow for which the transport rate is usually less; how- 
ever, since the damping effect also reduces the roughness, it is not possible 
— to determine the relative magnitude of the two effects from the field data ~~ ? 
alone. . Bottom configuration data of the kind obtained by Carey and Keller, 2 
A along with stage, slope, discharge and load data, would do much to clarify _ 
field data of Buckley, Lindley, Leopold and Maddock on the 
a tion of stream characteristics, report only total sediment concentration and 
7 do not give size distribution of the load or division between wash load and = 
material load. AS seen from Fig. 11, the effect of concentration on k, and 
hence upon friction factor, depends on the settling velocity w as well as the 
concentration, so that for ‘the same S coarse material with high settling - 
— has a greater effect on k and on the reduction in friction factor. By 
" the same token, fine materials such as silt and clay have low settling veloci- 
_ ties, and hence small damping effects even for large concentrations. ‘ie 7 
 - loads in streams such as the Indus, Colorado and Rio Grande, which show a 
correlation between mean concentration and velocity (or friction factor) (6, 13) 
predominantly fine sediments which should have small damping effects, 


for the same discharge. According to laboratory results it is unlikely that all 
of this is caused by and id changes bed are 


2 probably involved. 
duction in friction coefficient as the load increases; however, it will be re- q 
~ called that Kantlack(6) reported that Swiss engineers found that the. velocity — 
- was reduced by an increase in the load, and that Latham(3) reported a simi- 
; lar behavior in a storm sewer. ‘The latter might be explained by postulating 
_ that dunes were formed when sediment was being transported, thus causing a 
4 high friction factor, and that during clear flows sediment was cleaned out of | 
the sewer, leaving the water in contact with the smooth walls, , which resulted 
1 ce in low friction at low flow. ‘The explanation. of the behavior « of the Swiss : 
im rivers may be associated with a phenomenon observed by Brooks(27) in ~) 
4 flume at low transportation r rates. rs He observed that 1 there was a tendency at at : 
a low: velocities (i.e., low transportation rates) for the friction factor to in- 
a, crease with the velocity to a maximum, and then to decrease as the venechiil 
increased further. The experiments which showed these results were made © 
in a flume 33.5 in. wide. Two similar sets of experiments made by Brooks | 
in the 10.5-in. flume used in the present experiments did not yield a maxi- _ 
‘mum for f, but showed fp, ‘decreasing as the velocity increased for all values 
- the velocity. Therefore, the explanation offered above, based on this one 7 
t of experiments, must be considered as speculative. 


al Four sets of runs were made in a flume 10.5 inches wide in an effort to _ 
effect sediment on me friction factor from 
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STREAMS 
that of the bed configuration. First, a natural run was established over a 


loose bed of fine sand , and then, after r draining the water, the bed was chemi- a 
cally solidified in its naturally formed configuration. Loose sand was re- 
moved from the system then and a flow of clear water was run over the Bg ng tg 
stabilized bed, with the usual measurements being taken. n. Following a te 
small increments of sand were added, noting the change in the friction factor — 
which resulted. From these se experiments the principal conclusions are as as ne. 

Pee Suspended sediment reduces the friction factor of a flow. In a limited 


series of laboratory experiments, the amount of this reduction — 
“2.5 The changes i in the friction factor due to va variation | in bed ene 


ae factor may change by a factor of 3 to 6 in the ? laboratory flume between | 


of velocity al at constant depths. con 


Evidence that the load reduces the friction factor 
by damping the turbulence or interfering with its production near the 

bed where the sediment concentration is the highest and the | rate of 
production is the greatest. 


, As observed by previous investigators, the von Karman constant k is 


substantially reduced by the suspended sediment. A 


| 5 In the field it is impossible to distinguish between the two above- : 
mentioned factors affecting alluvial channel roughness. However, it 
2 has been established that the reduction due to the suspended sediment _ M 
ee of importance only for streams carrying a very high suspended load ~ 
a flat bed, and is of minor importance when there are on the 


_ The research work upon which this paper is based was sponsored by the © 
Missouri River Division, Corps of Engineers, U. S. 
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_" wave erosion tests were re performed ina hydraulic w wave flume o on pi 
“materials shipped from Yakima Project, Washington. The tests were made to 


“from the with fine base material to make the tests. 
Rk The most stable | cover blanket was formed by screening the talus material p 
on a 3/4-inch screen. The minus 3/4-inch material being placed over the _ 
base material with the plus 3/4-inch material on th the surface. Motion 
_ tures were taken for a record of the 14 tests. be 
Some factors which must be considered in the of unlined and earth 
lined canals are those which may cause erosion and, consequently, develop i 
= into an extensive maintenance problem. They include: excessive tractive b 
| forces which are caused by high velocities and, consequently, result in ero-- 
sion over a large part of the wetted perimeter; and those forces caused by _ e; 
surface waves which produce erosion on the banks at the water surface. . ae 
extent and rate of erosion produced by these forces which exist in many un- g 
lined and earth lined canals designed and constructed by the Bureau of Recla- 2 PY a 
mation, depend upon the type and gradation of the soil material through which 
we The Kinnewick Main Canal, feature of the Bureau of Reclamation’s — 
‘Project in the state of Washington, was constructed in an area where the ma- a 
terial is fine and silty and has little cohesion. It is described as a endl 


fe Discussion | open until October 1, 1959. To extend the closing date one month, on 
Pp? written request must be filed with the Executive Secretary, ASCE. Paper 2021 ill . 
of the Journal of the Proceedings of of ‘the ae = 
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When dry, it is } easily crumbled into a powdery dust. Over 70 percent of the 
‘silty material passes the No. 200 United States standard screen. _ This type of 
4 material, combined with the prevalent wind in the area, provided the condi- _ 
: i tions for possible erosion on the canal banks due to waves created by the wind, 
‘The ‘prevailing direction of the wind is generally west to northwest and is in 3 
‘the same general direction as a large part of the canal. The many curves" e 
_ which were necessary on the canal, constructed in “hilly te topography, are es- Z| 
_ pecially vulnerable to erosion caused by yy surface waves. . It was decided that 4 | 
ib 7 some type of protective cover over the fine material for the Kennewick Main 
ts Canal was essential, particularly on the curves. A hydraulic model study was 
a to help determine the adequacy of the available materials to be used as 5 
cover blanket, and determine the method of would give the 


_ of from 500 to 275 cis, . The unlined and earth lined sections have 

_ a 14-foot bottom width and side slopes of 2:1. The depth varies from 7. 61 
feet to” to 5.66 feet, | corresponding to the discharges cited above. Construction a 
of the canal provides irrigation for approximately 14,500 irrigable acres s of 
new land and additional irrigation for approximately 4, 600 acres of | land in i 
‘the ; area that had been under irrigation. The design capacity of the canal will 


permit eventual irrigation of of approximately 6,100 acres of additional land. a 


a gravel wave absorber was placed to. reduce reflections from this area. vi 

_ A wave machine, Figure 3, consisted of a bulkhead hinged at the bottom of 

me the flume and moved back and forth about the hinge by an arm connected to of 
a large pulley wheel. The pulley wheel had holes drilled at various radii from 

the center, thus giving a stroke to the!wave-making bulkhead which could be a 

an varied from 2-1/8 to 11 inches. A 3/4-hp motor connected to a set of varia- _ 
ble speed pulleys completed the wave-making apparatus. With this ‘setup, it 


Kae was possible to change the frequency of waves from 6 to 315 waves per OATS 


©, 


= 


s of wave would have been ina peerrne canal where the waves would. 


pare the various test setups Ww. with this eunditien and thereby select the best 
ra - test arrangement. . With the va various setups on the wave- -producing machine, 
it was possible to obtain wave heights up to 0.75 foot from trough to crest. | 
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— Ahad its oricin in old lake heds 
— 
oe 
| 
2 
— 9 River. between the towns of rosser d Kennewick (Ficure 1) 
| 
— 
Which Was Teet wide, 0 feet deep, and approximately 24 leet in length. 
7 } _-—s make a maximum number of tests with the material which had been shipped ‘7 
the field to the Bureau’s Hydraulic Laboratory in Denver, a plywood 
a _ divider was placed on the slope, dividing it into two sections. The 1-foot- 9 
_ a aati wide area on the side of the plywood divider next to the glass observation _ | . 
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desired which w was 45 maintained constant during a test. Water 4 


_ For each test, both still and ‘motion p seniaiis a taken as an aid in re- 
. cording the data from the model. The test procedures were designed to de- — 
a termine which cover material available on the project would give the best Ling a 
i) protection for the fine silty base material and to determine which would be * 
the most efficient way to place the material to give the best resistance to 


On the project there were two materials available which could be used as © 


cover blanket for the Kennewick Main Canal. Both were equally 
nly and in adequate supply for the canal job. Both of these materials fell laggely 
the Bureau of Reclamation filter criteria for the base material, ( 1,2) 
i ‘a as shown on Figures 4 and 5. The placing of a cover blanket over finer base — 

J 3 material on a canal bank is similar to the filter problem for the foundation of 


terial. The following criteria are used ” the Bureau of Reclamation for a 


filter material be more e pervious than ‘base material 


a 
assure there will be no hydraulic pressure built up to disrupt the filter and re 


| , engineering structures to allow drainage and prevent leaching of the base ma-— 


“Adjacent structures, q 
By b. The voids of the. in-place filter aeiilinasia must be small enough to pr pre- 3 
4 vent base material particles from penetrating the filter which causes ‘logging MY 
7 c. The layer of ae protective filter must be _— thick to ‘paliabei a 


good distribution of all particle sizes throughout the filter and also to 


adequate: insulation for the base material action is 


a Results of a large series of tests resulted | in the following ioe PM 
for a graded filter. A graded filter material is defined in this paper as - : 


a _m material having a comparatively broad range of particle sizes. Graded ma- a : 
aS & may have concave, convex, S-shaped or straight line gradation curves, 3 


and may be defined as “poorly” and “well-graded” material, , depending oa 


their gradation curve shape. Uniform grain-size filters have a similar set 


at For graded filters, it is necessary that the ratio of the 50 percent grain 
- size of the filter material to the 50 percent grain size of the base material — 


¥ 
The 
grat 
— 
— USION ali Leas indi together with 1 
— e-type water surface indicator, 
| In the early tests, a type of wave which should be produced cle ¢ 
for all tests. After recording the various offi 
(0 various strokes and frequen > of fi 
by adjusting the wave machine to various s ng 
standard wave patterns were set, and the — 
was 
‘ 
of m 
Mate 
— | 
perc 
mate 
oH 16 U 
— teria 
the t 
0.008 
= 
a 
— Test 


WAVE 


50percent size F.M. 
50 percent size B. 


a ‘The ratio of the 15 perc rain size of filter material to the 15 


4 addition the following requirements for a graded fi fil sho uld be met: a 4 


The filter material should pass the 3-inch screen for minimizing parti - 


“cle segregation and bridging during placement. Also, filters must not have — ie 
4 more than 5 percent : minus No. 200 particles to prevent excessive movement — 


fines in the filter and into drainage pipes causing clogging. 
= 


ae . The gradation curves of the filter and the base material should be ap- 

proximately parallel in the range of finer sizes, ‘because the stability and 
: proper function of protective filters depend upon skewness of the gradation © al 


_ curve of the filter toward the the — giving a support to the fines in the base 


| 


filter material adjacent to the drainage Pipe should have sufficient 


‘Ther maximum size of the perforations o: or - joint openings of the drainage pipe | 


Was selected as 1/2 of the 85 oS percent grain size of the filter | material. — = i 


designing filters for for materials containing larger than 


No. size, the base manorial should be on the basis of the 


The gravel, 4, was rounded river-deposited material. ‘Twenty 
~ percent of this material was smaller than the No. 30 United States standard _ 
sieve, and 80 was smaller than the 1- 1/2- inch size. the 


aes the limits of the graded filter criteria for the base material. The other 
f 7 material available for a blanket cover was an angular talus found at the base — ; 
- of local mountain slopes. This material was obtained from two locations. me 
_ The combined material gives a gradation of 10 percent smaller than the No. | 2 
—6U United States standard sieve, and 80 percent smaller than the 1- 1/2- = 


_ Opening. . The talus material was better graded than the rounded gravel ma- . 
terial; however, it also fell outside the limits for a graded filter for 
the base material, as shown on Figure 
on A gradation of the base material, Figure 6, showed 6 percent Sein than 
i 0. 005 mm, 27 percent finer than 0. 037 mm, 74 percent finer than 0.074 mm, mm 


Wave Tests 
Gravel Moistened and Well Mixed 
The first test was conducted on the gravel material, Figure ‘six inches _ 


material | moistened and well mixed w was on the 2: 
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feet, and the wave was at 32 waves per ‘minute. This 


duced a wave having a height from trough to crest of 0.32 foot. —— 
after the wave machine was started, the gravel and sand material began to 


= and roll down the slope. The sand which had adhered to the gravel par-— = 
ticles began washing from between the particles, reducing the friction through- 


os the layer, and the force of the water flowing out of the blanket when each 4 
wave receded from the slope was enough to start the particles to slide and 
was exposed. photograph in ‘Figure 7b shows Test 1 after the test had 
been in operation for 45 minutes. Complete failure near the water surface 
lowered a few times, allowing the sand particles to wash from between the 
gravel particles, that a greater coefficient of friction could be obtained 

_ throughout the gravel. ‘This cover layer then might have been more stable. 
The fine sand particles acted as ball bearings, allowing the — — 


‘Test 2—Gravel Separated on 3/4-inch Screen 
The gravelly material used in Test 1 was dried and aeparated ona 3/4- 
4 pm. screen. The resulting analysis curve and Bureau of Reclamation filter 


had occurred. It seems possible that if the water surface were raised and q 


4 
- eriteria are shown in Figure 8. The analysis curve for the fine base materi- _ 


al, minus 3/4-inch material, and plus 3/4-inch material are shown. Filter 
criteria limits for the base material and the minus 3/4- inch material are’ 
s shown on the same graph. The graph shows that the minus 3/4-inch material » e 
_ falls roughly within the limits of the graded filter criteria for the base ma- q 
terial, except that there is very little material between the No. 4 and No. 30 : 
_ sieve sizes. The 3/4-inch material falls within the limits of the wre fl 
ter criteria for the minus 3/4-inch materia 
A 4-inch layer of the material smaller than 3/4 -inch size was placed sm 
the 2:1 slope and covered by a 4-inch layer of material larger than 3/4-inch — 
in size. The water was raised to a depth of 2.7 feet, and a wave frequency _ 
D 7 set at 32 waves per minute. . The resulting wave height was 0.32 foot. The 
BF model, ready for operation, may be seen in Figure 9a. After 1 hour’s opera- 4 
tion at these conditions, there was practically no observable damage on the — ‘a 
a cover blanket. “The wave frequency 1 was then increased to 42 per minute, with 
iy on water depth maintained the same. This increase in frequency caused a 
height to increase from 0.32 to0.43 foot. = 
bi Pe _ After operating the model for 23 hours at 42 waves per minute, there iii! 
- only a small amount of beaching. It was believed the resulting damage, shown 
on Figure 9b, was not sufficient to cause failure in a canal installation. , There 
was some ‘movement of sand in the gravel below the 3/4-inch layer. . The ‘sand 
La he appeared to move down the slope through the gravel in the area where the two 
: layers were in contact. . As may be seen by studying the photographs (Figure * 
9), the sand moved in the top 1-1/2 inches of the 4-inch layer of minus 3/4- ~ 
ast inch material. Some of the sand in the material smaller than 3/4 -inch | in size 
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Test 3- -Gravel Placed Dry and Raked 


An attempt was made to find a less enpensive of placing | cover 
material and stabilizing it with a raking or harrowing process. It was se 
_ that with harrowing the gravel material, the finer material would tend to set- 
_ tle to the lower portion of the cover and the coarser material migrate toward 
“the surface; in other words, a separation with the coarse on the 


a wave height of 0.32 foot for a short time. This created only a small ; amount i. 
of movement of the finer material. . The wave machine was then stepped up | > 
42 waves per minute, with a resulting wave height of 0. 43 foot. The depth of 
water V was again held constant at 2. ‘7 feet. _ Asc can n be | seen on n Figure 10b, the | 


water surface as the machine was first turned on, and a small amount of | 
beaching | continued ‘until the e wave machine was turned off 24 hours later. _ z. = 
_ A photograph showing the model after it had been in operation for 19 hours if? 
may be | seen in Figure 10b. The final beach slope, after the model had been ine 


14:1. This slope seemed to be 
rather stable at the end of the 24-hour period. 


Test 4— -Gravel Separated on 1-inch 
In Test 4, the gravel was separated on a 1-inch screen, and 4 inches o 


“with 4 inches of the 1- inch material. ‘The ‘resulting mony curve 
‘plotted with the base material, and filter criteria are shown in Figure 1. 1 
The model was filled to a depth of 2.7 feet, turned on, and operated at 42 4 
waves per minute for approximately 45 minutes. _ Although the surface ma- 
terial withstood the 0. 43- foot waves 8 very well, some sand conte from the 
‘Figure 12a shows a photograph of the model after 45 of 
se The area from which the sand moved in the lower layer is seen in the photo- — 
graph. At this condition, the surface material was very stable. To create a 
ff 5 more severe erosion condition, the water surface was then raised to a depth | 
of 3.2 feet. The wave machine was operated at 42 waves per minute, and the © Paes 
resulting waves increased from 0.43 foot to 0.5 foot from trough to crest. 
With this larger wave, there was some movement of the surface matetiel; ee 
however, it was not considered serious. Increased movement of the sand & 
from the material smaller than 1-inch in size was noted. 
i A study of the photograph in Figure 12b reveals that the sand had a in noi 
the lower layer to the full depth of 4 inches. It was obvious that the compara- a 
tively larger | voids in the plus 1-inch material allowed a flow of water through 
the gravel great enough to cause the fine sand to leach out through these voids. 
The larger voids would not reduce the force of the water enough by the tim ; 
reached the fine sand material in the lower layer. The fine base material 
placed under a protective blanket similar to that which was placed for Test 4 


would leach out, as shown by 
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To find the optimum screening size, tests were continued, and ii in n this test Be , 


the gravel was s separated © on a 1/2-inch screen. For Test 5, the > resulting — 

analysis curve with filter criteria and the bed material curve are shown in i 

Figure 13. Four inches of material smaller than 1/2-inch in size was placed _ 

on the 2: :1 slope 2 and covered d with 4 inches of material | larger than 1/2 -inch in 

size, Figure 14a. The: material in this test was screened while > damp, and + a 

was noticed that the sand adhered to the coarse gravel particles even after a ; 
considerable shaking in the screening machine. same condition occurred 
in Test 1. Water was again raised to a depth of 2.7 feet, and the w wave machine ~ ; 
was operated at 42 waves per minute, creating waves with a height of 0.43 foot. “a 
After a few waves had washed up the ‘slope, a slide occurred. | This slide —p- 
peared to start in the area where the water from the waves flowed out of the e if 
gravel material and continued to the toe of the slope. A photograph oe 

_the condition of the slope after the ‘slide may be seen in Figure 14b. “ie 

= Motion pictures taken of this test, including the slide, are included in be 

record of the hydraulic model tests. é From viewing the motion pictures, it 

appeared that the the sliding surface started at the surface of the plus 1/2- - 
material, progressed through the plus 1/2-inch ‘material, then through the 

minus 1/2-inch material, and ended at the board surface at the toe of the 


es = As the material had been placed in the same manner as for the other 


test, except that it was damp ar and the sand remained coated to the rock, it is 
_ believed the sand reduced the friction between rocks in the top layer and the wae - 
force of the water flowing out of the blanket when a wave down the 
_ After the slippage, the material was removed, rescreened ina dry condi-_ eto 
tion, and replaced in the same two layers, as for Test 5. The material placed Mii Md 
in this dry condition was called Test 6. 3 The model was then backfilled to ee | 
depth of 2.7 feet, the wave machine was set as before and operated for Ba Pe 
hours. No slippage occurred, and the fine material appeared stable; there _ oa # 
was little movement of sand. - Some of the surface material (larger than 1/2- alts: 
inch) moved, forming a beaching slope. This 1/2-inch and larger material 
had a tendency to move more readily than t the 3/4- -inch and larger or the a oe 
inch and larger, which had been tested previously. | ‘The photograph, “showing — 
the condition of the model after 23 hours is in 15a. 


me Tests 7 and 8 were conducted in an attempt to determine the effect of . 
‘sand when it remained coated to the rock, causing failure. For Test 7, the ae ie 
gravel material was screened while damp on the 3/4-inch screen and placed, tee 

as in Test 2. A4-inch layer of material smaller than 3/4-inch in size was 


covered by a 4-inch layer of material larger than 3/4-inch in size. The only | ies 


difference between Tests 2 and 7 was that the material in Test 7 was screened ws ie 4 


cles. _ The flume was filled to a depth of 2.7 feet, and the wave machine was 


damp, thereby | leaving a coating of sand particles on the larger gravel parti- — ; a 
“started at 42 v waves per minute , resulting in a wave height of 0. 43 foot. ee re “ 


‘me. slippage 1 was less pronounced, however, than when the ‘material was he Sd g 


"separated on the 1/2-inch screen. The coarse material aig near sur- 
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1 mi minute operation, the material then became stable. 


- ther operation produced little change. _ The test appeared to confirm the fact 
that the sand coated to large material contributed largely to the —— 

For Test 8, the was placed, as in Test 1, well mixed, but instead 
_ being in a damp condition, it was dry. An 8-inch layer was placed on the 2: 1 

_ slope, the model was filled to a depth of 2.7 feet, and the wave machine > | 

started at 42 waves per minute, creating 0.43-foot wave height. 

$inches of beaching at the water surface occurred in the first 2 minutes, and 

- after about 3- 1/2 inches of beaching had occurred, the gravel appeared to be > 

* quite stable. A photograph showing the ‘condition of the cover after 5 hours’ a 

_ operation is shown in Figure 19. The dry gravel Lean in a more stable 

4 condition than the damp gravel of of Test , eee 


i) 


Test 10— Talus Material Moistened and Well Mixed = o. 
ig A 6-inch layer of talus material moistened and well mixed was placed on a 
2:1 slope in the testing flume. The method of placing was similar to Test 1, 
in which gravel material was used. _ The flume was filled to a depth of 2. hf 
feet, and the wave machine started at 32 waves per minute. | There was —— 
ment of the finer talus material, but not as much as that observed in the wl 
a material of Test svay «pb 20 minutes of operation, the frequency of the 


Ar 


“tively. _ At the higher frequency and wave height, there was considerable ero-_— 
sion of the talus material which was washed down to form a beaching slope | of 


_ 3.66:1. The model was operated for 2-1/3 hours at 62 waves per minute. — - : 
is A photograph showing the condition of the talus material after 3 hours of * 7 


‘§§ minute, and the water was raised slowly to a depth of 3.2 feet, then oeeed : 
toa depth of 2.2 feet. Raising and lowering the water surface caused Various 
wave heights for 62 waves per minute. scour observed showed consider- 
Biswedigs material had moved down the slope, but at no point was the 2:1 slope ex- 
- posed. A photograph showing the condition of the model after raising and © 
_ lowering the water surface is shown in Figure 16b. The model was eee 
| Ses 2 - total of 5 hours. | Test 10 indicated that at the angular talus material | was - 


“Test _Talus on n 1/2- inch Screen 
_ The talus material was dried and screened on the 1/a-inch s screen, eri : 


- giving a wave height of 0.32 foot. At 32 waves per minute, there was svery ag 

little movement of the surface material. _ The material in the penton = 
smaller | than 1/2- -inch in size also appeared very stable. 
AS no damage occurred at 32 waves per minute, the frequency of the waves 
was increased to 42 and then to 62 waves per minute, and the model was Se ae 
for period, The wave height was 0. 0. ‘5 foot, 
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respectively. At the higher wave frequencies, the surface material ‘moved > 


freely, and approximately 1-inch was eroded in forming the 


slope. However, the material in the bottom layer below the 1/2- inch in size F 
- appeared very stable and very little movement occurred. Ma | 


‘Test 9—Recommended Cover—Talus Separated on 3/4 -inch Screen 


i the base material analysis curve and the Bureau of Reclamation filter — | 


ites pre The resulting analysis curve is shown in Figure 18. Before placing 
rl the talus material on the 2:1 slope, 3.6 percent moisture, by dry weight, was F 
dded to the material smaller than 3/4-inch in size to simulate field condi- 


e ie tions. Four inches of the material smaller than 3/4-inch in size was placed 
ig = the 2:1 slope and covered by 4 inches of material larger than the 3/4 -inch 


machine» was started. ‘The model \ was operated at a frequency of 32 waves 
minute, creating 0.32-foot waves for 15 minutes. the material ap- 


size, Figure 15b. The model was filled to a depth of 2.7 feet, and the wave 4 cf 
peared to be very stable, the frequency of the wave machine was increased t to 
— waves per minute. % The model at this frequency < and wave height (0.43 foot) 
was operated for 30 minutes. ‘There v was a small amount of movement of the + 
a fine material. The wave frequency was then increased to 68 waves per min- i 
with a resulting wave height of 0.5 foot. 


‘The 0.5- foot wave approaching slope r may be seen in the photograph, 


this" « in the layer to : stop after about 
re hours of operation. Rar model was operated at 68 waves per minute for ap- 


proximately 5 hours. The photograph, shown in ‘Figure 20b, gives the condi - 
“i a ‘tion of the slope at the end of the 5-hour - period. A comparison of ae 15 


20b very the appearance of the material. qq 


A 6-inch layer of the fine base “material (analysis curve Figure | 


_ material. The 18 percent moisture is the optimum moisture content =< a 
hy oe maximum compaction. The material was placed on the 2:1 slope in layers of 
Ca approximately 2 inches, and each layer was tamped to near maximum com- 4 
bg paction with a wooden tamper. . Before placing the base material on the 2: a 
board slope, it was roughened with a layer of metal lath to prevent any ten- 7 
ency for the base material to slide on the plywood board slope. 1 ee ie 
_ The water was raised slowly in the flume to a depth of 2.7 feet, Figure 21a, 
There was some settlement of the material on the slope, and near the water 
surface, cracks could be observed. The wave machine was started and 
- operated at 32 waves per minute. The wave height varied from 0. 25 to 0.40 
foot as the material 1 washed down | the ‘slope. During this « operation, complete 


failure of the base material slope occurred, ‘Immediately on starting the 


machine, the earth material began to fail. The photograph taken 

is} _ minutes after operation started may be seen in Figure 21b. The density cur- 
rent of fine material flowed down the slope to the bottom of the flume. After 

- 12 minutes of operation, the 2:1 slope was completely exposed. The model 


was operated for 1 aay: after which it was turned off and slowly drains . 
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-inch Talus Material 


by dry and compacted in thin layers with a wooden tamper. The 
talus material was dumped on the compacted base material, simulating the 
‘ie way that it would possibly be placed in the field, Figure 22b. The model was 
ig ‘slowly filled to a depth of 2.7 feet, and the wave machine was started at 32 

Waves per minute, giving a wave height of 0.32 foot. 
- Immediately after the model was started, the fine base material began F 


_ leaching er the talus cover, the material being leached out formed a 
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(a) Test 12--after_ 1 hour of operati 
at 2.7 ft. depth--32 waves per 
height 0. 26 


Test 13--ready for operation 


Base material covered with 
plus “inch h talus” 


| im. < = 


inch screen, % Four inches of talus material smaller than 3/4-inch in size was | 


ensity current down the s in the water immediately ‘above 
e talus cover and continued on to the flume floor . The density current may 
be seen in the photograph of Figure 23a which was taken 7 minutes after the 
-model was turned on. After 30 minutes of operation at 32 waves per minute, 
_ the wave machine frequency was increased to 42 waves per minute. The in- — 
_ crease in frequency increased the wave height to | 0.43 foot. At this iis 
frequency, the leaching of base material increased. 


—— photograph taken 90 minutes after operation started, shows the ee z 
a material had washed within approximately 1/2 inch of the board slope (Figure. 

‘ es) . The model was operated an additional 17-1/2 hours at 42 waves per 
a i minute. During this time, the upper part of the slope failed completely. ‘The 

was from beneath the cover. j 


Test 14— Recommended Cover 2 and Base 
: i A final test was run using the talus material as neubiathes cover for the 
: fine base material. The moisture in the base ‘material was increased to op- 
re 3 timum moisture content, and a layer 4 inches thick was placed by the tamping 


a operation explained previously. The talus material was screened on the 3/4- 


mig _ placed by dumping over the fine base material, and on top was placed 4 inches 
of talus material larger than 3/4-inch in size, Figure 24b. The analysis _ 

iy curves are shown in Figure 18. The total thickness of the test section nal . 

we 12 inches. _ The model was filled to a depth of 2.7 feet and operated for 28 
hours ata frequency of 42 waves a minute, giving a wave height of 0.43 foot. 


am As shown by Figure 25, the test section withstood the wave very} ‘well, and 
A No movement occurred 


5 
; ; = wave height; 24 hours at 2.2 foot depth, 42'\ waves per minute, 0. 42 foot a 


"wave height; 4 hours at 2.7 foot depth, 22 waves per minute, 0.6 foot wave 
height; 24 hours at 2.7 foot depth, 51 waves per minute, 0.62 foot wave height; 
_ and 24 hours at 3. 3.03 foot t depth, 22 waves per minute, 0. 6 foot wave height. 
may ‘be seen from a study of Figure 26a, | the material withstood 
operation very well. At the frequency of 51 waves per minute, there was — 
— movement of the surface material. — During the complete test, the talus 


- material below 3/4- -inch in size was very y stable, and no leaching | of the fine . 


q 6-foot high waves at 22 waves per in the most | damaging 
5 _ condition where leaching action was concerned. However, no leaching of the 


_ The section was operated to failure after the above test, the water surface 
yas raised to a depth of 3.03 feet, and the wave machine set for 45 waves per § 


th minute, with a resulting wave height of 0.75 foot. . For this: condition, the ma- 
terial failed. _ Failure started on the surface and progressed to the bed ma- 


terial. Total failure occurred in 12 minutes, as may be seen from the photo- : 


‘graph of Figure 26b. The large waves were large enough to wash the larger 
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Results Conclusions of Model Test 


% action than did the same mean size angular particles of the talus “material = 


ity in two layers, the larger material protecting the finer material. . The la 


_ the average sized yee ee on the surface layer, the more stable was the —_ 


surface 1 was to movement t by wave action 
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we The material that fell within the Bureau of Reclamation filter criteria 


als by wave action. The study indicated that a protective material must be E 


a 3 enough to prevent leaching of the base material through it and that par- =a 


ticles on the exposed surface of the protective . material must be of vecemeal = 
_ size and weight to resist rolling down the slope from wave action. _ 
an The study further indicated that the angular talus particles form a more > I, 
- stable protection blanket than the rounded gravel particles. tte 
es _ The results of Tests 9 and 14 showed the best protective cover was was the = 
talus 1 material separated on a 3/4-inch screen and placed in 2 layers: (1) 4 q 
layer | of talus material smaller than 3/4 inch placed over the fine base ma 
a and (2) a layer of talus material larger than 3/4 inch placed on the — 
surface. The resulting analysis | curves in the Bureau of Reclamation Siter 


criteria are shown in ene 18. The test showed that the layers | of cover 


ae 


7 
model tests were conducted in the Hydraulic Laboratory 
ni of | the Bureau of Reclamation, Denver, Colorado. Most of the tests were made a 


s 6 6by P. F. Enger of the Sediment Investigations Group. The tests were per- 


formed for the Diversion Works of Bureau 


e Use of Laboratory Tests to Develop Criteria For 
Filters” by K. P. . Karpoff, Proceedings ASTM, Volume 55, 1955. _ Also is- — 


a sued as Bureau of Reclamation Earth Laboratory Report No. EM-425. pt 
2. “Soil Mechanics in Engineering Practice” by Karl Terzaghi and Ralph “a 
Peck, John Wiley and Sons, Incorporated, 1948. 
“Land Drainage and Flood Protection” 1” by Bernard A 
“Ground Water, Its Development, Uses and Conservation” by E. W. 
Bennison, Edward E. Johnson Company, Incorporated, 1947, 


. “Flood- Erosion Protection Highway Fills” by Cc. J. ‘Posey, 
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DRAUL Ics DIVISION 


| HYDRAULIC MODELS OF THE § ST. LAWRENCE POWER PROJECT 


features of the St. Lawrence Power Project and in the development of a plan 
for controlling the river construction. A description is given of the 


In the design of the of the St. ‘Lawrence 
Power Project, a fundamental part was played by hydraulic models, In all, 
s: eleven models were constructed for the Power Project, some being compre- 
. hensive river models and others reproducing complete structures or parts he 
wi structures. The models were constructed and tested at the Hydraulic Model 
_ Laboratory of the Hydro-Electric Power Commission of Ontario for the us 
- of the Commission and the Power Authority of the State of New York, who 
jointly responsible for the design and construction of the 
% In the models, the final detailed design of all channel enlargements in the pat 
International Rapids section was developed. Also evolved and demonstrated _ al 
Was a satisfactory plan of river control which would permit the construction _ i 
_ of the channel enlargements and the structures without adversely affecting -* i 
_ Tiver levels and flows, and without interfering with existing navigation. Inad- 
Ps dition, the optimum location of temporary and permanent structures was indi- Sy i 
5 cated, and assistance was given in predicting flow and velocity conditions dur- 
ing construction. hydraulic design of certain features of the structures 
_ was also developed, In the Power Project Model investigation, many interest- 


ing and complex problems were encountered, It ” the purpose of this paper oy 


| Note: “Wibsansion open until October 1, 1959. To > extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 2022 is +53 
ie part of the copyrighted Journal of the Hydraulics Division, Proceedings 8 of the the eke 
_ American Society of CivilEngineers, ‘Vol. 85, No. HY 5, ‘May, 1 1959, 
‘Hydr. The -Electric Power Commission of Ontario, Toronto, 
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Description of Project 


York, and Prescott, Ontario, downstream to Cornwall, Ontario, a distance at: 
Be 45 miles. In nature a ‘fall of approximately 92 feet occurred in this il 
7 _ distance, the main concentrations being at: the Galop Rapids, 9 feet; the Plat 
‘ Rapids at Ogden Island, 12 feet; and the Long Sault Rapids, 30 feet; the o F 
: mainder of the fall occurring at lesser concentrations and as river slope. ‘The 
| flow in the International Rapids Section is uniform as a result of 
the natural regulating effect of the Great Lakes. The long term average flow 3 
y is. ‘slightly greater than 240,000 cfs, while the maximum and minimum monthly 
mean recorded flows between 1860 and 1954 were 315,000 cfs and 166,000 cfs 
= principal yee the Power Project may be be briefly described and 


“Located at the outlet of the the side of Barnhart 
dit Island, a short distance upstream from Cornwall, the powerhouse is some 
 & 120 feet long, about 160 feet high and contains 32 generating units with a 7 
. rated capacity of about 2,500,000 hp at 81.0 feet head. Six ice sluices 
also are provided, Together with the Long Sault | Dam and 13. miles: of dykes, 
: this power dam enables the level to be raised over 80 feet at this point. — ‘Down- 


stream from the powerhouse to Polly’ 8 Gut, channel enlargements are being 


; 3 The Long Sault Dam is required to contain the power pool by closing the 2 


- for the project and must be capable of discharging the full flow of the river. 
“i e ‘Thirty sluices, each 50 feet wide, are provided, and each with a gate 52 feet 
is wide and 30 feet high. The dam is about 2,400 feet re - and is about 110 feet 


et 


— 


four  sluices or each gate. es 


quois Control D: 


main river channel south of Barnhart Island. This dam provides the spillway ; 


an Canal to permit continuance of the flow to the Aluminum Company plant ‘during ; 
construction of the project and domestic supply thereafter. The isteke con- 


ables the St. Lawrence River flows to be controlled and regulated. The dam ¥ 
is some 2, 340 feet long, 70 feet ‘high, and contains (32 sluices, each 501 
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Channel Enlargements Morrisburg to Sparrowhawk Point 
channel enlargements | are ‘required | in this reach of the river 


provide velocities not exceeding those at which an 1 ice cover may form by a 
_ packing during winter operation. © These ) enlargements also must provide 3 an 
adequate course for 27-foot navigation. 


In 

In this reach, very ‘extensive iain enlargements are required to provide 


While 1 not with the Power Project, certain navigation 
ae features in the international Rapids section are interrelated and will be © » 4 
- mentioned here. The Canadian St. Lawrence Seaway Authority has improved 
_ the channel north of Cornwall Island and portions of the channel south of ub Be 
_ Cornwall Island. The United States Seaway Development Corporation has en- 
larged the remainder of the south channel around Cornwall Island and has 
: _ constructed the Snell Lock, the Eisenhower Lock and the Wiley- -Dondero Canal 
ss upstream from the locks. Twenty-seven-foot navigation will enter the head- 4 
g a water pool by this means and after leaving the Wiley-Dondero Canal will sail 
| up the flooded river channel to pA and the enlarged river channel 


= 


pets the enlarged Chimney Point and the natural river there- 
|The: ion 
problems associated with the St. Lawrence Power 
ra ae were many and varied, and all were accentuated by the large scale of mall 
ee _ Project and the continuously large flow in the river, For many of these _ 
ae problems it was considered that the construction and testing of hedwenite — 
models would yield the best, and in some cases the only reliable method for a 
their satisfactory solution. Accordingly, it was decided to carry out a 
ik - comprehensive program of model testing to provide the maximum assurance 
that the Project would be constructed and would operate as planned. Also 


_ from past experience, it was considered probable that substantial economies 


: might be realized in view of the large scale of the project. Construction of 
first St. ‘Lawrence models commenced in 1953 and testing began later that 


a the emphasis shifted to the steps in —eeeeceel and the — of 


information useful to construction. 


their hydraulic elements, and the of 


capacity both during and after construction. 
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@) The e development 0 ofa plan of river control dur ng g the construction 
(4) The determination of velocities and currents in the navigation channels 


For the investigation of ‘the problems eleven hydraulic models were construct- 


ed, of these models were distorted-scale river models and reproduced 
altogether over 40 miles of the St. Lawrence River. The remainder were — : 
_ natural scale models, three of these being large > comprehensive models o of ail 
_ principal structures. With regard to the large : river models, it might be ob- 


served that one of the principal reasons for their construction was that inthe 
International Rapids section many islands divided the flow and abrupt changes 


in direction occurred at some critical points. It was considered that by no a) 
_ other means could reliable information be obtained under flooded conditions © 


g the flow distribution through the multiple channels, the velocity distribution — 


across the sections, the actual velocity in the navigation channels, and the de- 


To Se eae over 40 miles of the St. Lawrence River in the models, it ; 5 
a essential that accurate river bed contours and topographic data be ob- a d 
e tained for their construction, and that river levels, flow distributions a and 
currents be available for model verification. _ The collection of such data over 4 
the length of such a large river posed a major problem, 
Fortun 01 rovii 

a general coverage of subaqueous | contours and other topographic | data and» A 

_ these were used as a base. To verify the information, echo soundings were | 
i arried out at sections about one mile apart over the length of the river to be — 
reproduced in the models. Where discrepancies occurred, or in critical 
_ areas, extensive programs of additional echo soundings were carried out. In 
the Long Sault Rapids section, soundings were made from helicopters as the 


most feasible method. | Aerial photographs | of the whole reach were made to 


With regard to hydrometric information, records from several | long term 
- lock gauges were available, and these were augmented by over 30 automatic 


™ tained by establishing over 100 temporary gauges. Intensive reading of these 


a gauges at selected flows provided the necessary information, Fortunately, | 


these gauges were established at the time of the peak flow in 1952, and as the ‘ 
flow subsequently subsided a considerable range in levels was obtained. a 4 : 
a Current meterings | to determine velocities and flow distribution around sy af 
islands were made, and float. observations indicated the current pattern, 


Islington } Model ‘Testing Laboratory 
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: a verification purposes a much denser c verage was required and this was ob- ; 
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All the models for the St. Lawrence Power Project were constructed and ‘ 
tested at Ontario Hydro’s ‘Hydraulic Model Laboratory at Islington, just west 
‘of Tore, Ontario. ‘This laboratory, constructed specifically for model 
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“extended for the St. investigation. +The building is 500 feet 

_ long and 70 feet wide over above three-quarters of its length. Over the re- 

_ maining one-quarter, the width is doubled to 140 f feet, as 3 Shown on n Fig. 2, Be 

i which also indicates the plan « of the laborato: 

a of the models. The laboratory is serviced by a travelling gantry crane over r 
mon of its length, , and d accurately | levelled rails provide > vertical control = 
travelling point gauges. Two constant-level head tanks ‘provide the flow for 

the models through distribution mains; and sumps generally encircling the | Ps ? 

: building provide the return for recirculation. _ A special use of Bailey ‘bridging ; 
material provides floor apace throughout. 


Lawrence model test investigation. these models have many features in 
- common, they will ” described and discussed collectively rather than indi- 


id In Fig. 1 is shown a plan of the International Rapids motion on which is 
os 


— marked the limits of the individual models. Model No. I extends from . 
—}| Ogdensburg to Leishman’s Point and reproduces about 16 miles of the river. 
| This model contains the Galop rapids reach and the Sparrowhawk Point— a 
§ Iroquois reach where very large channel enlargements were required, and in in- 
et also the Iroquois Dam and the Iroquois 27-foot navigation lock. - Model 


No. 2 reproduces approximately 8 miles of the river from above Point Three ae 


Points to below Canada Island and contains the Ogden Island reach where im - 
"portant channel enlargements were necessary. Model No. 3 | extends from 

above Croil Island to Cornwall, a distance of about 14 miles, and contains the = 

Long Sault Dam and Power House. Model No. 3 was used for the dewatering 

studies of the Long Sault Dam, and it was considered necessary to extendthe 
ft model upstream beyond Croil Island where the river is ina single channel in ; 

order to obtain the right flow distribution in the multiple channels in the Long 

Sault area. In Fig. 3 is shown a photograph of Model No. . 3 which is typical of 7 

the various. river models ‘constructed. Model No. 9, reproduces the area from. 

above the Powerhouse to the upper end of Cornwall Island and includes the 

area where channel ¢ enlargements are necessary to obtain an €co- 

nomic head regain at the plant. ‘It will be noted that these four river models 4 

the entire Rapids reach from Ogdensburg to just above 
‘Cornwall, with the exception of an 8 mile reach in the Goose Neck | Island area ¥ 
where problems warranting model study were not considered to « exist. _ It was 
considered desirable to construct several models rather than one long ol 
“tinuous model, as work could proceed simultaneously on all ‘Without 
the tests on one reach with tests c on ot other reaches. 

The for river models 1, and :500 horizontally and 
X 100 vertically. The choice of the horizontal scale ratio was largely dictated 
iy laboratory space requirements consistent with other requirements. ate 
was made of the optimum placement of the models in the 
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ia.. Fig. en . River Model No. 3 showing Powerhouse and Dykes i in 1 Foreground = 


laboratory, and as will from | Fig. 2, a larger scale would 


have required a very considerable increase in laboratory space. This scale 
was also considered satisfactory from accuracy and workability s standpoints. 
Due to the horizontal scale selected and the large river width relative to 
its depth, distortion of the vertical scale was considered necessary for simi- 
ilar performance of model and prototype. The measurement of model veloci-— 
ties and depths was of great importance in these models, and it was desired 
F that the vertical scale be sufficiently large to enable these quantities to be © 
re measured with the required precision. % However, as the > vertical scale in 
- ereased, the difficulty of producing the required model roughness also in- 
a creased. A distortion ratio of 5:1 was selected, thus producing a vertical | 5 
; a scale of 1:100, which was considered the largest feasible e in view of the rough- 
se ness requirements, and adequate from the measurement standpoint. From — 
a past experience this distortion ratio had proven successful in other models of 
— similar nature. The scale distortion in these ‘models required we hes the 
i a Iroquois Dam and the Long Sault dam structures be distorted as well. In tests 
where distortion of these structures affected the results, their true head- ail 


as - of these structures, and the correct levels produced i in the  distorted-scale 
oo The scales chosen for Model No. 9 of the power hot use tailrace were 1:160 
horizontally and 1:80 vertically. It was considered that for this model the _ 
Seales should b should be and distortion | a dete detailed economic 
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‘Construction of River Models 


—— 


_ which has been found to be satisfactory. Over the entire model von . 
"plywood templates were € erected ona two-foot grid running both laterally and 
longitudinally. The tops of the templates w were accurately cut to the surface 
elevation of the ground or river bed. In the areas where no channel im- — 
provements were expected, sand was placed between the templates to within _ 
two inches of their tops, and then a layer of concrete was placed on top of _ - 
| sand to produce the final surface. Between templates, a ‘ “outer” on a travel- > 
ling point gauge traced the contours to the correct elevation. In areas where E> 
| channel changes were _ expected a different procedure was used. A concrete 
| “table” was poured just below the minimum grade expected. A concrete block © 
| was the then poured on top of the table with the surface finished as ny could be a 
above. _ Provision was made for lifting out these blocks so that they could be ci 
"replaced with others containing desired In Fig. 4 is shown 
photograph of the models | under construction. 


Model Ro hness 
ciel Roughness 


| bed roughness to reproduce accurately prototype levels and veloci- 

ties. Artificial roughness beyond that provided by the concrete channel beds 
was necessary due to the scale distortion of the models. _ Anticipating this 
problem, a “Roughness Model” was constructed and tested prior to the com- : 
pletion of the river models themselves. In Fig. 6 is shown a photograph of 
this model which represented, to the same scales and distortion as the river i 
models, a channel 20,000 feet long, 2,000 feet wide and with a normal depth alll ; 
40 feet. These dimensions were chosen as they were of the same order of ee: 
depth and width as the river model channels. a In this model various means of 
producing roughness were tested, including treatment of the concrete bed it- ihe 
self, the addition of various thicknesses of expanded metal screening, and the — . 
addition of one-half-inch wide vertical metal strips of various heights and at ce 
various spacings . The resulting roughness was measured intermsofthe — 
coefficient “n” in the Manning formula, and it is felt that a few of the results 
obtained may be of interest. At the normal depth, the ‘mean roughness 

| ficients obtained for the various treatments were as follows: 


Single thickness ‘metal: screening; n: = = 0. 028 
(3) Double thickness expanded metal screening; n = 0. a 
s (4) Metal strips, 8 inch centres, height = 1/2 depth; n = 
= (5) Metal strips, 6 inch centres, height = 1/2 depth; n = 
a (6) Metal strips, 4 inch centres, height = 3/4 depth; n= 
mur. ) Metal tal strips, 2 inch centres, height = 3/4 depth; n = =6 16 


F 
rom an analysis of these measurements and others obtained at various depths 


it was concluded that the metal strips would provide the best type of artificial 
roughness for models 1, 2 and 3. It was found also that the roughness pro- _ 
duced by the strips was a function of the strip spacing and the ratio of strip | e 
— to water viahaaeall As it wa was desired to have an unobstruct 


" study was desired and as a considerable change in direction of the flow oc- — | i ao 
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‘River Models enter Constructs | 
“for float 
be produced by one- -half- inch wide strips spaced o on; a six- -inch staggered “av a 
toa height above the bed of one-half the depth of mean flow conditions. Actu- q : 
ally, slots were provided in the model bed on a three-inch staggered gridto __ + 
provide flexibility in increasing or decreasing the roughness. 
_ When the models were verified it was found that the basic roughness pattern 
could only be considered as a general guide, as the “form-factor resistance”, aa 


i.e. that caused by the shape and alignment of the river channel andthe Ss ee 
‘Presence of islands, constituted an appreciable portion of the total resistance. _ ‘ 

At such points, and where rapids occurred, it was necessary to appreciably a 4 
vary the strip spacing to provide correct levels. In such cases the strips and 
“the slots on the three-inch grid proved to be a very satisfactory and flexible nt b 
means of adjusting roughness. For Model No. 9, with less distortion and fab: 

larger scales, it was found that metal screening over the ‘semesete bees was 
Sufficient to produce the desired ed roughness. 


Miniature Guevent Meter 


In all the river model tests, an accurate measurement of v velocities in the 

§ model channels was of the utmost importance as the criterion for the channel eS 
enlargement design was the production of a specified mean velocity. Similar- _ 
ly, the determination of flow distribution through multiple channels and iar 
islands upon the ability to measure velocities. To enable these ve- 
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‘Fig. 6. Model ‘Showing Metal Strip Test Pattern 
propelier-type e current meter was developed. _ The requirements for 
: _ were that it be as small as possible, offer. minimum resistance to the water, 
= register accurately model velocities between 0.1 feet per second and 2 feet — 
per second, and be remote indicating. In Fig. 7 is shown a photograph of the 
Bis current meter developed. Propellers of various pitches and having different 
ug a numbers of blades enabled velocities as low as 0.07 feet per second and as _ 
. _ high as 4.5 feet per second to be measured reliably. A feature of the meter © 
is the remote electronic velocity indicator. Previous electrical indication of 
velocity required a commutator which inevitably produced a frictional drag. _ 
This was eliminated in this meter by devising a means of utilizing changes in 
electrical resistance of the water itself around the propeller. Basically the 
_ instrument consists of an electronic meter which counts the propeller ro- P 
tations by 1 means of the changes in resistance of the water between the blades 
_ and a small electrode, marked X in Fig. 7, located within 0.05 inch of the orbit 
of the blade tips. _A second electrode, marked Y and remote from the blades, 
- incorporated to balance any changes in water conductivity. An. electronic 
ih _ circuit counts resistance changes, and the number of changes in a given time 
_ period is a multiple of the rate of rotation of the meter and a function of the 
water v velocity. Ino operation, the meter is set up in the model channel and an 
electric cable carried the impulses to the counter which is set up in a con-— 4 
venient location. The operator determines the velocity from the the electronic — 


; = and an automatic timing device. A towing channel is is provided in the 


laboratory for periodic calibration 
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‘Measurement of Surface | Stream Lines and ind Velocities 


technique been at for the e production of surface 
a stream lines from which surface velocities may be determined. _ The technique 
7 “aerial” ’ photograph is taken of the model below. . Next, in semi-darkness, 
“number of uniformly-spaced surface floats each containing a lighted candle ae 


=. 


are! released in the model channel at the upstream limit of the photograph. | es 


a ‘iT he camera shutter is simultaneously opened and the lighted floats trace — ae 
pay surface stream lines on the previously exposed photograph of the channel. . An 
attached to the camera, provides a rotating segmented disc which 
>. interrupts the field of view once every second. This produces a series of | 


_ dashes instead of a continuous stream line, the ° length ¢ of tl the dashes being pro- 
portional to the surface 


of such a 
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“The work on the over a of years, ant 


= Design of nae. 


“2 model tests, detailed designs were developed which entailed the removal of iy 

i some 64, 000, 000 cubic yards of material by land- ‘based and dredging i. 
; equipment. The criteria for the design of these enlargements v were stipulated 
in the International Joint 8 of Approval for 


‘to below Lotus Island, , designed to give a maximum mean velocity in- ss 
any cross- -section of the channel which will be used for ‘navigation not : 
pr four feet per second at any time, also between Lotus Island _ 
and Iroquois Point and from above Point Three Points to below Ogden Ee 
4 Island, designed to g give a maximum mean velocity in any cross- “section 
he exceeding two and one-quarter feet per second with the flow and at me 
i regulation of outflow and levels of Lake Ontario in accordance with Plan 2 : 
_ Of Regulation No. 12-A-9, as prepared by the International Lake Ontario % 
a Board of Engineers, dated 5 May, 1955. _ Downstream from the power- oe 
- houses channel enlargements will be carried out for the | purpose of re- ce 


ducing the tailwater level at the powerhouses. ” 
addition to these criteria, navigation interests a channel | having a 


minimum width of 600 feet with special widening treatment at bends, and with ~ 
minimum depths of 29. 5 feet ¢ and ~ 5 Sent upstream and downstream from 


from ‘backwater computations. These were then tested in the 
and the velocities compared with the criteria. A very intensive series of ae 
_ tests then were undertaken to refine these designs and test alternatives, to — 4 
oy obtain the most economical channels which would satisfy the governing cri- 4q 
_ ® teria. As the final design was approached, the critical water surface profiles 5 
corresponding to the design criteria were also obtained from the model. _ ‘The q 
design that emerged from these tests is shown in Fig. 1, where hatching indi- 
_ cates the areas of channel enlargement. ‘It is of interest to note that in the "a 
&§ oa Galop Island reach, Model No. 1, improvement of the South Galop channel was — 
- - found to be of considerable economic value in reducing velocities and exca- _ 
: vation in the main Galop enlargement which contained the navigation channel. Ts 
a Similarly at Ogden Island, Model No. 2 p proved that economies amounting to idl ry 
Many millions of dollars could be effected by diverting more flow to the south ; 
i: ee of the island by dry excavation, and thus reduce costly and difficult wet exca- 
__- vation in the Plat rapid channel to the north. In the Sparrowhawk Point - alae 
:- -Toussaints Island area it was found that the — of the river was oe 
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hat excavation somewhat in excess ss of that required by the criteria was neces- 
sary to provide satisfactory ice-packing velocities. The navigation channel 
was located in the model to conform with the requirements, and velocity — a 
measurements were made throughout to ensure compliance with the criteria. %y 
_ Float pictures indicated the direction of the currents. Comprehensive water 
- surface profile tests indicated the water levels and depths that would occur 
_ throughout the Project. In Model No. 9 a very complete model analysis indi- 


cated the economic tailrace enlargement, in | Fig. 


a In conjunction with an investigation of the foundation conditions for the 7 
Iroquois Dam, a series of tests was run in Model No. 1 to establish the a 
= location and design of the dam from the hydraulic standpoint. Tests were 
- made . first on the dam at the location and to the design originally shown on the 
FA plans approved by the International Joint Commission in its Order of Approval. 
2 _ ‘The dam as shown in these plans was located at Iroquois Point, was curved in © 
plan, and had 40 sluices each 50 feet wide, The tests indicated ese 
conditions, as an abrupt change in direction of the flow occurred at the dam e 
_ making a considerable portion of the sluice area ineffective. After testing six 


ie plans, it was concluded that much more satisfactory conditions would “an 


result with the dam moved downstream to Point Rockway and with a straight 
rather than curved alignment. At the new locations, channel enlargements _ 
7 were necessary upstream to meet the ice-packing velocity requirement. 
| . _ Much-improved hydraulic conditions were found with the dam at the down- 
stream location and with the associated channel enlargements. In Fig. 8 are 
shown stream lines with the dam in the original and final locations. - the 
abrupt change in direction of the flow at the dam has been eliminated, thus ; 
a permitting a reduction in the number of sluices from 40 to 32 while providing 
‘ an even greater discharge capacity. Many | other advantages also resulted; — 
. half the dam could be built in the “dry” at Point Rockway thereby greatly _ 


z levels at the original location; head losses in the Iroquois reach were sub- ae 

: stantially reduced resulting in a benefit of many millions of dollars in ultimate — 

7 power generation; the foundation conditions were much more favourable at the — 
“new location, which, together with the fewer number of sluices, resulted in a wy 


ls g greatly reduced dam cost and an overall saving of close to a million dollars a 
in total construction cost of d dam and channel % 


‘Diversion Channels for Long Sault ist 


" During the first ‘Stage. of construction of the Long Sault Dam, it was planned _ 
‘ a the southern portion extending from the United States mainland to Long ee 
‘Sault Island would be constructed, and the South Long Sault channel se § 


e through “Cut Ay shown in Fig. the . During the second stage, the North Long 
Sault channel would be cofferdammed just upstream from the Long Sault 
. “= and the flow diverted through “Cut F” diversion | Channel and through ~ 
. deep diversion sluices left in the first stage of the dam. The design of these = 
diversion channels was carried out in Model No. 3. The criteria for their de- 
_ Sign was that sufficient discharge capacity must be provided at all times by | 
ote channels so that at levels | upstream would be satisfactory for the ciiieence 
of 14- -foot navigation and for the diversion of water through the Massena power 
Canali. . Inthe of “Cut a very thorough | study was made to develop its 
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- most economic location and dimensions. This channel, which had to carr 
+ almost the whole St. Lawrence River flow, involved the removal of about 
B 6,000,000 cubic yards of material. In Model No. 3 studies were made of the - 
. sequence and timing of the complete series of construction operations re-_ 
quired in the transitions from one stage into another, and a procedure w al 


d 
Control ed 


it was required that the Project be constructed without 
fecting either upstream or downstream interests and that existing 14-foot 
navigation continue without interruption throughout the construction period. “ 
= In effect this required the general maintenance of pre-project levels and out- L 
flows from Lake Ontario and pre-project levels in the International Rapids 
. section at the locks and in the river used by navigation. These ptr 
_—- a considerable problem in developing an overall plan which would 
_ coordinate the construction of the channel enlargements and structures to pro- a 
duce these conditions. In the development of this plan, the three large river ee : 
models, Nos. 2 and 3, were used. The large natural-scale models of the 
Bs Iroquois and Long Sault Dams, Models 8 and 5, were also used in determining 5: 
+ detail the effects of construction at these dams. Utilizing all these — 
a step-by-step chronological program of construction was carried out in the 
= models and the effects of the various operations noted. In the earlier stages, 
the timing of various operations was determined such that the effect of one — = 


_ operation would be balanced by the effect of another. In the later stages, the 


control required at structures was detailed. coordinated was —— il 


The Sault Dam Models 


left generally to river bed elevation in each sluice to pass the full 1 river r flow 
_ in the second stage. The permanent sluice gates were adapted to act as di- 
| Pd ersion gates in these openings for control during the second stage. In the | : 
= controlled by a gate, to pass the full river flow during the third stage. : 
res In the third stage the ogee sections were completed in the thirteen | southern — 
_ Sluices. _ Following this operation, the port gates were closed to raise the — ty 
“i water to ‘the operating level, and the flow controlled by the permanent sl sluice 
In the hydraulic design of the dam, several protoms arose 


construction and ultimate operation periods, and operation period 
a required, the design o of an adequate energy dissipator. oar % 
’ (2) The hydraulic loading on all the gates to ensure adequate > hoisting c we 
The discharge capacity o of the and the 


q — 
at 
Sault Dam was constructed in three stages. During the first 
stave hat nortio ntaininoe the southern el 
at 
a 
— 
4 
The need for energy dissipation downstream from the dam during — 


cription of Mo 
P he Models were built to investigate these various problems. . The first, 


Model No. 10a, was a model of the dam consisting of two piers locat- 


in construction, i.e. diversion sluices, closure ports, or final ogee profile. * 
This model, shown in Fig. 9, was constructed of plexiglass to a natural 1:50 
scale and tested in a plexiglass sided flume, 30 inches square and 26 feet long. 


_ Way was constructed in such a manner that it could be converted to all stages 4 


ae pator design developed, and exact models of the various gates were made and — 
; tested. The results of these gate tests are reported in a later section. are 


parts of Long Sault and Barnhart Islands, and a portion of the United States — 


mainland. The river channel reproduced extended some 1,000 feet encod 
: i and 2,500 feet downstream from the dam. - The model was constructed entirely 


in concrete and the dam was made with removable sections such that all 


‘stages in construction could be reproduced. The energy dissipator Acibapet: 


in the sectional model was also constructed downstream from the dam. A 


photograph of the model is shown in Fig. 10. The primary purpose of this 2 
a model was to check the efficacy of the energy dissipator developed in Model 10 3 


over the full range of conditions ina comprehensive setting. 


In this model the energy dissipation problem was investigated and the dissi- a 


a8 The second model, No. 5, was a comprehensive model of the Long ‘Sault 
Dam and reproduced, at a natural scale of 1:80, the entire Long Sault Dam, & 
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past experience, it had been found that in many cases asbeneen similarity sath. 
_ cannot be achieved in in sectional models due to the restricted channel repro-— 


— 
= 


of the diversion sluices and ports seers were affected by high velocity of ap- 


From the model investigation a very severe energy dissipation problem sas 


| was found to exist at the Long Sault dam both during construction and under | 
final operating conditions. Because of the nature of the river bed, which aan 
slopes downward rather steeply downstream from the dam, the tailwater depth a 
found to be low under all operating conditions, and insufficient to forma 
F hydraulic jump. Preliminary tests indicated that, without energy dissipating 4 


works, the lack of adequate tailwater depths resulted in downstream velocities 
5 which would endanger the safety of the structure even with the relatively low _ 


_headwater levels that would exist during the construction periods. The design 
was unique in that the dissipating ‘works must provide protection for ra 
three entirely different flow conditions which would be encountered during “ie 
construction and under final operation of the dam. 
i After many trials in the plexiglass- -sided flume, where over thirty different 
wae were tested, the Recommended Design was developed, and was veri- ae 
fied in the comprehensive Model No. , 5. It was found that this dissipator main- 
tained a forced hydraulic jump un under all flow conditions, that 


of the order of 75 feet per second with no dissipator were reduced le Bag 
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Fig. 10. Comprehensive Model of Long SaultDam — 


the three flow conditions under v ‘which the dissipator had 
ca _ Although the Recommended Design was the best developed, rock seinem 
: downstream from some of the sluices made it advisable to raise the bottom of 


ie ‘the bucket by at least 5 feet. To meet this ; requirement the Revised Recom- 
mended Design was and was shown to be almost effective as the 


4 


almost half of the dam in the “dry” on the United States shore. During this ot 
period, , the channel enlargements on that shore were also completed in the ing 
“dry”, a rim dyke extending a short distance out into the river contained the _ 
a construction area. When this phase was completed, the rim dyke was pro- ae 2003 
gressively removed and the river diverted through the completed portion 
a ‘ ‘the dam. Simultaneously, cofferdams were progressively extended out into 
the river to dewater the area for the second half of the dam. This work had 
be carefully planned to maintain pre-project levels upstream. 


«ttt It was decided that the following problems during construction and oper-_ 


capacity during each stage of construction. Also required wasa 

—— of the velocities and currents that would e exist at the ‘coffer- 

— CC particularly during the construction of the second stage coffer-_ 
(2) The | determination of energy dissipation requirements just downstream 

from the dam. during the later stages of construction, when headwater 

 — would be raised and tailwater levels would be at or below natural 


In the final location or the Dam, it was build 


3) The determination of the hydraulic forces acting on tl nthe sluice gates, to . 


oe (4) The development of a a discharge calibration of the cial for construction 


and operation conditions. The pattern of gate to 
most satisfactory 1 river conditions was 


= % 
models were c constructed of the Iroquois ‘Dam the «foregoing 
problems. _ The first was a sectional model constructed of plexiglass to a 
natural scale of 1:50 consisting « of two piers located to form a central Mm. 
- sluice with a half sluice on each side. This model, designated No. 10b, was Ot 
_ installed and tested in the testing flume, and was used to investigate energy 
dissipation requirements and to design the energy dissipator. Exact models — 
of the sluice gates were also built and tested in the model to determine the i 
_ hydraulic loading on the gates. The results of this .emte loading t test are re- 
The second model, Model No. was a large 
}. ; ‘scale model of the Iroquois Dam, which included the whole dam and itsas- _ 
: sociated channel enlargements, and the downstream 1 approach to the Iroquois: 


Lock. The model reproduced the r river ‘upstream for some 5,300 feet a and Be: 
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downstream 1m about 3, 300 river the mo model was constructed 

entirely in concrete, while the dam itself was built of plexiglass. The coffer- 

o cells for the second stage cofferdam were also built for construction coe, 
in and the cofferdams for the second stage. 
-” As in the case of the Long Sault Dam, it was desired to check the per- — 

formance of the | energy dissipator ina 4 comprehensive setting. Also, the mm = 
model was required to determine in detail the velocities and currents to be a 
= during construction of the  cofferdams during the « dewatering of the 
_ second stage area. ¥ The discharge calibration of the dam had also to be ob- 

_ tained from this model; as under final operating conditions the a across" 


the very 1 and approach conditions affect the results. 


TO 


_ Under final operating conditions, the water-level differential across the Mee 
"Iroquois Dam is quite small, and downstream velocities are low. However = 
- during one period prior to raising the levels in the forebay of the Steet 
Island powerhouse, the river was discharged through the sluices of the first - _ 
“half of the dam with headwater levels raised and with tailwater levels at 
_ slightly below natural. From model tests of these conditions, downstream cil 
locities as ‘high as 45 feet per second were ‘measured, and the need for energy 
"dissipation for this period was demonstrated. _ Several schemes were investi- 
_ gated in the flume and the Revised Recommended ‘Design was developed. By 
_ its action, velocities downstream from the apron of the dam were reduced ua 
approximately 17 feet per second. Additional tests indicated that the pro- __ 


tection of ten sluices would be adequate for this period, and the dissipator was 


‘omprehensive Model of Iroquois Dam 
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_ constructed to protect this number of sluices. Details of the design developed 


A 


From Model No. 8 a complete was obtained for the 


dam during all stages of construction. . This model, however, proved almost 
indispensable in predicting the optimum gate opening pattern to be usedto 
produce satisfactory conditions | downstream. This was particularly so in the 
period of construction when the headwater level was raised, but tailwater __ 
levels were still natural. - Under these high head conditions, up to 15 feet, the 
ad issuing velocities were necessarily high and unavoidably affected the veloci- 


_ ties and currents downstream in the approaches to the Iroquois Lock. iF From 3 
= the model, however, a pattern of part- gate openings was devised which woe & 


aa was ; planned that diversion of flow through the Massena power canal 7 

- would be continued during the construction period and also for a limited period 
after the raising of the St. Lawrence power pool. In the design of the intake 
‘it was necessary to ensure that there was sufficient gate capacity at normal 
river levels, and to determine the performance of the intake structure under | 

_ raised water level conditions. It was desired to investigate energy dissipation 
and to determine the hydraulic loadings on the intake gates under high-head 


= of Model 


and a portion of the. existing power canal to ensure correct approach and exi 
Conditions, _ The intake was constructed of plexiglass and the model bed of © 
concrete. Models of i service gates were also constructed, 
_ From the model, the best alignment of the structure in relation to the 
channels was determined to minimize scour in the downstream channel. . 
Limited energy dissipation was found necessary, and this was achieved by z 
ramping up the floor of the sluices near their outlet to direct the jet toward 
the surface, and directing the jet by an extension of one pier. Interesting 7 
sults were aia 3 from the gate tests and these are discussed ina later — 
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; —@) The best location for the piezometer | taps for the Gibson n method of flow 


.™ The hydraulic design of the ice sluices and the calibration of the ice 


_ Model of Powerhouse | Water Passages 


house passage from headwater to tailwater. ‘This model, to a 1:36 

natural scale, was constructed entirely of plexiglass to enable observations to to a 
be made of internal flow conditions. Exact models « of the headgates were made i 
of sheet metal and installed in the model, and were used to determine the nail 
vertical hydraulic loading on the gates. The results of the headgate tests are 
given in a later section. Fig. 14 is a photograph of the model which shows the 
considerable number of Gibson test piezometer connections tested. At this aan 
powerhouse three large intakes —s each turbine, and the objective of the — 
and the best locations for the ‘connections in each c cross- -section 
to record the true average pressure. The model was found to be most useful = 


bd 
4 
4 
‘oa 


May, 


Model of Powerhouse Siuic es 


At the powerhouse, le ic a 
downstream. Two adjacent sluices, 75-foot wide, are located at 
pe end of the powerhouse; while two adjacent sluices, 50-foot wide, are located | 
at the middle of the plant. Each sluice is controlled by a drum-type gate. ee 
"When operating, these sluices discharge ice and water through the powerhouse 
from headwater to tailwater in a closed water passage which does not run full. 
sg mu To. study the behaviour of ice being discharged through these sluices, and to - 
: develop a design which would operate safely, Model No. 4) was constructed of iS 
_ plexiglass to a natural scale of 1:36. Paraffin wax was used to represent ice — 
and Fig. 15 ‘shows a a photograph of the model in action. Earlier designs t tended 
to trap the ice inside the passage at the foot of the slope with the ice churning i 
heavily against the structure. However, by adjusting the sluice profile and a mm 
selecting the right floor and ceiling elevation in the lower portion, it was 84 1 
_ found possible under all conditions to discharge the ice and water at high ve~ 
: locity through the ‘passage an and i into the saroapbee without the the hydraulic a 


made by the ¢ gate ‘manufacturer. a calibration 
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the gates, Model No. 11 was eneneeand of the upper portions of the ic ice » chute 
“to a natural scale of 1:36. Two 75- foot wide sluices and two 50-foot wide - 
were > reproduced and tl and the discharge calibration n obtained. — 
a of the model described earlier it was mentioned 
in tests were run to determine the vertical hydraulic loading on the water- - 
controlling gates of the various structures. It is considered desirable 
‘port: on all these results in this section, as it is of interest to compare the 


results obtained. In all seven different project gates were tested to determine | = 


operating conditions. of the gates were t to be completely submerged ir in 
operation while the remainder were not. The submerged gates were the | * 


intake service and gates. Unsubmerged gates tested, were for the 
Iroquois Dam and for the diversion sluices and spillways of the Long Sault _ ay 

Dam. ‘Because of the interesting nature of t the results, a summary of the test We _ 
findings has been included in this paper. . It should be noted that in all the test 
results the effects of buoyancy, gate weight and friction in the checks a 
been eliminated, and the loads shown are therefore due to vertical hydraulic — OM 
action | only. ‘The apparatus \ used in the gate tests is shown in Fig. 9, where or 

| the gate is angyennes from a hoist through a spring balance used to measure 

_ the vertical loading. . Friction between the gate and the checks is eliminated ese 

4 by pulling the gate upstream just free of the checks. This is accomplished . a 
a epecial harness which a horizontal | pull and also that the gate is 


Perhaps the teats were on the temporary gates, 
_used to close the tunnel ports in the Long Sault Dam. _ While these gates were 
the smallest of those tested, their weight being about 25 tons each, it was oe 
found that they were subject to the most severe loading | conditions. As shown 
| in Fig. 16, investigations on the preliminary design indicated that a maximum 
“load of 154 tons could be expected if the holes required for eventual ae 
_ of the ports were plugged, and 124 tons if the holes were left open. By alter- 
ing the shape of the bottom section to 45° and then to 60° as shown in Fig. 16, 


it was to reduce the maximum to 73 t tons and 54 = 

| 


are shown in are where the grouping iesodh upon whether the gates ch 
submerged or not. While it is not practicable to analyse the —— detail, “4 
unsubmerged gates conventionally ‘installed with skin plates, 
such has at t the Iroquois Dam and Long Sault Dam spillways, it would appear Ps 
that neither downpull or uplift is a significant force. The long Sault Dam di- % 
version sluices, being an adaption of the final sluice gates, were not con- aa ah 
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that of the gate, with ‘the result that water collected on on 1 the _— down- 
stream members and thus increased the indicated downpull. 
oa three submerged gates were generally similar in size and illustrate 
the effect of some of the varying factors. At Massena the skin plate was on 4 
Ses side of the emergency gate and on the upstream side of the — 
service gate. The large hydraulic downpull on the emergency gate, and the — 


the skin plate influences greatly the hydraulic loading. . These loads may iow | 
_ a considerable bearing on the required hoisting capacity or the ballast re- 
; quired to ensure that the gate will close. The powerhouse headgate was tested 
under turbine r runaway ‘conditions with the other two | gates in the supply 


passages closed, and the large resulting downpull shows the effect of the | in- = 


crease in head across the oy 


a 
HYDRAULIC. ACTION ONLY 
No WITH BUOYANCY EXCLUDED 
| AIR_ REACHES DOWNSTREAM 
SIDE OF GATE AT POINT X 


20 40 60 80 100. i20 140" 
“HYDRAULIC “DOWNPULL ~ ¥ 


| 
Zz 
a 
— 
— 
: 
MEASURED VERTICAL HYOR DAM; 
— 
DIVERSION PORT GATES. AT 


> 


MASSENA 


DIVERS 


INTAKE 


“eg 
20° 


20 
UPLIF 


MERGENCY” 


100 


a 


40 


ud 


” 
Ww 


‘TO 


TESTS 


0 


RESULTS 


‘VERTICAL 


T 


GA 


LOADS 


AKE 


LIC 


HYDRA 


- 


‘LONG 


‘DAM AND 


1G. 17 


‘IROQUOIS 


INT 


= 


@-- 


= 


4 
May, 1959 


the review, it is apparent that a very 
he been made by the Power Authority of the State of New York and the ‘Hydro- si 
Electric Power Commission of Ontario to test thoroughly all aspects in the — 
hydraulic design of the St. Lawrence Power Project. It is considered that the 
= have proven | to be e indispensable ir in providing designs | which carry the =| 
maximum assurance of successful performance and at the same time make é 
possible maximum economies. In the case of this Project, economies were 
effected which exceeded by many times th the cost 0 of the models. — Information Le a 
= the models enabled each step in construction to be taken with confidence, 
in the knowledge that its effect had been tested in advance and found to be 
within the criteria. _ Another considerable advantage of the models, » was th the e 
ability to demonstrate in advance | proposed plans and procedures to super- _ — 
visory boards and other interested agencies. In it is felt that the models 


have played an essential part in the ect. 
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 tItis apparent that the carrying out of a comprehensive m ‘model 
ca testing program for a project such as the st . Lawrence requires the full co- ei 
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of ‘Mr. D. G. Harkness, supervised the program, and the able and re- 
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Clark’ s discussion yn of his pa paper. 

example between the instantaneous hydrograph and the template routing, are 7 
realistic and accordingly lead to a better understanding of flood-wave phe- 7 a 

nomena. His caution that lag is in reality translatory storage is well taken. a aes 
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DALRYMPLE,! la. M. ASCE. ~The author wishes to thank 


‘Turner for his careful review of the paper. The discussion presents con lacy 
eae material that adds to the subject of the paper, and thereby improves -  &§ 

i= _ ‘The author would hesitate to assign a frequency to a given flood that ne- I 

_ cessitated a long extrapolation of knowledge. Generally, a 50 year recurrence 

_ interval is about as far as is justified, except in a few instances where long “o-. =: 

_ term knowledge is available, as in New England. ; & project now being com- 

- pleted by the Geological Survey has uncovered considerable historical flood 
information in New England, and promises to form a basis for preeeppeied flood pe 
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Discussions L Alterman and N. 


ore 
ALTERMAN. 1_This paper paper is a welcome addition to the but 
paneer literature o on the subject, as evidenced by the large number of Refer- 
The addition of semi- empirical data of horizontal forces on barge ‘models 
is a timely and very important contribution to problems arising in the —- 
_ increasing number of sea- ~structures executed. T The a article is partly of inter- 
est for those who prefer a theoretical approach and partly for those who prefer *3 
‘For compensation, it might be of interest that the writer has made use ~ ~ 


the theoretical formulae, such as (6), (7) of the article, or rather r an inte- an 


_ gration of those, and the results have shown a dramatical application: Pee 
es series of 320 steel piles, 16" diameter, driven 33 feet into the sea- -bed i 
P ed 


wwe 


have served as a temporary jetty. Once the superstructures had to be dis- 
? mantled, the writer calculated the maximum waves which these pies could — 2 
_ withstand freely, as winter storms approached. 


| | = The computation showed that 16-ft. waves, with a period of 15 seconds =a 


would cause breaking stresses to the piles. "Orders were given accordingly — 

to cut-off the » piles as quickly as feasible, before the stormy period would 

_ Before all piles could be cut off, a first wore came, with waves 12 feet 
_ high and 15 seconds period striking the piles. An inspection of the =: 
_ showed no detrimental effect. _ The Resident- -Engineer thus questioned i 

_ accuracy of the formulae and slowed down the work. _ A check on the formulae 

_ for this case proved that these waves could not cause over-stressing of the 

piles, because the bending- -moments were proportional to the square of old ae 

— next storm actually brought 15 feet waves, while only 4 out of the 320 =: 
| piles were left. All four piles were bent at the sea- -bed level beyond repair. — aa 
a The above full scale experiment may add evidence to the applicability of aie 
é fed On the other hand, it would be desirable to complete the article by intro- 
ducing data relating t to horizontal elements nig 
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A useful formula for the forces exerted by waves on an immersed cylinder fi 


‘Thus, the resistance to be taken the cylinder is: 


; 


(d) Horizontal elements of structures etc. 


centre of cylinder below mean water level 


-e. the cylinder is not t close to the free 


Or, if we the constant 


Cte ai? 


‘This simple formula shows the main difference between the action of the _ 
waves ona plate or a vertical cylinder to a horizontal cylinder. In the latter 
case the | wave action is shown to increase with the 4th power of the radius and 
inversely with the square of the depth. explains why a small diameter 
; pipe may be laid without any support on the sea- -bottom, while a larger one 4 
suffer more than an equivalent vertical pile. (4) . Also. a submarine vessel 
s. Tt, “bounce” in shallow water while in deep water, since the action 1 decreases: . 
with the ‘Square of the depth, it will be quite steady. 
i Also, a flat plate will be affected by forces equivalent to the exposed area 
(see e.g. formula (8) of the article), while a submerged horizontal ~~ 
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‘Submarine Outfall Pipe; Technion, Israel Institute e of Tech- 


‘Alterman, Within Water: Symposium on Water Conduits; Ministry 
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_ Association of and in Israel. Second e edition 1956. 


Alterman, A Test ona Sea- -Outfall Pile > 


Engineers & Architects in in Israel, Vol. XIV, No. 4 (Nov. Dec. 1956). 


 N. HAMLIN. —The authors’ tests. are a a welcome : addition to the store o 
5, ee available to the designer of offshore ‘structures. . To the writer’s _ 
knowledge these tests were the first to measure wave forces methodically on 
Separate fixed bodies representative of submerged barges in the absence of | 
+ The ‘significant conclusions from the tests and the authors’ analysis thereof 
are first, the predominant wave forces on the barge-like models were inertial 
and ‘second, values of of mass coefficient needed to determine these forces from 
; theoretical calculations compared favorably with values estimated from purely | 
| theoretical considerations or determined experimentally by oscillating models a 
ie In order to estimate e wave forces on a rectangular barge, the ang an 
practice has been to use a mass coefficient based upon data reported by Yu = : 
the authors’ Reference 19. An overall comparison of wave force data from 
model tests of a complete drill rig, together with preliminary data from the 
authors’ tests which were made available to the writer’ s company as one of a 


significantly lower than those from the Yu formula. The comparison = A 
to is reported in the writer’s discussion of a 1957 paper before the 

Z Ss. N.A.M.E. 2 The final horizontal barge forces reported in the present paper 

are considerably higher than the preliminary data indicated and lead to mass © 
coefficients which now agree favorably with the Yu formula. (The final ol 
are about it 15% ‘higher tl than formula values. This tends further confirm a the 


“Yo data as providing means for estimating mass coefficients for | wave an } 


on submerged rectangular barges with acceptable engineering accuracy | 
aa _ However, the results of tests with a flat plate are surprising in that the — cir 
= mass derived from the tests was ‘high compared to theory. In. addition, eee 

there appeared to be a substantial drag force acting on the plate which could — 
only be explained by a drag coefficient higher tl than has been found fromtests _ 
“in steady flow. Regarding the flat plate tests, a possible source of scatter in q 
- the data and of high transient forces might be the tendency of shed ——— } 


; e. This. would seem to be more of a a 
drag effect rather than inertial. Have the authors any evidence in support of 
such an explanation? In this connection, it is of interest that the flat plate . 
1. Bethlehem Steel Co., Shipbuilding Div.. -, Central Technical Dept. ot, Quincy, 


Rechtin, Steele and Scales, P. 633- -681, 1957 S.N.A.M.E. Transactions. — 
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test aus reported by Yu tends . o confirm ‘the theoretical value of CmV wit an 
infinitely long plate \ which, as the authors state, is equal to the volume of a o 
circumscribing cylinder. “Incidentally, with regard to the authors’ reser-_ 
ae concerning Cy, values obtained from high frequency oscillations, the ig 
i | Yute tests in water are reported | for frequencies ranging from about .29 to. = 


c.p. +.p.s., which a1 are lower than the wave frequencies used in the authors’ Pag 4 


toe _ There are presently two interpretations of Eqs. (4) and (5) inuse. Many 

; have worked with wave theory, including the authors, take the symbol = | 
(d+ z) to mean the instantaneous particle height above the bottom. The writer 

of the opinion that (d + z) should mean the average height, in other 

_ words, the height of the center of the particle’ s orbit. . The latter interpretation — f 
enables ‘one to depict orbital motion and other features of the wave in a more > 

i realistic manner. An apparent confusion is indicated, but is of little conse- | 

Fe quence in analysing the authors” tests, inasmuch as both interpretations aa 


equal but opposite upstream and downstream horizontal inertial forces. 
x _ However, assuming (d + z) to represent the average particle height, eine, 


downward for bodies fixed in the wave but near the surface. The tests with _ 

the barge models appear to confirm this tendency, but also show an unexpected 
trend for br upward for forces exceed downward when the models were near the 


The to be congratulated for a clear test 
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4 SNOWMELT RUNOFF* 


> H. C. RIGGS,! A. M. ASCE.— —The authors propose a method of computing 7" 


the rate of runoff from melting snow on the premise that, for practical use, — 
me runoff from melt must be related to generally available ; meteorological data. 
a relation between 


a relation probably would be one derived from data obtained on a one 
plot on which melt water and the meteorological factors were accurately a 
measured. The second part of the problem is that of routing the melt water 


toa stream channel and thence downstream to the gage site, ne a 


is that the reliabilities of the several procedures of on 
meteorological data , distribution of runoff from one day’s melt to streamflow eg 
for | ‘several days, altering of discharge | data on the beginning day of melt 1 to 
affect a trial-and-error solution) cannot be evaluated separately. pape 

_ The authors’ results are incomplete for the purpose of verifying t their pro- 
cedures. Stream | discharge is not related to melt potential alone; the e amount 

of available water in the snow pack and the precipitation during the melting at “a 
period must be included in order to define a usable relation. The wide di- ae - 
vergence of the curves of the Fig. indicates" me need for | additional 


lo 


inches. This equation was an acceptable fit in only 7 of the 15 years tried 


| and therefore, does not constitute verification of the method. No equation 
. - this type, however reliable, could be used for forecasting because results are 


‘given, not in quantitative forms, but asa percentage of cumulative runoff 


_ The most ‘reliable data available should be used in developing and testing a 
method. it is unfortunate that the streamflow records in the selected 
® basin are not of highest accuracy. However, the use of stream discharges re i . 
vised | by the authors does not lend confidence to their results. ‘The “ ice —_ é 
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Hoyt(1) in 1913. Preseht practice is described by Moore. (2) 
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use only when the causative factors can be predicted in advance or when there 
lon an appreciable lag between occurrence of the causative factors and the 
a ing runoff. Compliance with one of these necessary conditions needs to be 


- 


1, U. S. Geological Survey Water-Supply Paper 337. 


—_— Streamflow | under Ice Conditions, ASCE Proc. Paper 1163, pe 
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